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ABSTRACT 



Nucleate pool boiling is an essential part of the vast cooling systems today's 
combatant ship combat systems are dependent upon. Understanding the mechanisms 
that influence heat transfer in tube bundles in a liquid pool is the stepping stone for 
improving these cooling systems. This thesis attempts to bridge the gap between single 
tube performance and bundle performance by studying the effect of a lower heated tube 
on the heat transfer from an upper tube in a simple' two tube bundle. This study con- 
cludes that a nucleating lower tube (regardless of the spacings tested between tubes) has 
a significant positive (i.e. improvement of heat transfer) influence an upon upper tube. 
This is especially evident for a smooth tube where any hysteresis effects are completely 
eliminated when the lower tube nucleates at a heat flux of 10 kW/m 2 or greater. Fur- 
thermore, the only influence for the pitch-to-diameter ratios tested was at the highest 
heat fluxes for the smooth tubes where a p/d of 1.8 was found to give the maximum heat 
transfer. No such maximum was obtained for the enhanced tubes. 
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NOMENCLATURE 



Note that properties and measured parameters of the lower (auxiliary) tube are dif- 
ferentiated from the upper tube with A added to that property or parameter (i.e. ACp 
is the specific heat of the lower (auxiliary) tube). Though this may cause confusion with 
area dimensions, this nomenclature is preserved because it is used within programs 
SETUP72 ( appendix C) and DRP72 (appendix D). Note also that like tubes were used 
in upper and lower positions such that spatial dimensions are the same and are not dif- 
ferentiated between tubes. 

A area 

Ab tube outside surface area of active boiling section 

Ac cross sectional area of the tube 

Cp specific heat 

ACp specific heat of lower (auxiliary) tube 

D diameter 

Di tube inside diameter 

Do tube outside diameter 

Dl diameter of the position of the thermocouple 

D2 outer diameter of the boiling tube 

g gravitational acceleration 

h heat transfer coefficient 

Ah heat transfer coefficient of lower tube 

I current 

A I current for lower tube 

Is output voltage of AC current sensor 

AIs output voltage of AC current sensor for lower tube 

k thermal conductivity of liquid 

Ak thermal conductivity of liquid associated with lower tube 

kc thermal conductivity of copper 

L active boiling tube length 

Lu non-boiling tube length 



x 



Nu 

P 

Pr 

APr 

Q 

AQ 

Qf 

AQf 

Qh 

AQh 

q" 

Aq" 

Ra 

T 

Tavg 

ATavg 

Tc 

Tf 

ATf 

Tn 

Tsat 

Two 

ATwo 

V 

AV 

Vs 

AVs 

a 

Aa 

P 

A P 



Nusselt number 

tube outside wall perimeter 

Prandtl number 

Prandtl number associated with lower tube 

heat transfer rate from boiling surface 

heat transfer rate from boiling surface of lower tube 

heat transfer rate through one non-boiling end 

heat transfer rate through one non-boiling end of lower tube 

heat transfer rate from cartridge heater 

heat transfer rate from lower tube cartridge heater 

heat flux 

heat flux from lower tube 
Rayleigh number 
temperature 

average wall temperature at the thermocouple location 

average wall temperature at the thermocouple location for 
the lower tube 

critical temperature 

film temperature 

film temperature of fluid associated with lower tube 
temperature of the thermocouple location 
saturation temperature 

outer wall temperature of the boiling test tube 

outer wall temperature of the lower tube 

voltage across the cartridge heater 

voltage across the cartridge heater for the lower tube 

voltage output by AC-DC true RMS converter 

voltage output by AC-DC true RMS converter for lower tube 

thermal diffusivity 

thermal diffusivity associated with lower tube 

volumetric thermal expansion coefficient 

volumetric thermal expansion coefficient associated with 
lower tube 
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uncertainty in measurement and calibration 
superheat 

superheat associated with lower tube 
dynamic viscosity 

dynamic viscosity associated with lower tube 
kinematic viscosity 

kinematic viscosity associated with lower tube 
density 

density associated with lower tube 



I. INTRODUCTION 



As the United States continues to struggle with the dilemma of increasing require- 
ments for CFC's, while at the same time joining the world in concern over ozone de- 
pletion, the exigency of finding replacements for the Navy's high ozone depletion 
potential fluids is of the highest order. The urgency of this search was made more 
pressing by the reduction of the time allowed (ordered by President Bush) to meet 
Montreal Protocol deadlines to purge use of CFC refrigerants. A particularly damaging 
fluid is R-114, used primarily in centrifugal chilled-water air-conditioning plants onboard 
ships. In order to best determine a short term or "drop in" replacement, data of R-114 
nucleate pool boiling characteristics is required for comparison. The literature provides 
many studies of single tubes in nucleate pool boiling; however, studies of multiple tubes 
are more scarce. Most studies of multiple tubes with varied tube spacings were con- 
ducted in a manner to simulate a bundle, i.e. all tubes in the pool were studied with the 
same applied heat flux. There is little work in the literature on the effect of a lower 
heated tube on the heat transfer from an upper tube (i.e. a simple bundle) and none with 
R-114. This type of data would begin to bridge the gap between single tube nucleate 
pool boiling and bundle effects, and would complement the search for a suitable re- 
placement for R-114. 

To increase the available data for comparison with possible replacement refrigerants, 
and to further investigate the effect of tube spacing, the following objectives of this thesis 
were established: 

1. Modify the existing single tube pool boiling apparatus to accommodate a simple 
two tube bundle, including program modifications to facilitate instrumenting the 
second tube. 

2. Operate the apparatus to prove repeatability with single tube data. 

3. Obtain convection and boiling data over several tube pitches and heat flux settings 
for both a smooth tube and an enhanced tube. 

Of particular note is that this study was designed as a follow up to the single tube 
work of Sugiyama [Ref. 1]. The format of this study, including the correlations and 
programs used, are adaptations from Sugiyama's work in order to make the results di- 
rectly comparable. 
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II. MECHANISMS 



A. SINGLE TUBE BEHAVIOR 

Experimental heat transfer behavior has been fairly well predicted from single 
smooth cylindrical tubes in an infinite pool for the natural convection region by 
Churchill and Chu (Ref. 2] and Churhchill and Usagi [Ref. 3] and, to a lessor degree, for 
the boiling region by Rohsenow [Ref. 4] and Stephan and Abdelsalam [Ref. 5], Single 
enhanced tubes have also been widely studied with nucleate pool boiling enhancements 
ranging up to 15 times the performance of smooth tubes. The largest enhancements 
have been obtained for re-entrant cavity surfaces as reported by Yilmaz and Westwater 
[Ref. 6], Marto and Lepere [Ref. 7], and Wanniarachchi et u/.[Refs. 8, 9]. In these 
studies, enhancements were largely attributable to increased surface area in the con- 
vection regime, and to stable vapor sites (i.e. vapor trapped in re-entrant cavities), which 
provided a high density of active nucleation sites at relatively low values of wall super- 
heat in the boiling regime. However, to date, no comprehensive model or correlation 
has been presented which can predict such enhancements (Thome addresses this in dis- 
cussion of nucleate pool boiling correlations [Ref. 10]). 

B. MULTIPLE TUBE BEHAVIOR 

The behavior of a particular tube in a multiple tube environment will be greatly in- 
fluenced by its neighboring tubes. In both the natural convection and nucleate boiling 
regimes, fluid heated by a lower tube will rise due to buoyant forces, and an upper tube 
will tend to be affected, specifically by the heated liquid plume in the convection regime 
and by the bubble plume in the boiling regime. Also in a bundle the lower tubes may 
impart a 'drag' on the fluid as it traverses the bundle affecting the fluid reaching an up- 
per tube. The lower tube, because it is not impacted by a similar plume, can be assumed 
to behave as a single tube in a pool. However, if a bundle is in a confined pool, it may 
be the case that strong recirculation patterns cause an upper tube, in turn, to affect the 
heat transfer behavior from a lower tube. 

The effect of a lower tube on an upper tube has been studied by Sparrow and 
Niethammer [Ref. II] and Marsters [Ref. 12] in the convection region (using air), and 
by Fujita et al. [Ref. 13] and by Hahne, Qiu-Rong and Windisch (Ref. 14] in the 
nucleate boiling region (using refrigerants). These effects arc primarily the result of two 
contradictory influences, increased fluid velocity in the convection or boiling bubble 
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plume and increased fluid temperature. In the natural convection region, fluid heated 
by the lower tube will have a velocity due to buoyant forces when it arrives at the upper 
tube. In effect the upper tube is no longer in a true natural convection regime but is 
beginning to experience some forced convection (i.e. mixed convection), thus increasing 
its heat transfer coefficient. At the same time the heat from the lower tube that causes 
this buoyant plume to rise has increased the temperature of that fluid within the plume, 
thereby decreasing the temperature difference between the upper tube and the fluid and 
decreasing the heat transfer capability. Two of the parameters affecting the strength of 
these influences are tube separation and the heat flux setting of the lower tube. The 
temperature of the fluid within the buoyant plume arriving at the top tube will be slightly 
lower than its temperature when it left the lower tube due to convection to the sur- 
rounding fluid. This drop in temperature will tend to increase with increasing tube sep- 
aration, thereby increasing the AT between the upper tube and plume. This implies that 
the upper tube's heat transfer capability at large pitch may be expected to improve over 
that at small pitch. In addition, with Sugiyama's [Ref. 1] single tube work, convection 
plumes were turbulent (Grashof numbers were greater than 10’), implying that plume 
velocity will not change with tube separation. However in reality at high enough 
Grashof numbers, the formation of eddies and other currents in the pool would tend to 
cause a decrease in the plume fluid velocity at the upper tube and thereby a decrease in 
the heat transfer coefficient over that of a smaller pitch. These two effects, increasing 
AT with increasing pitch and decreasing heat transfer coefficient at large pitch suggest 
that an optimum pitch is possible. These influences are confirmed by Sparrow and 
Niethammer's work with air [Ref. 11] and by Marsters [Ref. 12] with air, for which they 
found a dependence on the x-based Grashof number (based on the position of the tubes) 
for the heat transfer from the upper tube. 

In the nucleate boiling regime, the above two effects are also important. The fluid 
entrained in the bubble plume from the lower tube (and the bubbles themselves) impinge 
upon the upper tube at a much higher velocity than with a convection plume (due to 
greater buoyant forces associated with the bubbles). This greater velocity of the plume 
and the bubbles would have the same influence as in convection above and may be suf- 
ficient to strip the thermal boundary layer from the upper tube, thereby lowering heat 
transfer resistance and increasing the heat transfer coefficient. Enhancements due to the 
impingement of plume and bubbles is most significant at lower heat fluxes when active 
nucleation site density on the upper tube is low. This provides more 'contact' area be- 
tween the upper tube and the rising fluid/ bubble mixture. As flux is increased on the 
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upper tube, the active nucleation site density also increases. At high heat fluxes, boiling 
is so vigorous that the rising fluid is prevented from gaining any appreciable 'contact' 
with the upper tube. Heat transfer coefficients may then be expected to increase with 
increasing tube spacing (as in the convection regime) with a possible optimum. How- 
ever, enhancement due to the lower tube may be expected to decrease when high heat 
fluxes are applied to the upper tube. The trends of these influences were verified by 
Fujita et al. [Ref. 13 J and Hahne et al. [Ref. 14J. 

An additional influence in the boiling regime may be the effect of 'seeding' from the 
lower tube when it is boiling. At low and medium heat fluxes on the upper tube (active 
nucleation site density is low), vapor bubbles from the lower tube may impinge and in- 
fluence enter an inactive cavity on the bottom half of the upper tube. This may prompt 
an otherwise inactive site to become active earlier than it would have done for a single 
tube. This influence would also be more effective at low fluxes (low site density) on the 
upper tube. At higher fluxes, the number of active sites is higher and the probability of 
an added site becoming active due to this seeding process is significantly lower. 
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III. DESCRIPTION OF EXPERIMENTAL APPARATUS 



A. GENERAL DESCRIPTION 

The apparatus used is an adaptation of the apparatus used by Karasabun [Ref. 14: 
pp 24-41], Reilly [Ref. 15 : pp. 30-44], and Sugiyama [Ref. 1: pp. 8-32] in previous studies 
of pool boiling of R- 114/oil mixtures from single tubes. Karasabun provides complete 
details of the original configuration of the apparatus. The current apparatus was es- 
sentially that used by Sugiyama for single tubes with modifications made to accommo- 
date two tubes in the evaporator. Additionally, the original R-12 cooling system was 
operated in conjunction with a R-502 system to increase cooling capacity for the in- 
creased heat load of two tube operation. An additional smooth and high flux tube were 
manufactured in order that both a smooth and an enhanced tube set could be studied 
in the two tube configuration. The unused auxiliary variac power supply was employed 
as the power supply for the second tube. A two tube data reduction program DRP72 
was developed by adapting Sugiyama's single tube data reduction program DRP7. 

The apparatus is labeled on a general schematic in Figure 1 and consists essentially 
of seven components: 

1. An evaporator, a boiling vessel assembled using a Pyrex-glass tee with 
aluminum/teflon endplates. 

2. A condenser, assembled using a similar Pyrex-glass tee with aluminum endplates. 

3. A reservoir for R-114 liquid storage. 

4. A cooling subsystem composed of an 1/2 ton R-502 and a 1/4 ton R-12 refriger- 
ation plant, a 30 gallon water/ ethylene glycol sump, and two positive displacement 
pumps. 

5. A vacuum pump. 

6. A data acquisition and instrumentation system. 

7. An aluminum framework with plexiglass siding within which components 1-3 were 
housed. 

The apparatus was designed for reflux operation. Vapor generated in the evaporator 
followed a path through an aluminum L-shaped pipe to the condenser. R-114 
condensate then returned to the evaporator via copper tubing by gravity. The 
watcr/ethylene glycol sump of the cooling subsystem was maintained between -10 and 
-17 °C. This was accomplished using two separate cooling schemes. The R-502 refrig- 
eration plant cooled the mixture via- a counter current heat exchanger tlirough which the 
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Figure 1. Schematic of Apparatus 
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water/ethylene glycol mixture was circulated by an 8-gpm turbine-type pump. The R-12 
refrigeration plant cooled via an evaporator, constructed of coiled copper tubing, located 
directly within the sump. The second 8-gpm turbine-type pump circulated the mixture 
from the sump through copper condensing coils within the condenser via a control valve 
(VC, see Figure 1), and then returned to the sump. 

An alternate mode of operation was used to evacuate the boiler to facilitate tube 
change out. R-114 was boiled off from the evaporator and condensed as before, but in 
this case R-114 condensate was directed to the reservoir (by changing certain valves), 
instead of returning to the evaporator. Once tube changeout was completed the 
evaporator was simply refilled from the reservoir by gravity. 

This study was devoted to pure R-114. In order to ensure that no oil could intrude 
into the system, the oil path between the oil reservoir and the evaporator (used for 
R-l 14/oil mixture tests) was removed and the access to the evaporator plugged. The 
major components of the oil subsystem (reservoir and graduated cylinder) were left in- 
tact for use in future studies. 

B. BOILING TEST SECTION 
1. Evaporator 

The evaporator consisted of a Coming Pyrex-glass tee (3 in. nominal interior 
diameter) with aluminum endplates. The endplates were coupled to the glass tee by 
Coming cast iron removable flanges. The glass vessel was mounted horizontally with 
the side-arm of the tee vertical. The major modification in adapting the apparatus from 
a single tube to two tubes was to the evaporator endplates. New endplates were manu- 
factured with thermocouple and liquid flll/retum penetrations moved to the periphery 
and a three inch hole centered on the axis of the evaporator. Teflon inserts were- man- 
ufactured to fit tightly into this hole with varied spacings for the two test tubes. Though 
the teflon endplates could be rotated to accommodate a variety of tube geometries, only 
vertical spacing was addressed in this study ( i.e. one tube vertically below the other). 
The addition of another tube to the evaporator also required the liquid level to be raised. 
The new liquid level was only 10 mm from the top of the Pyrcx glass tee and 10 mm 
above the top of the upper tube (note that for all spacings the upper tube was fixed in 
position relative to the evaporator and the lower tube was varied). This severely reduced 
the vapor path between the liquid surface and the top of the evaporator to the exit tee. 
A comparison of the modified and previous apparatus liquid levels, tube positions and 
fittings is shown in Figure 2. 
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Comparison of Modified and Previous Apparatus 
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Pyrex-glass was used for both the evaporator and condenser due to its greater 
strength compared to ordinary glass. It was necessary for the apparatus to withstand 
pressure differences (between apparatus interior and atmospheric)of up to 20 psig due 
to possible R-114 vapor pressures at high room temperatures during summer months. 
In addition, a glass vessel provided several other advantages over a metal one: 

• The transparent vessel allowed easy visual observation and videotaping of the 
boiling phenomena occurring inside. 

• The smooth interior surface minimized any nucleate boiling at this inner surface. 



A relief valve (set at 20 psig) was mounted in the aluminum piping between the 
evaporator and condenser to prevent the safe working pressure (30 psig) of the Pyrcx 
glass tee from being exceeded. A sketch of the glass tee and the cast iron flanges is 
shown in Figure 3. A sketch of an endplate and a teflon block insert is shown in 
Figure 4. 

2. Test Tubes 

A general schematic of the tubes is shown in Figure 5. The tubes protruded 
through (and were supported by) the teflon inserts in the endplates of the evaporator. 
The teflon inserts were sealed by viton O-rings, two O-rings between the teflon inserts 
and each tube, and one O-ring between the insert and the aluminum endplate. 

Two types of tube were used, a smooth hard-copper tube and a High Flux 
porous-coated tube. The tubes were heated via a stainless steel, 240-volt, 1000-watt 
(nominal), Watlow Firerod cartridge heater, 6.35 mm in outer diameter, 203.2 mm in 
actual length with an actual heated length of 190 mm. The heater was tightly inserted 
into a copper sleeve, which functioned as a mounting device for the thermocouples. The 
sleeve was then inserted in the tube. Both the heater and the sleeve were tinned with soft 
solder prior to being inserted into the tube and the whole assembly heated in a furnace 
to bond the assembly together. This minimized thermal resistance and provided a uni- 
form heat flux along the length of the tube. The heater and sleeve were inserted such 
that the central 190 mm were heated radially by the heater. This was taken as the active 
boiling length. A correction was applied for the heat lost from the two ends of the tube 
in the pool. 

The tubes were each instrumented with eight thermocouples soldered into 
channels which were machined in the outer surface of the sleeve. These thermocouples 
were located at various axial and circumferential locations as shown in Figure 6. The 
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a) Corning Pyrex Glass Evaporator (D x Ei = 
L = 178 mm, Li = 127 mm) 
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Figure 3. Sketch of Pyrex Glass Vessels and Flanges 
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Figure 4. Sketch of Endplates and Teflon Inserts 
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thermocouple grooves were axially machined to the nearest end of the sleeve to provide 
a path for the thermocouple leads. Type-T Teflon coated copper-constantan 
thermocouple wire was used both in the tubes and throughout the apparatus. Dimen- 
sions and for the tubes are given in Table 1. 



Table 1. TEST TUBE DIMENSIONS 



Tube Type 


Dl(mm) 


D2(mm) 


Di(mm) 


Do(mm) 


L(mm) 


Lu(mm) 


Kwall 

{Wjm.K) 


Smooth 


12.44 


15.88 


12.70 


15.88 


190.0 


76.20 


344 


High Flux 
(95/5 Copper- 
Nickel) 


12.95 


15.82 


13.20 


15.82 


190.0 


76.20 


45 



C. CONDENSER SECTION 

The condenser was assembled similarly to the evaporator. The same size Pyrex-glass 
tee was used as the main vessel with aluminum endplates. The glass vessel was mounted 
vertically with side-arm tee horizontal to receive the R-114 vapor. A helical condensing 
coil made of 9.5 mm copper tubing was inserted in the Pyrex glass tee. Swagelock fit- 
tings were used to connect the condenser coil to the coolant tubing through the 
condenser endplates. The coil was fabricated to an approximate outside diameter of 76 
mm, providing an active condensation length of approximately 4.5 m. 

The vacuum pump was connected to the top of the condenser via two isolation 
valves to remove any noncondensable gases that collect there. A tubing connection was 
placed in the bottom endplate of the condenser to enable R-114 condensate to drain 
back to the evaporator by gravity. Valve and tubing lineup was such that condensate 
could also be directed back to the R-114 reservoir. The condenser was connected to the 
evaporator by an L shaped aluminum tube two inches in diameter. The pressure relief 
valve, mentioned earlier, and a bourdon pressure gage were mounted on this aluminum 
tube. 

D. COOLING SECTION 

1. Coolant Sump 

A 0.15 m \ rectangular sump for the water-ethylene glycol mixture was made of 
13 mm Plexiglas sheet. The sides were glue jointed using a methylene chloride solution 
with additional support provided by small screws. The tank was placed on a wooden 
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Figure 5. Schematic of the lest Tube 
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Positions of the Thermocouples 
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SECTION A-A (DIMENSIONS IN MM ) 



platform to separate it from the concrete floor and all sides were insulated with 22 mm 
sheet insulation. The coolant was a 52% mixture of ethylene glycol/distilled water 
producing a solution freezing point of approximately -25 °C. 

2. Refrigeration Plants 

A 1/2 ton R-502 and a 1/4 ton R-12 refrigeration plant were used to chill the 
ethylene glycol/water mixture. Each consisted of a hermetically sealed compressor as- 
sembly, an air cooled condenser, a receiver, a filter-dryer, a pressure regulator, a tem- 
perature control switch and a thermostatic expansion valve. A counter current heat 
exchanger was used as the evaporator for the R-502 plant. An evaporator for the R-12 
plant was constructed of 9.5 mm copper tubing which was coiled and placed directly in 
the coolant sump. Both plants were controlled by a temperature control switch with 
each plant having its own thermostatic expansion valve. The plants were adjusted to 
maintain coolant temperatures at about -17 °C. Figure 7 shows a schematic of the 
R-502 refrigeration plant. Figure 1 shows placement of the R-12 refrigeration plant 
evaporator. 

3. Pump and Control Valve 

The two 8 gpm, 115 V Burks turbine type, positive-displacement pumps were 
floor mounted next to the coolant sump. One pumped the ethylene glycol/water coolant 
mixture through the R-502 counter-current heat exchanger. The other pumped the 
coolant mixture to the condensing coil via control valve VC (see Figure 1). Valve VC 
was the operators primary control device to maintain saturation pressure at a desired 
value by varying the coolant flow rate (and hence the condensation rate) in the 
condenser. A bypass line was also installed on the discharge line from the pump. Dis- 
charge bypass was controlled by valve V9. The bypass line served to avoid overloading 
the positive displacement pump when small heat loads were placed on the condenser 
(low coolant flow rates). The bypass return line returned to the sump via a penetration 
in the top of the tank, where the coolant return stream from the condenser returned as 
well. 

E. R-l 14 RESERVOIR 

The R-l 14 reservoir was an aluminum cylindrical vessel, 230 mm in diameter and 
254 mm in height. The reservoir was equipped with a transparent plastic tubular sight 
glass to monitor liquid level. The reservoir was located within the apparatus enclosure 
at a vertical location above that of the evaporator and below that of the condenser to 
facilitate both gravity flow from the condenser to the reservoir and from the reservoir' 
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Figure 7. Schematic of R-502 Refrigeration System 
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Air-Coo Jed Condenser 



to the evaporator. The reservoir could be replenished from an outside source (storage 
cylinder) via a connection on the vacuum line from the condenser. Vapor from the 
storage cylinder was condensed in the condenser and drained into the reservoir. The 
arrangement of the reservoir is shown in Figure 1. 

F. APPARATUS ENCLOSURE 

The apparatus (apart from the cooling subsystem) was contained in an aluminum 
framed, plexiglas enclosure. The enclosure consisted of a rectangular welded frame ap- 
proximately 1 m x 0.5 m x 0.6 m. The frame was divided into upper and lower halves 
with the lower half of the frame straddling and enclosing the coolant sump. The upper 
half of the frame was enclosed at the top and bottom with aluminum sheet. The four 
remaining sides were enclosed with 13 mm Plexiglas. The opposing 0.5 m x 0.6 m 
plexiglas sides were equipped with hinges to allow access to the interior. The refriger- 
ation plants, pumps, and heat exchanger were located on the laboratory floor adjacent 
to the enclosure. Data acquisition equipment was located in a separate cabinet also 
adjacent to the enclosure. The enclosure provided some insulation from ambient con- 
ditions and provided a safety barrier between personnel and the apparatus should one 
of the glass vessels fail. 

G. INSTRUMENTATION 

1. Power Measurement 

A 240 volt AC supply was used as the power source for the apparatus and was 
controlled via a variac. Output from the variac ranged from 0-220 V, 0-5 A, adjustable 
by the operator, to obtain a desired heat flux at the tube surface. Power input to the 
tube heater was measured with an AC current sensor and a voltage sensor (both sensors 
output in volts). The voltage sensor output was also processed through an AC-DC true 
R.M.S. converter which provided a proportional output signal (also in volts). The cur- 
rent and R.M.S. voltage output were input to the data acquisition system. Calibration 
of power measurement was checked by comparing voltage and amperage measurements 
against a digital voltmeter and ammeter. The power measurement system is shown 
schematically in Figure 8. 
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Figure 8. Schematic of the Power Measurement 
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2. Temperature Measurement 

Several temperatures were monitored throughout the apparatus. These were 

• Average test tube wall temperature (using the 8 thermocouples positioned as men- 
tioned above). 

• Pool temperature (three thermocouples located at three different positions 
longitudinally within the pool but at approximately the same height) 

• Coolant sump temperature (one thermocouple) 

The pool thermocouples were inserted into special housings that penetrated the 
evaporator endplates (Figure 9). The main body of the housing was manufactured from 
stainless steel (low thermal conductivity) to minimize errors in pool temperature meas- 
urement due to conduction from the surroundings through the housing. The tip of the 
housing was manufactured from copper to take advantage of copper's high thermal 
conductivity and to minimize the temperature difference between the pool and the 
thermocouple. 

Copper-constantan thermocouples were used for all measurements. 
Thermocouples were read directly by a Hewlett-Packard 3497A data acquisition system 
controlled by a Hewlett-Packard 9826 computer. The average temperature for each 
thermocouple was obtained by scanning its output 20 times over 5 seconds and taking 
an average. 

H. DATA ACQUISITION AND REDUCTION 

All sensor outputs (thermocouples, current sensor and R.M.S. voltage) were ana- 
lyzed by a Hewlett-Packard 9826 computer and the data stored using the iterative data 
collection/reduction program DRP72 (Sugiyama's single tube program DRP71 [Ref. 1] 
was used during initially repeatability experiments). The program was controlled (by the 
operator) by keyboard interaction to prompt the system to take desired steps. Raw data 
was stored on the computer hard drive while a printout of reduced data was provided 
on a Hewlett-Packard Inkjet printer. Channel assignments for the various sensor inputs 
to the data acquisition system are shown in Table 2. 



19 




Figure 9. Sketch of Thermocouple Housing 
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Table 2. HP 3497A CHANNEL ASSIGNMENTS 



Channel 


Assignment 


00-07 


Upper test tube wall temperatures 


08-15 


Lower test tube wall temperatures 


16-18 


Pool liquid temperatures 


19 


Coolant sump temperature 


20 


RMS voltage tube heaters 


21 


Upper tube current sensor 


22 


Lower tube current sensor 
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Following data acquisition, the data were reduced utilizing the following 
procedure^ 

1. Input the name of the user- specified file to be stored on hard drive. 

2. Select number of tubes to be powered and from which tube data is to be collected 
(upper or lower). 

3. Select tube type (all dimensions of the boiling test tubes are then automatically 
selected). 

4. Set desired saturation temperature (°C) of the pool (for these tests, 2.2 °C was used, 
corresponding to approximately 1 atmosphere). 

5. Input desired heat flux setting. 

6. Set desired heat fluxes by adjusting appropriate variac. 

7. Set desired saturation temperature by adjusting flow of coolant through condenser 
coils with control valve VC. 

8. Once saturation temperature is achieved, wait for steady state conditions (minimum 
of 5 minutes) prior to taking data. 

9. Prompt data acquisition unit to scan all channels listed in Table 2. All channel 
readings are made in volts and stored in user specified fields. 

10. Compute parameters from these voltages (i.e. temperature and power). 

1 1. Compute the heat transfer rate from the cartridge heater. 

12. Compute the av erage wall temperature of the test tubes and calculate the wall 
superheat (Twall — Tsat) 

13. Compute the physical properties of R-l 14 using the property correlations (see 
Sugiyama[Ref. 1]) at film temperature ((Twall — Tsat)/2). 

14. Compute the natural convection heat-transfer coefficient of R-l 14 from the un- 
heated ends of the test tubes. 

15. Compute heat losses from the unheated ends of each tube. 

16. Calculate heat flux from the heated length of each tube. 

17. Calculate the heat-transfer coefficient of the R-l 14 from the heated length 
of each tube. 

18. Store heat flux and wall superheat for each data set in user specified field. 

19. Plot the data using available software. For this study, plots were produced on the 
NPS mainframe computer. 

Sample calculations of the above procedure are given in Appendix A. 



1 This procedure was essentially identical to that followed by Sugiyama (Ref. 1| and therefore 
Ills procedure is largely reproduced. 
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IV. EXPERIMENTAL PROCEDURE 



A. ASSEMBLY 

The mechanical assembly of the tubes within the evaporator was very straight for- 
ward and will not be discussed in great detail. However the preparatory steps following 
assembly are important from the point of view of ensuring a leaktight apparatus to 
prevent loss of R-114 from the system and the intrusion of noncondensable gases into 
the system. Additionally it should be noted that vertical alignment of tubes was ensured 
by using an air bubble level. 

1. Pressure Test of the Apparatus 

Following reassembly of the evaporator, a pressure test was conducted on the 
apparatus by pressurizing the system (with air) up to 10 psig. A soap-water mixture 
then was sprayed on all joints and leaks were detected by bubble formation from the 
soap-water mixture. Any leaks were then repaired. The apparatus was repeatedly 
pressure tested until there were no "visible” leaks detected. 

2. Vacuum Test 

Following a successful pressure test, the apparatus was evacuated (using the 
vacuum pump) to approximately 27 inHg. The apparatus was then allowed to stand 
overnight. If an appreciable vacuum loss was observed (i.e. a drop of more than l inHg. 
in a 24 hour period), the whole process was repeated, beginning with a pressure test, 
until all leaks were corrected. 

3. Filling the Evaporator 

After the apparatus was successfully vacuum tested, the evaporator was filled 
from the R-l 14 reservoir. The pressures were first equalized between the evaporator and 
the reservoir by cracking open valve V7. R-114 liquid was then drained by gravity into 
the evaporator through valve V6 until the liquid reached the desired evaporator level, 
approximately 10 mm above the top of the upper tube. Once the. evaporator was filled, 
the reservoir was isolated from the rest of the apparatus by closing all connecting valves. 

4. Sensor Connections 

Thermocouple leads and any other sensor connections were not reconnected 
until the evaporator was successfully filled without leaks. This prevented unnecessary 
wear and tear on the heater leads and thermocouple wire due to the large number of tube 
changes required. 
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5. Degassing and Acquisition Channel Check 

Once all connections were made, program SETUP72 checked the output on all 
channels. Any erroneous output was promptly investigated and the necessary cor- 
rections made. Before any heat was applied to either test tube, the system pressure was 
reduced to the desired saturation pressure by opening flow of the coolant mixture to the 
condenser. Each tube in turn was boiled vigorously at the highest heat flux setting 
(approximately 80 kW/m 2 ) for about 10 minutes in order to degas the refrigerant and to 
remove excess air from the tube surface cavities. Noncondensable gases then collected 
at the top of the condenser and were removed by the vacuum pump. During the 
degassing procedure, program SETUP72 was again used to observe thermocouple 
readings and power sensors to ensure proper operation. Any faulty thermocouples were 
marked to ensure they were excluded during data runs. Following degassing, the appa- 
ratus was secured and the tubes allowed to soak overnight in the R-114 pool to provide 
good surface wetting. 

B. NORMAL PROCEDURE 

The following procedure was used to obtain heat transfer data2: 

1. Prior to system startup, the refrigeration plants were operated for approximately 1 
hour to reduce the coolant sump temperature to about -15 °C. 

2. The data acquisition/control unit, computer and variac panel were switched on. 

3. The computer program SETUP72 was loaded and run: 1) All data acquisition 
channels were then rechecked. 2) A power output of approximately 2.5 Watts was 
input into each tube to check operation of the heater and power sensors. 3) The 
average temperature of the refrigerant was slowly reduced to 2.2 °C by circulating 
a small amount of coolant through the condenser. 

4. The data acquisition program DRP72 was loaded and run. 

5. The desired heat flux setting for the lower tube was input to the program and the 
lower tube variac adjusted to give the desired heat flux (within ± 500 IV/m 2 ). 

6. The control valve VC was adjusted to regulate the flow of cooling liquid through 
the condenser to maintain a constant saturation temperature at a given heat flux. 
Desired versus actual saturation temperatures were monitored continuously by the 
program until they agreed to within ± 0.1 °C. 

7. For each data point, conditions in the evaporator were allowed to stabilize for at 
least five minutes prior to taking data. The following raw data were measured and 
stored in a user specified field: local test tube wall temperatures, pool liquid tem- 
peratures, sump temperature, current sensor readings and voltage sensor readings. 



2 Ibid 
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8. Two data points were taken at each desired heat flux and saturation temperature. 
The following processed data were recorded as a printout: wall temperatures of the 
test tubes, liquid pool temperatures, vapor temperature, sump temperature, wall 
superheat, heat-transfer coefficient and the heat flux. 

9. For each data set, the above procedure was repeated from step 5. Various heat flux 
steps were used to obtain uniform steps on a log-log scale. For increasing flux, 
smaller steps were used up until the incipience of boiling in order to accurately de- 
termine the point of incipience. 

C. DATA LISTING 

Data were filed using the following file name system: 

• Files were given 10 character names (ex. DAT0531I52). 

• The first three characters DAT simply refer to data. 

• The next four characters date the file (0531 = 31 May). 

• The eighth character indicated an increasing or decreasing run (I = increasing, D 
= decreasing). 

• The ninth character indicated tube type ad listed in program DRP72 (4 = smooth 
tube, 5 = High Flux tube). 

• The tenth character indicated which program wad used DRP71 or DRP72 (1 = 
DRP71, 2 - DRP72). 

• The 1 1th character when present indicated successive runs for the indicated tube 
and indicated program for that day. 

Table 3 is a listing of data runs. 
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Table 3. LISTING OF DATA RUNS 



Data File 


Tube Type 


Purpose 


DAT0122I51 


High Flux 


Cal; repeatability; single tube 


DAT0123I52 


High Flux 


Cal; repeatability; single tube 


DATO 125152 


High Flux 


Repeatability; p/d =2 


DATO 12712 


High Flux 


Repeatability; p/d =2 


DATO 130 152 


High Flux 


Repeatability; p/d =2 


DAT0130D52 


High Flux 


Repeatability; p/d = 2 


DAT03 12142 


Smooth 


Data; repeatability; p/d= 2 


DAT0312D41 


Smooth 


Data; repeatability; p/d = 2 


DAT03 13141 


Smooth 


Repeatability; p/d = 2 


DAT0313D41 


Smooth 


Repeatability; p/d =2 


DAT03 17142 


Smooth 


Repeatability; p/d= 2 


DAT0317D42 


Smooth 


Repeatability; p/d =2 


DAT0320I41 


Smooth 


Repeatability; p/d = 2 


DAT0320D41 


Smooth 


Repeatability; p/d =2 


DAT0330I42 


Smooth 


Repeatability; p/d =2 


DAT0330D42 


Smooth 


Repeatability; p/d =2 


DAT0331I42 


Smooth 


Data; p/d= 2; 500 W/m 2 on lower 
tube 


DAT0331D42 


Smooth 


Data; p/d=2; 500 IV/m 2 on lower 
tube 


DAT0401I42 


Smooth 


Data; p/d =2; 1 klV/m 2 on lower 
tube 


DAT0401D42 


Smooth 


Data; p/d = 2; 1 k IV/m 2 on lower 
tube 


DAT0404I42 


Smooth 


Data; p/d = 2; 1 klVIm 2 on lower 
tube 


DAT0404I421 


Smooth 


Data; p/d = 2; 1 klV/m 2 on lower 
tube 


DAT0404D421 


Smooth 


Data; p/d=2; l klV/m 2 on lower 
tube 


DAT0406I42 


Smooth 


Data; p/d = 2; 10 klV/m 2 on lower 
tube 
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Continuation of Table 3. 



Data File 


Tube Type 


Purpose 


DAT0406D42 


Smooth 


Data; p/d = 2; 10 klV/m 2 on lower 
tube 


DAT0407I42 


Smooth 


Data; p/d =2; 25 klV/tn 2 on lower 
tube 


DAT0407D42 


Smooth 


Data; p/d= 2; 25 klV/m 2 on lower 
tube 


DAT0408I42 


Smooth 


Data; p/d = 2; 3 klVjm 2 on lower 
tube 


DAT040SD42 


Smooth 


Data; p/d = 2; 3 klV/m 2 on lower 
tube 


DAT04 10142 


Smooth 


Data; p/d = 2; Tubes run up to- 
gether 


DAT0410D42 


Smooth 


Data; p/d = 2; Tubes run down 
together 


DAT0415I42 


Smooth 


Data; p/d= 1.8; 10 kWInt 2 on 
lower tube 


DAT0415D42 


Smooth 


Data; p/d= 1.8; 10 kWIm 2 on 
lower tube 


DAT04 151421 


Smooth 


Data; p/d= 1.8; 25 kW/m 2 on 
lower tube 


DAT0415D421 


Smooth 


Data; p/d= 1/8; 25 klVIm 2 on 
lower tube 


DAT04 17142 


Smooth 


Data; p/d = 1.8; no power on 
lower tube 


DAT0417D42 


Smooth 


Data; p/d= 1.8; no power on 
lower tube 


DAT041S142 


Smooth 


Data; p/d= 1.8; 1 klV/m 2 on 
lower tube 


DAT0418D42 


Smooth 


Data; p/d* 1.8; 1 klV/m 2 on 
lower tube 


DAT0419I42 


Smooth 


Data; p/d* 1.8; 3 kWfm 2 on 
lower tube 


DAT0419D42 


Smooth 


Data; p, d= 1.8; 3 klV/in 2 on 
lower tube 


DAT0423I42 


Smooth 


Data; p/d = 1.5; no power on 
lower tube 
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Continuation of Table 3 



Data File 


Tube Type 


Purpose 


DAT0423D42 


Smooth 


Data; p/d= 1.5; no power on 
lower tube 


DAT0424I42 


Smooth 


Data; p/d= 1.5; 3 kW/m 2 on 
lower tube 


DAT0424D42 


Smooth 


Data; p/d= 1.5; 3 kW/m 2 on 
lower tube 


DAT0425I42 


Smooth 


Data; p/d= 1.5; 1 klV/m 2 on 
lower tube 


DAT0425D42 


Smooth 


Data; p/d= 1.5; l kW/m 2 on 
lower tube 


DAT0426I42 


Smooth 


Data; p/d= 1.5; 10 kW/m 2 on 
lower tube 


DAT0426D42 


Smooth 


Data; p/d= 1.5; 10 kW/m 2 on 
lower tube 


DAT0427I42 


Smooth 


Data; p/d= 1.5; 25 kW/m 2 on 
lower tube 


DAT0427D42 


Smooth 


Data; p/d= 1.5; 25 kW/m 2 on 
lower tube 


DAT0430I52 


High Flux 


Data; p/d= 1.5; 1 klV/m 2 on 
lower tube 


DAT0430D52 


High Flux 


Data; p/d= 1.5; 1 klV/m 2 on 
lower tube 


DAT0501I52 


High Flux 


Data; p/d= 1.5; 10 klV/m 2 on 
lower tube 


DAT0501D52 


High Flux 


Data; p/d= 1.5; 10 klV/m 2 on 
lower tube 


DAT050 11522 


High Flux 


Data; p/d= 1.5; 25 klV/m 2 on 
lower tube 


DAT0501D522 


High Flux 


Data; p/d= 1.5; 25 kW/m 2 on 
lower tube 


DAT0502I52 


High Flux 


Data; p/d= 1.5; 3 klV/m 2 on 
lower tube 


DAT0502D52 


High Flux 


Data; p/d= 1.5; 3 klV/m 2 on 
lower tube 


DAT0503I52 


High Flux 


Data; p, d= 1.5; no power on 
lower tube 
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Continuation of Table 3. 



Data File 


Tube Type 


Purpose 


DAT0503D52 


High Flux 


Data; p/d= 1.5; no power on 
lower tube 


DAT0506I52 


High Flux 


Data; p/d= 1.8; no power on 
lower tube 


DAT0506D52 


High Flux 


Data; p/d = 1.8; no power on 
lower tube 


DAT0507I52 


High Flux 


Data; p/d=. 1.8; 1 klV/m 2 on 
lower tube 


DAT0507D52 


High Flux 


Data; p/d= 1.8; 1 kW/m 2 on 
lower tube 


DAT0508152 


High Flux 


Data; p/d = 1.8; 3 klVjm 2 on 
lower tube 


DAT0508D52 


High Flux 


Data; p/d= 1.8; 3 kWjm 2 on 
lower tube 


DAT0508I521 


High Flux 


Data; p/d= 1.8; 10 kW/m 2 on 
lower tube 


DAT0508D521 


High Flux 


Data; p/d= 1.8; 10 AlF/m 2 on 
lower tube 


DAT0509I52 


High Flux 


Data; p/d= 1.8; 25 kWjm 2 on 
lower tube 


DAT0509D52 


High Flux 


Data; p/d= 1.8; 25 klVjm 2 on 
lower tube 


DAT0509D52 


High Flux 


Data; p/d= 1.8; 25 kWjm 2 on 
lower tube 


DAT0515I52 


High Flux 


Data; p/d= 2; 25 kW}m 2 on lower 
tube 


DAT0515D52 


High Flux 


Data; p/d= 2; 25 kWjrn 2 on lower 
tube 


DAT05 16152 


High Flux 


Data; p/d= 2; 10 kWjm 2 on lower 
tube 


DAT0516D52 


High Flux 


Data; p/d= 2; 10 kW/m 2 on lower 
tube 


DAT05 17152 


High Flux 


Data; p/d = 2; 1 klV/m 2 on lower 
tube 


DAT0519I52 


High Flux 


Data; p/d= 2; 3 kW/m 2 on lower 
tube 
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Continuation of Table 3 



Data File 


Tube Type 


Purpose 


DAT0520I52 


High Flux 


Data; p/d= 2; no power on lower 
tube 


DAT0529I52 


High Flux 


Data; p/d= 2; tubes run up/down 
together 


DAT0530I52 


High Flux 


Data; p/d= 2; lower tube position 
plugged 
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V. RESULTS AND DISCUSSION 



It should be noted and it is stressed that results are expressed in terms of enhance- 
ment of the upper tube over a single tube and that the results are applicable to this ap- 
paratus. These results may not be reproducible outside of this apparatus. 

A. REPRODUCIBILITY 

To investigate repeatability, tests were conducted with both a 'single' High Flux and 
smooth tube (with only the upper tube operating). Runs were conducted using both 
program DRP71 (from Sugiyama's single tube work [Ref. 1]) and DRP72 with no heat 
flux applied to the bottom tube. This ensured that the new program DRP72 performed 
correctly. An additional set of runs were made using the pair of smooth tubes with a 
heat flux, of 1 klV/m 2 applied to the lower tube using program DRP72. The runs for all 
these comparisons are listed in Table 4. A pitch-to-diameter ratio of 2 was used for the 
runs listed in Table 4. A plot of heat flux vs. wall superheat (Thetab) for these runs is 
shown in Figure 10 through Figure 12. Figure 10 shows the 'single' smooth tube data 
together with the single smooth tube data of Sugiyama [Ref. 1]. The three runs agree 
very well and are within 10% of Sugiyama's data; any discrepancy is attributed to the 
presence of the unheated lower tube modifying the flow over the upper tube. Clear ev- 
idence of hysteresis is seen as with Sugiyama's work; as the heat flux is decreased , a 
closed loop is formed indicating that as the nucleation sites die out, the transfer mech- 
anism returns to one of natural convection. Figure 1 1 shows the 'single' High Flux tube 
data together with the single tube High Flux data of Sugiyama [Ref. 1], Again good 
agreement is obtained except at theTowest heat fluxes and lowest AT's where discrep- 
ancies of up to 50% are noted. However, due to the very high uncertainty in this region 
(marked on Figure 11), primarily due to uncertainty in the wall temperature measure- 
ment due to the fabrication procedures, the agreement is still considered reasonable. It 
should also be noted that run DAT0530 was conducted with plugs in place of the bottom 
tube (as opposed to the bottom tube having zero applied heat flux) and also shows ex- 
cellent agreement. This indicates quite clearly that with the lower tube off, the upper 
tube behaves almost exactly like a single tube. Note also that the hysteresis 'loop' is 
open (again similar to the work of Sugiyama) indicative of the fact that there are still a 
significant number of active nucleation sites operating with the High Flux tube even at 
heat fluxes as low as 1000 IV/m 2 . Figure 12 shows the three smooth tube runs con- 
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ducted with a IklV/m 2 heat flux applied to the bottom tube. Without discussing the 
significance of the results here (see section B.), it can be seen that repeatability is good 
and within 5% for all three data sets. Based on Figure 10 through Figure 12, data re- 
peatability and agreement with past single tube work carried out on the same apparatus 
is considered good. 



Table 4. DATA RUNS USED FOR REPEATABILITY 



Data Set 


Tube Type 


Lower Tube Setting 


DAT0313I41 


smooth 


no power on lower tube 


DAT0313D41 


smooth 


no power on lower tube 


DAT0320I41 


smooth 


no power on lower tube 


DAT0320D41 


smooth 


no power on lower tube 


DAT0330I42 


smooth 


no power on lower tube 


DAT0330D42 


smooth 


no power on lower tube 


DAT0401I42 


smooth 


1 kWjm 2 


DAT0401D42 


smooth 


1 kW/m 2 


DAT0404I42 


smooth 


1 kW/m 2 


DAT0404D42 


smooth 


1 kWIm 2 


DAT04G7I42 


smooth 


1 kWjm 2 


DAT0407D42 


smooth 


1 kW\m 2 


DAT0123I52 


High Flux 


no power on lower tube, increas- 
ing run only 


DATO 130152 


High Flux 


no power on lower tube 


DAT0130D52 


High Flux 


no power on lower tube 


DAT0530I52 


High Flux 


lower tube removed and position 
plugged, increasing and decreas- 
ing run recorded in this file 



B. INFLUENCE OF THE LOWER TUBE 

To investigate the influence of a lower upon an upper tube, each upper tube (for 
both types of tube surface) was run up and down (starting in the convection region) 
through various heat flux settings for three different pitch-to-diameter (p/d) ratios and 
five different fixed heat flux settings on the lower tube. Note that in each test, the same 
tube surface was used for the upper tube and the lower tube (i.e. the surfaces were never 
mixed). The p/d ratios used were 1.5, 1.8 and 2. A p/d ratio of 2 was selected to match 
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Figure 10. Repeatability Comparison for Smooth Tube, No Heat Flux on Lower 
Tube 
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Repeatability Comparison for Hi Flux Tube, No Heat Flux on Lower 
Tube 
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Figure 12. Repeatability Comparison for Smooth Tube, Heat Flux of 1 k IV I in* on 
Loner Tube 



the vertical pitch, not only in Naval centrifugal flooded evaporators, but also in the 
bundle apparatus in the Heat Transfer Laboratory. The fixed lower tube heat flux set- 
tings (LTHFS) used were 0, 1, 3, 10 and 25 kW/m 2 (there was one test conducted with 
a lower tube heat flux setting of 500 W/m 2 ). To evaluate the influence of the lower tube, 
the data has been presented in two ways: the first by common p/d ratio and the second 
by common LTHFS. In every case, plots of heat flux vs AT are shown. 

1. The Effect of Heating the Lower Tube 
a. Smooth Tube p/d Ratio of 2 

Figure 13 shows data taken for a p/d ratio of 2 for the six LTHFS labeled 
in the legend representing klV/m 2 and CC representing the Churchill and Chu correlation 
[Ref. 2). This labeling was used for all the 5 subsequent plot legends. An extra run was 
made at this pitch with a LTHFS of 500 W/m 2 for comparison. The natural convection 
regime shows very small enhancements with increasing LTHFS of up to 3 kW/m 2 . Of 
note is that for LTHFS of 0, 0.5, 1, and 3 kW/m 2 , the lower tube was in the convection 
regime at the beginning of the increasing run. For the 10 and 25 kW/m 2 case the lower 
tube had nucleated immediately and was nucleating for the entire run. The natural 
convection curves are still very close to values calculated from Churchill and Chu's cor- 
relations [Ref. 2] and relatively parallel with one another. It can be seen that for a 
LTHFS of 500 W/m 2 , the top tube behaves almost as if it were a single tube (i.e. agrees 
with a LTHFS of 0). In this case, the lower tube remained in the natural convection 
regime for the entire run and shows that if this is the case, little to no enhancement of 
the upper tube is obtained at any heat flux. The point of incipience is different for each 
LTHFS and proved to be a fairly random event for the smooth tubes, occurring around 
12 kW/m 2 . In all cases except the 0 and 500 W/m 2 , the lower tube was "seeded" by 
bubbles from the R-l 14 return line (impinging upon the lower tube) when at a high flux 
(> 20 kW/m 2 ) on the upper tube. The lower tube then continued to nucleate after being 
'seeded' throughout the remainder of the run. The LTHFS of 0, 0.5, 1 and 3 klV/m 2 
showed enhaacement upon nucleation, jumping to a common nucleate boiling curve. 
This entire run showed a very definite hysteresis effect and was very similar to 
Sugiyama's [Ref. 1] single tube data. Following nucleation, the four lower LTHFS fol- 
lowed an identical nucleate boiling curve for the remainder of the increasing run. On the 
decreasing run, the two lowest LTHFS showed a gradual return to the convection curve 
indicating a gradual deactivation of nucleation sites as heat flux was reduced on the 
upper tube. This coincided with no nucleation on the lower tube and also agreed well 
with Sugiyama's [Ref. 1) single tube results. The remaining LTHFS departed from this 
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curve (on the decreasing run) and showed significant enhancements (almost an order of 
magnitude over a single tube at the lowest heat fluxes). The enhancement increased with 
increasing the LTHFS until at LTHFS above 10 kWjm 2 , additional enhancement was 
minimal. The curves followed by LTHFS greater than 1 kWfm 2 during the decreasing 
run were also parallel to the increasing run natural convection curves. This suggests 
strong convection effects on the upper tube when the lower tube is nucleating and heat 
flux setting on the upper tube is below 10 klVIm 2 . This agrees well with the data of 
Fujita et al. [Ref. 13]. It should be noted that all LTHFS followed the same curve above 
an upper tube heat flux setting of approximately 20 kW/m 2 showing no effect of LTHFS 
at high flux settings on the upper tube. This too is in agreement with Fujita et al.. Of 
particular note is that the two highest LTHFS follow the same curve for both increasing 
and decreasing runs, indicating complete elimination of hysteresis prominent in single 
tube data. The above discussion indicates a significant effect of the lower tube upon the 
upper tube when the lower tube is nucleating. There are almost negligible enhancements 
when both tubes are in the convection regime. 

b. Smooth Tube p\d Ratio of 1.8. 

The smooth tube jjins with a p/d ratio of 1.8 , shown in Figure 14, show 
very similar results to that for a p/d of 2. As before for LTHFS of 0, 1 and 3 kW/m 2 , 
both tubes begin in the convection regime for the increasing run: again no significant 
enhancement is seen in this region, although these curves are in better agreement with 
Churchill and Chu [Ref. 2]. For a given LTHFS, onset of nucleate pool boiling (ONB) 
occurs at different heat fluxes than in Figure 13, demonstrating the random nature of 
this phenomena. A common nucleate boiling curve is followed by all LTHFS runs 
above an upper tube heat flux setting of about 20 kW/m 2 . Hysteresis effects are similar 
for low LTHFS as are the elimination of hysteresis for the two highest LTHFS. The 
lower tube was seeded by bubbles from the R- 114 return line at high upper tube heat 
flux for all LTHFS (except 0). When compared with Figure 13, the same enhancement 
factors are seen at all but the highest heat fluxes, where the 1.8 p/d shows a slight en- 
hancement over the 2 p/d. This will be shown more clearly on a later graph. 

c. Smooth Tube p\d of 1.5 

The same trends and enhancements seen in the previous two p/d ratios are 
also seen in Figure 15. LTHFS of 0, 1 and 3 are very similar to the previous two p/d 
ratios with no appreciable enhancement in the convection region. The lower tube 
started in the convection region for LTHFS of 0, 1 and 3, and was nucleating throughout 
for LTHFS of 10 and 25 k IVjm 2 The lower tube nucleated via 'seeding' at high fluxes 
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Figure 13. Comparison of Lower Tube Flux Settings for a p/d Ratio of 2 for Smooth 
Tube 
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Figure 14. Comparison of Lower Tube Flux Settings for a p/d Ratio of 1.8 for 
Smooth Tube 
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on the upper tube as for the previous two p/d ratios. ONB for the upper tube continued 
to be a statistical phenomenon. Similar hysteresis patterns are also seen. At the highest 
heat flux settings, the 1.5 p/d shows a degradation compared to both the 1.8 and 2.0 p/d. 
This will again be shown more clearly on a later graph. 
d. High Flux Tube p/d of 2 

A plot of data runs made for the High Flux tube with a p/d ratio of 2 is 
shown in Figure 16. The convection region looks very similar to that of the smooth 
tubes in that they agree fairly well with Churchill and Chu [Ref. 2]. Also, there is no 
appreciable enhancement in the convection region for increases in LTHFS from 0 to 3 
kWjm 2 . For LTHFS of 0, 1 and 3, the lower tube remained in the convection regime 
until nucleation of the top tube, as with the smooth tube cases. However for all the 
High Flux tube runs, with the exception of the 0 LTHFS, the lower tube nucleated at 
the same time as the upper tube when both tubes started in the convection regime. This 
appeared to be due to the explosive nature of nucleation with these enhanced surfaces 
which caused vapor to impinge upon the lower tube thereby causing immediate 
nucleation on the lower tube as well; there are probably significant shock waves set up 
within the pool that add to this nucleation of the lower tube as well. Once the upper 
tube nucleated, the three LTHFS curves jumped to a common nucleate boiling curve, 
and demonstrated a significant temperature overshoot, similar to that found by 
Sugiyama [Ref. 1] for a single High Flux Tube. These LTHFS runs followed the same 
nucleate boiling curve for the remainder of the increasing run and for the entire de- 
creasing run; the form of the hysteresis loop agreed very well with Sugiyama’s single tube 
work [Ref. 1]. 

The lower tube nucleated from the beginning of the increasing run for a 
LTHFS of 10 and 25 klV/m Unlike the smooth tubes, the hysteresis loop is not com- 
pletely eliminated for these two LTHFS. These curves demonstrate that the upper tube 
experiences several effects during the increasing run. The segment of the curve between 
500 and 1000 W/m 2 is parallel to the natural convection curves but well to the left (i.e. 
well enhanced), indicating strong convection effects due to the bubble plume from the 
lower tube sliding over the upper tube. This segment of the curve agrees very closely with 
that for the smooth tube in the convection region for the same pitch and the same 
LTHFS. This indicates that the bubbles sweeping over the upper non-nucleating High 
Flux tube have the same enhancement effect as on an upper non-nucleating smooth 
tube, i.e. it makes no difference which tube surface is used because the majority of 'he 
heat transfer is associated with the bubble plume sweeping over the upper tube from the 
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Figure 15. Comparison of Lower Tube Flax Settings for a p/d Ratio of 1.5 for 
Smooth Tube 
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nucleating tube below. This supports Cornwell's bubble sweeping hypothesis at low 
heat fluxes [Ref. 17]. For upper tube heat fluxes above 1000 Wjm 2 , the curves for both 
the 10 and 25 kW/m 2 LTHFS show almost vertical slopes, indicating partial nucleation 
and mixed convection, eventually; at high enough heat flux as they join the same com- 
mon nucleate boiling curve as before with LTHFS of 0, 1 and 3 klV/m 2 All decreasing 
runs for all LTHFS followed the same curve, and agreed closely with both the 0 LTHFS 
and with single tube data (Sugiyama [Ref 1 ]) for decreasing heat flux. There is some 
difference between the curves on the decreasing run at low fluxes; however these differ- 
ences are within the experimental uncertainty which is largest in this region. The above 
suggests that there is no enhancement 

of the upper tube heat transfer coefficient gained by heating the lower tube at a p/d of 

2 . 

e. High Flux Tube p/d Ratio of 1.8 

Figure 17 for a p/d of 1..8 is almost identical to that for a pitch of 2. The 
same effects (as discussed for Figure 16 above) apply equally to Figure 17. Convection 
and nucleate region curves fall on top of each other with the upper tube nucleating at 
almost the same heat fluxes. The first segment between 500 and 1000 W/m 2 for LTHFS 
of 10 and 25 kW/m 2 agrees very closely with the nucleate boiling region for the smooth 
tube for the same pitch. The reasons for this are also discussed above, and demonstrate 
the repeatability of this effect. The excellent agreement between the plots in Figure 17 
and Figure 16 suggests that there is no effect of tube spacing at any heat flux. This will 
become clearer in later figures. 

/. High Flux Tube pld Ratio of 1.5 

Figure 18 shows the plot for a pitch of 1.5 to be almost identical to the 
plots for pitches of 1.8 and 2. The differences in the plot are in the same areas as in the 
previous pitch. The upper tube nucleates at different fluxes than in the previous two 
pitches, again demonstrating the somewhat random nature of ONB. In all other respects 
the plot is identical to the previous two, further reinforcing that there is nothing to be 
gained from an upper enhanced tube by the presence of a nucleating lower enhanced 
tube except if the upper tube is in the natural convection region. 

2. The Effect of Tube Spacing 
a. Lower Tube Unheated 

To show the effect of tube spacing, plots have been made for all three 
pitch- to-diameter ratios at each LTHFS used. For the base of an unheated lower tube, 
the data are shown in Figure 19 for both the smooth and High Flux tubes for all three 
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Figure 16. Comparison of Lower Tube Flux Settings for a p/d Ratio of 2 for High 
Flux Tube 
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Figure 17. Comparison of Lower Tube Flux Settings for a p/d Ratio of 1.8 for High 
Flux Tube 



44 



QDP (W/SQ METER) 

10000 100000 




Figure 18. Comparison of Lower Tube Flux Settings for a p/d Ratio of 1.5 for High 
Flux Tube 
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values of p/d tested. This plot clearly shows the hysteresis loop for the smooth tubes 
(discussed earlier) and the earlier incipient point associated with the High Flux tubes. 
It should be noted that on a log/log plot, the temperature overshoot for the High Flux 
tubes looks much larger than that for the smooth tubes. In fact it is the other way 
around, with a smooth tube temperature overshoot of around 25 K and a High Flux 
value around 15 K. In the natural convection region, all the data fall very closely to- 
gether. Following nucleation of the upper tube, the three curves for each tube jump to 
that tube's respective nucleate boiling curve for the remainder of the run; these curves 
agree well with Sugiyama's [Ref. 1] single tube data for both smooth and High Flux 
tubes. As expected for an unheated tube, there is no effect of tube spacing for High 
Flux tubes for any heat flux. However for the smooth tubes at high heat fluxes, there 
does appear to be a definite effect of tube spacing, with a p/d of 1.8 giving the best heat 
transfer performance. With the bottom tube unheated, it is unclear as to why this 
should be the case (except that its presence obviously effects the circulation patterns in 
the evaporator), but as will be seen, this trend is consistent even when the lower tube is 
heated. 

b. Lower Tube Heat Flux Setting of l k W/m 1 

* 

Results for both tube surfaces at three values of p/d for a 1 kW/m 2 LTHFS 
are shown in Figure 20. In the natural convection region, it appears that both High 
Flux and smooth tubes show a small increasing enhancement with increasing tube 
spacing. This agrees with Sparrow and Niethammers [Ref. 11] and Marsters [Ref. 12] 
work with air. It should be noted that for this LTHFS, both upper and lower tube were 
in the convection region at the start of the increasing run. Figure 20 also shows the 
lower nucleation heat fluxes associated with ONB for High Flux tubes. Once the upper 
tubes nucleate, all curves 'jump' to a nucleate boiling curve. The High Flux tubes follow 
the same nucleate boiling curve regardless of spacing for the remainder of the run. This 
suggests that for a High Flux tube which is nucleating, there is no effect of heating from 
a lower tube with a LTHFS of 1 klV/m 2 or less for any p/d ratio. The smooth tubes on 
the other hand, do show small enhancements throughout the nucleate region for the 
three values of p/d. These are similar to the data shown in Figure 19 when the bottom 
tube was present but not heated, and suggests that these effects are directly related to 
the tube spacing itself and not due to the fact that the lower tube is heated. In 
Figure 20, there does seem to be some effect of p/d at these lower heat fluxes on the 
boiling curve; furthermore, this appears to be opposite to the trend found at high heat 
fluxes, i.e. a p/d of 1.5 gives the best heat transfer. However, this is thought to be ex- 
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perimental scatter as it does not occur on later graphs, whereas the effect at high heat 
fluxes is very repeatable. 

c. Lower Tube Heat Flux Setting of 3 It Wjm 2 

The plot shown in Figure 21 for a LTHFS of 3 kW/m 2 is very similar to that 
for 1 klV/m 2 . The same enhancement trends are seen in the convection region in that 
both tubes show increasing enhancement with increasing pitch, with the exception of the 
smooth tube with a p/d of 1.8. This is believed to be due to scatter because the per- 
formance shown for both tubes for the other LTHFS 

in the natural convection regime tend to follow an increasing enhancement with pitch 
trend. It should be noted for this LTHFS that both lower and upper tube were in the 
convection region at the beginning of the increasing run. Once the upper tube had 
nucleated, then the lower tube also started to nucleate from the 'seeding' process men- 
tioned earlier. Following nucleation, the High Flux runs show the same results as in the 
previous LTHFS. The smooth tube still shows small enhancements at high heat fluxes 
as in the previous LTHFS, with a pitch of 1.8 giving the best enhancement. However 
in the boiling region below lOkW/m 1 , no enhancements are discemable due to p/d (as 
mentioned earlier) and the curves fall directly on top of one another. 

d. Lower Tube Heat Flux Setting of 10 kWjm 2 

It should be noted that for this LTHFS, the lower tube nucleated at the 
beginning of the run and remained that way throughout the run. The plot for a LTHFS 
of 10 kW/m 2 is shown in Figure 22 and clearly shows the elimination of hysteresis for 
the smooth tube. Thus the effect of the bubbles from the lower tube sweeping over the 
upper tube (even when the upper tube is in the convection region) is enough to 'simulate' 
full nucleation from the upper tube, again supporting Cornwell's hypothesis that at low 
heat fluxes, the majority of the heat transfer from upper tubes in a bundle is due to 
convection effects from below. Figure 22 also clearly shows that there is no effect of 
spacing for smooth tubes except at high heat fluxes, as before (above 20 klV/m 2 ). For 
the High Flux tube it is clear that there is no effect of spacing in the decreasing portion 
of the curve as before. It is unclear what influence spacing may have on the increasing 
curve as any effects appear quite random. This suggests that the effect of spacing in this 
region is very small if not negligible. As mentioned earlier, the segment of the increasing 
High Flux curve between 500 and 1000 IVIm 2 is parallel and very close to the same seg- 
ment for the smooth tube curve. Note that for the smooth tube, the curve is the same 
for increasing and decreasing heat flux, suggesting that the top tube is not nucleating 
at these low heat fluxes and that the enhancement seen is due entirely to the sweeping 
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Figure 21. Comparison of Pitches for a Lower Tube Heat Flux Setting of 3 kW/m* 
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bubbles from below. This implies that both tubes experience very similar convective 
effects when influenced by a vigorously boiling tube below. This was discussed in more 
detail earlier. 

e. Lower Tube Heat Flux Setting of 25 kWjm 2 

Figure 23 shows results for a LTHFS of 25 AIF/m 2 and is very similar to 
that for 10 kW/m 2 . It should be noted that the lower tube in each case was nucleating 
from the beginning of the run. The only region where spacing has any effect is in the 
high heat flux boiling region (above 20 kW/m 2 ) for the smooth tube. As in the previous 
case, any effects of tube spacing on the High Flux tube can be considered negligible. 
Again, of particular interest is the segment of the increasing curves between 500 and 
1000 W/m 2 . Here the curves for the smooth and High Flux tube fall on top of one an- 
other, further implyimg that at low heat fluxes (with the lower tube nucleating vigor- 
ously) the convection effects on the two types of surface are the same. Of note is the 
excellent repeatability of the data. 

3. Both Tubes Subject To the Same Heat Flux 

In addition to the data already presented, two additional runs were conducted, 
(one for each type of tube) in which the upper and lower tube were run up and down 
together at the same heat flux setting; these tests were both conducted with a p/d of 2. 
It should be noted that the upper and lower tubes were in the convection region at the 
beginning of the runs. The results of these two tests (for the upper tube only) are shown 
in Figure 24. The High Flux tube data fall on top of previous data, showing no effect 
of the lower tube throughout the entire heat flux range. This was expected due to the 
'repeatability' of the High Flux tube shown in Figure 19 through Figure 23. The 
smooth tube shows a much earlier incipient heat flux for the upper tube. However, there 
is a large step in heat flux at this point due to the upper tube suddenly nucleating at a 
lower than expected heat flux (the cause of this nucleation is unknown and may be just 
a function of the randomness of this phenomenon mentioned earlier). Apart from its. 
early ONB, the smooth tube follows a very similar boiling curve to that shown previ- 
ously for a p/d of 2. At the lowest values of heat flux, the hysteresis 'loop' is still slightly 
open; this is because the lower tube still has a small number of active nucleation sites 
which are slightly enhancing the top tube. At 1 &capflux., the difference between the 
increasing and decreasing values of AT is similar to that found on Figure 20 for a 
LTHFS of 1 klVjm 2 . This also true for the other LTHFS, i.e.the curve on Figure 24 
coincides at one point with the boiling curves in Figure 19 through Figure 23 at the 
respective LTHFS. Of particular interest is that the foregoing results agree with and 
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Figure 22. Comparison of Pitches for a Lower Tube Heat Flux Setting of 10 klV/m* 
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Figure 23. Comparison of Pitches for a Lower Tube Heat Flux Setting of 25 klVfm* 
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complement certain results found from the bundle apparatus in the Heat Transfer labo- 
ratory. Not too many conclusions can be drawn from Figure 24 due to the lack of data 
taken at other spacings. This is one area which should be further investigated in the 
future. 
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VI. CONCLUSIONS 



1. There is no effect of tube spacing on High Flux tube performance at any heat flux. 

2. The effect of tube spacing on smooth tube performance is present (but small) in the 
convection region. At high fluxes (> 20 k Wjm 2 ), there is a systematic dependence on 
tube spacing with a p/d ratio of 1.8 giving the best heat transfer performance. 

3. A nucleating lower High Flux tube significantly enhances an upper High Flux tube in 
the convection region due to bubble sweeping. There is no effect of a heated lower tube 
on a nucleating upper High Flux tube. 

4. A vigorously nucleating lower smooth tube (> 10 kWjm 2 ) eliminates all hysteresis ef- 
fects on an upper tube. A partially nucleating lower smooth tube (<, 3 k Wjm 1 ) signif- 
icantly enhances (by an order of magnitude) an upper smooth tube. 

5. For low upper tube heat flux settings (< 1000 Wjm 2 ) and vigorous nucleation on the 
lower tube, High Flux and smooth tubes exhibited similar performance demonstrating 
that the majority of the heat transfer is due to bubble sweeping from below. 
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VII. RECOMMENDATIONS 



1. Continue present work, with other spacings and 'bundle' configurations. 

2. Adapt the current apparatus to provide a larger evaporator to accommodate more 
than two tubes. This would probably also provide better conditions for videotaping the 
phenomena within the evaporator. 

3. Resume experiments operating both tubes at the same heat fluxes for other spacings. 

4. Pursue modeling of the enhancement effect, taking into account the mixed effects of 
convection (from a lower tube). Studies of a smooth tube placed above an enhanced 
tube might shed light on this. 

5. Measurement of temperature in the plume from the lower tube to study the effect of 
spacing in both the natural convection and boiling region. 

6. Acquisition of a Hewlett-Packard to IBM converter would greatly increase the ability 
to produce results from processed data within a short time frame. 
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APPENDIX A. SAMPLE CALCULATIONS 



Data set DAT0529I52 was chosen for sample calculations. The saturation temper- 
ature was 2.25 °C, the heat flux on the upper tube was 35590 Wjni 1 . 

A. TEST-TUBE DIMENSIONS 

Do = 0.01582 m 
Di = 0.0132 m 
Dl = 0.013 m 
L = 0.190 m 
Lu = 0.076 m 

B. MEASURED PARAMETERS 
V = 147.55 V 

l — 2.44 A 

Tl = 5.88 °C 

n — 4.90 °c 

73 = 5.38 °C 

74 not read, defective 

75 =5.77 °C 

76 = 7.0 6°C 
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77 = 5.66 °C 



7^ = 5.48 °C 
Tsat= 2.25 °C 
Ac = 45 fF/m . AT 

C. OUTER WALL TEMPERATURE OF THE TEST TUBE 
p = n • Do = n • 0.01582 = 0.0497 m 

Ac - n(Do 2 - £>/ 2 )/4 = rc{(0.01582) 2 - (0.0132) 2 }/4 

Ac = 5.97 x 10” 5 m 2 

Qh=VI= 145.55 • 2.44 = 360.02 W 

Tavg — ZTnlm = 5.73 °C 

where m = the number of operational thermocouples used 
Two = Tavg — Qh» { ln(£>o/£>l)/(2 • n • L • kc )} 

Two = 5.73 - 360.02{ ln( ^ )/( 2 • tt • 0.2032 . 45)} 

Two = 4.48 °C 

e= Two- Tsat = 4.48 - 2.25 = 2.23 °C 
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D. PROPERTIES OF R-l 14 AT FILM TEMPERATURE 
77= (Two + Tsat)/2 = (4.48 + 2.25)/2 = 3.36 °C 

H = exp[ -4.4636 + (101 1.47/77)]. 10" 3 

H = exp[ -4.4636 + (1011.47/(3.36 + 273.15))] . 10~ 3 = 3.76 x 10" 4 N • 5/m 2 
7c = Critical Tempera ture(° R) = 753.95 °R 
TJ[°R) = 498.05 °R 

7=1- T/(R)ITc(R) = 1 - 498.05/753.95 = 0.34 
p = 581.77 + 984. 1 5y 1/3 + 263.02; + 279.99y l/2 + 17.94; 2 
p= 1522.7 kg\rn 

v = p/p = 3.7605 x 10” 4 / 1522.72 = 2.47 x 10" 7 m 2 /5 
k = l.lx 10 -2 — {2.61 x 10" 4 . TA°Q) 
k = 7.1 x 10~ 2 - (2.61 x 10 -4 . 3.36) = 7.01 x 10 -2 W/m • K 

Cp = 400 + 1.65 . 77+ (1.51 x 10 -3 • 7/) - (6.68 x 10 -7 . TJ 3 ) 

Cp = 959.86 Jjkg . K 
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a = kip • Cp = 7.0224 x 



1522.72.959.856 



= 4.80 x 10~ 8 m 2 ls 



/? = - ( AplAT)lp = 1.0201 1 x 10 _3 (1 IK) 

Pr = v/a = 5.15 

E. HEAT FLUX CALCULATION 

Using Churchill and Chu s correlation for a cylinder in natural convection, the heat 
transfer coefficient from one non-boiling end can be calculated as: 



m = {(/i • p)l{kc • Ac)}' 12 

and h was computed through an iterative process beginning with an assumed h of 
190 W/m 2 • K. The resulting h and m were: 

A - 110.78 Wjrn.K 
m = 45.28 

<2f= (h'p'kc* Ac)' 12 • 0 . tanh(m . Lu) = 0.27 W 

F. HEAT FLUX THROUGH ACTIVE BOILING SURFACE 

Q = Qh - 2 • Qf= 360.02 - 2(0.27) = 359.47 W 

Ab = n • Do • L = n • 0.01582 • 0.2032 = 1.01 x 10~ 2 m 2 



h = -~ 0.6 + 0.387 




g • B « Do 3 • 0 • tanh(m • Lu) 

v • <x • Lu • m 



1/6 



[1 + (0.559//V) 9 ' 16 ] 8 ' 27 



where 
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q = QlAb = 359.47/1.01 x 10~ 2 = 35595 W/m 2 
h = q'/Q = 35595/2.23 = 15955 Wlm 2 • K 

The following results were produced by the data acquisition and reduction program 
DRP72: 

q = 35590 Wlm 2 

h = 15969 Wlm 2 . K 
0 = 2.23 °C 

The calculations were exactly the same for both tubes. 
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APPENDIX B. UNCERTAINTY ANALYSIS 



For comparison purposes, uncertainty was determined using the same methods 
(those of Kline and McClintock [Ref. 18]) as Sugiyama in single tube work [Ref.l Ap- 
pendix E]. Four data points (two for the smooth tube and two for the High Flux tube) 
were selected for the uncertainty analysis. These points were selected such that uncer- 
tainty could be determined in the high and low heat flux regions for each tube. The 
points selected were data sets 19 and 35 from run DAT0529D52 and data sets 1 and 13 
from run DAT0425D42. The following is a sample calculation of uncertainty for data 
set 19 of run DAT0529D52 (i.e. low heat flux setting on the High Flux tube). Measured 
and calculated parameters were obtained as in Appendix A, sample calculations. All 
uncertainties are expressed as a percentage of the calculated parameter. Results are shown 
in Table 5 

A. UNCERTAINTY IN INPUT POWER. 

Qh = VI 



Is = 0.37 V SI =±0.025/1 

where <5 = uncertainty in measurement and calculation 
Vs — 1.75 V SV = ± 0.05 V 

I = 1.92 . Is = 0.7U 
V =25 •Vs = 43.78 V 

SQhlQh = ((5 VI Vs) 2 + (Sills) 2 ) 1 ' 2 
SQhlQh = ((0.05/ 1.75) 2 + (0.025/0. 37) 2 ) 1/2 



SQhlQh = 7.33 
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B. UNCERTAINTY IN SURFACE AREA 



Ab = 7i • Do • L 

Do = 15.82 mm SDo = 0.1mm 

L = 203.20mm 5L = 0.1mm 
SAblAb = (( 5 Do! Do) 2 + (<5L/L) 2 ) 1/2 

5 Ab/Ab = ((0.1/15.82) 2 + (0.1/203.2) 2 ) t/2 
SAb/Ab = 0.63 

C. UNCERTAINTY IN WALL SUPERHEAT 

AT = Two — Tsat 

Tsat = 2.27° C 5 Tsar = 0.01 °C 

Two = 7avg - QA[( ]n(DojDl))l(2 • *.£,. kc )] 

Tn = thermocouple readings 



T1 = 2.76°C 


72 = 2.61°C 


73 = 2.67° C 


74(defective thermocouple) 


75 = 2. 74° C 


76 = 2.81°C 


77 = 2.76°C 
7avg = (L7«/7) 


78 = 2.74°C 
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where n = 1 to 3, 5 to 8 



Tavg- 2.73° C 

S.D. = ((1(7* - Tavg) 2 )ln)' 12 = 0.06l°C 
where S.D. = standard deviation 

Compared to Tavg, the logarithmic term in the equation for Two is small and is 
neglected for the uncertainty analysis. 

Two = Tavg = 2.71° C S Two = S.D. = 0.061 °C 
AT = 0.17° C 

6AT/AT- (( STwo/AT) 2 + (<5 Tsat / AT) 2 ) 1 ' 2 
SAT/ AT = ((0.061/0.368) 2 + (0.01/0.368) 2 ) 1 ' 2 
SAT/AT = 16.8 

D. UNCERTAINTY IN HEAT FLUX 

q = {Qh - 2 • Qf)lAb 

Qh = 11 .08 W SQh = 2.28 W 

Assuming the same proportion in the uncertainty for Qf (losses from the unheated ends): 
Qf = 1.61V SQf = 0.1 2 IF 

Qh-2- Qf = 27.9 W 

5q"lq" = {(<5 QhliQh - 2 . Qf)) 2 + (2 • QJ]{Qh -2.QJ)) 2 + (dAb/Ab) 2 }' 12 
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bqlq = {(2.28/27. 88) 2 + (2 . 0.12/27.88) 2 + (0.0063) 2 } 1/2 
bqfq = %.2\ 

E. UNCERTAINTY IN BOILING HEAT TRANSFER COEFFICIENT 
h = q'l&T 

bh/h = {( bq'/q ) 2 + (MT/AT) 2 } 1 ' 2 
= {(0.0821) 2 + (0.168) 2 } i/2 
18.7 



Table 5. UNCERTAINTY ANALYSIS FOR FOUR 1 


POINTS 


Parameters 


Smooth tube 
2685 Wlm 1 


Smooth tube 
38440 Wlm 1 


High Flux tube 
3070 W/m> 


High Flux tube 
35590 W/m 2 


SQhlQh (%) 


7.01 


1.98 


7.33 


2.14 


SAb/Ab (%) 


0.63 


0.63 


0.63 


0.63 


(5AT/AT (%) 


1.91 


2.52 


16.8 


8.08 


W q (%) 


7.79 


2.18 


8.21 


2.23 


Sh/h (%) 


8.02 


3.33 


18.7 


9.40 



Uncertainties listed in Table 5 show that the main source of uncertainty is in the 
measurement and calculation of AT. The uncertainty in AT is probably due to inaccu- 
racies introduced during the manufacturing process. For instance, the solder pools at- 
taching thermocouple wires to the interior sleeve may have impurities (e.g. small air 
gaps) which in close proximity to the thermocouple junction could affect the local heat 
resistance and thereby the temperature. Eight thermocouples were used to obtain the 
wall average and reduce this uncertainty. This also points out that further improvements 
in accuracy must be made in this area first. Graphical displays of uncertainty are shown 
in Figure 10 and Figure 11. 
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APPENDIX C. 



REPRESENTATIVE DATA SET 
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Oate : 



6 Apr 1992 



NOTE: Program name : QRP72 
Oisk number * 00 
New file name: OAT0406I42 
TC is defective at location S 
No defective AUX TCs exist 
Tube Number: 4 



Oata Set 
TIME: 

TC No: 
Temp : 

TC No: 
Temp : 
Twa 
2.62 
Thetab 
.356 



Number * 
16:09:03 



I Bulk Oil 7* 



0.0 



1 

2.61 

9 

16.87 
ATwa 
21 .59 
Htube 



2 

2.62 ; 
10 

27.84 z: 
T liqd 
2.28 

Qdp 



3 

.64 
I I 

:. 1 1 
Tliqd2 
2.21 



4 

2.62 

12 

26.61 



6 

2.61 
I 4 

2 7.90 



7 

2.62 

IS 

18.51 



1 . 71 6E+03 6.1 12E+02 



5 

-99.99 
13 

10.17 

Tvapr Psat Tsump 
2.28 -1.70 -16.5 

Athetab Ahtube AuQdp 
19.325 5. 093E+02 9.842E+03 



8 

2.66 

16 

23.24 



Oata Set 


Number 


a 2 


Bulk Oil 7. - 0.0 






TIME: 


16:10 


: 5 1 








TC No.- 


1 


2 


3 4 5 


6 7 


8 


Temp : 


2.71 


2.72 


2.73 2.70 -99.99 


2.63 2.68 


2.73 


TC No: 


9 


10 


11 12 13 


14 15 


16 


Temp : 


16.87 


27.86 


22.10 26.60 10.15 


27.89 18.46 


23.21 


Twa 


ATwa 


T1 iqd 


Tliqd2 Tvapr Psat 


T sump 




2.70 


21 .58 


2.31 


2.24 2.26 -1.69 


-16.5 




Thetab 


Htube 


Qdp Athetab Ahtube AuQdp 




.437 


1 . 384E+03 6, 


.048E+03 19.309 S.O74E+02 9.798E+03 



Oata Set 


Number 


* 3 


8u Ik 


Oil 7. - 0.0 






TIME: 


16:15 


: 44 










TC No: 


1 


-i 


3 


4 5 


6 7 


8 


Temp : 


3.42 


3.50 


3.49 


3.47 -99.99 


3.40 3.40 


3.S2 


TC No: 


9 


10 


1 1 


12 13 


14 15 


16 


Temp : 


16.86 


27.86 


22.07 


26.63 10.08 


27.94 18.50 


23.25 


Twa 


ATwa 


T liqd 


Tliqd2 Tvapr Psat 


T sumo 




3.45 


21.58 


2.28 


2.21 


2.24 -1.72 


-16.4 




Thetab 


Htube 


Qdp 


Athetab Ahtube AuQdp 





1.206 1 . 307E+03 1 . 577E+03 19.342 S.I00E+02 9.864E+03 



Oata Set 


Number 


* 4 


8ulk Oil 7. » 0.0 








TIME: 


16:16 


: 2 1 










TC No: 


1 


n 


3 4 5 


6 


7 


8 


Temp : 


3.43 


3.52 


3.51 3.48 -99.99 


3.40 


3.40 


3.54 


TC No: 


9 


10 


11 12 13 


14 


15 


16 


Temp : 


16.87 


27.85 


22.07 26.62 10.12 


27.93 


18.51 


23.24 


Twa 


ATwa 


T liqd 


Tliqd2 Tvapr Psat 


T sump 




3.46 


21 .59 


2.31 


2.22 2.27 -1.69 


-16. 


, 4 




Thet ab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




1.189 


1 . 324E+03 l. 


•S74E+03 19.314 5 . 1 0GE+02 


9 . 862E+03 


Data Set 


Number 


- 5 


Bulk Oil 7. =« 0.0 








TIME: 


16:22 


: 00 










TC No: 


1 


2 


3 4 5 


6 


7 


9 


T emp : 


4.47 


4.66 


4.64 4.60 -99.99 


4 , 45 


4 . 44 


4.65 


TC No: 


9 


10 


11 12 13 


1 4 


15 


16 
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T emp : 


16.92 


27.79 


21.99 26.55 10.08 


27.86 


18.44 


23.18 


Tua 


ATua 


T I lqd 


TliqdZ Tvapr Paat Taump 




4.54 


21 .52 


2.25 


2.18 2.23 -1. 


74 -16, 


,3 




Thet ab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




2.321 


1.270E+03 2. 


■948E+03 19.301 5. 


106E+0Z 


9 . 855E+03 


Oata Set 


Number 


- 6 


8ulk Oil S • 0.0 








TIME : 


16:22 


: 4 1 










TC No: 


1 


2 


3 4 5 


6 


7 


8 


Temp : 


4.51 


4.65 


4.62 4.57 -99.99 


4.45 


4.44 


4.64 


TC No: 


9 


10 


11 12 13 


14 


15 


16 


T emp : 


16.81 


27.79 


21.98 26.55 10.06 


27.86 


18.45 


23.17 


T wa 


ATua 


T I iqd 


TIiqd2 Tvapr Paat Taump 




4.54 


21.52 


2.29 


2.19 2.23 -1. 


73 -16. 


.2 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




2.305 


1 . 2796+03 2. 


, 948E+03 19.203 5. 


1 1 6E+02 


9 . 8656+03 



Oata Set 


Number 


- 7 


3uik Oil X - 0.0 






TIME: 


16:27 


: 34 








TC No: 


1 


2 


3 4 5 6 


7 


8 


Temp : 


6.63 


6.95 


6.89 6.82 -99.99 6.48 


6.48 


6.90 


TC No: 


9 


10 


11 12 13 14 


15 


16 


Temp : 


16.79 


27.80 


21.97 26.57 10.00 27.88 


18.43 


23.19 


T ua 


ATua 


T 1 iqd 


TliqdZ Tvapr Paat Taump 




6.70 


21 .51 


2.27 


2.19 2.20 -1.75 -16. 


,2 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




4.484 


1. 2696+03 5. 


.6896+03 19.296 5.1136+02 


9 . 367E+03 


Oata Set 


Number 


- 8 


8uik Oil X - 0.0 






TIME: 


16:28 


: 36 








TC No: 


1 


2 


3 4 5 6 


7 


8 


T emp : 


6.58 


6.95 


6.89 6.81 -99.99 6.55 


6.55 


6.95 


TC No: 


9 


10 


11 12 13 14 


15 


16 


T emp : 


16.89 


27.88 


22.06 26.65 10.20 27.96 


18.54 


23.27 


Tua 


AT ua 


T I iqd 


TIiqd2 Tvapr Paat Taump 




6.72 


21 .61 


2.30 


2.26 2.32 -1.66 -16. 


, 1 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




4.419 


1 . 292E+03 5. 


.7086+03 19.313 S.13IE+02 


9.909E+03 



Oata Set 


Number 


- 9 


8ulk Oil X - 0.0 






TIME: 


16:30 


: 1 9 








TC No: 


1 


2 


3 4 5 6 


7 


8 


T emp : 


6.64 


6.99 


6.94 6.86 -99.99 6.62 


6.60 


6.97 


TC No: 


9 


10 


11 12 13 14 


15 


16 


Temp : 


16.39 


27.90' 


22.08 26.67 10.19 27.99 


18.55 


23.30 


T ua 


ATwa 


T I iqd 


TliqdZ Tvapr Paat Taump 




6.77 


21.63 


2.36 


2.28 2.33 -1.63 -16. 


1 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




4.443 


1.285E+03 5.71 0E+03 19.304 5.1326+02 


9 . 906E+03 


Oata Set 


Number 


- 10 


Bulk Oil X - 0.0 






TIME: 


16:36 


: 47 








TC No: 


1 


n 


3 4 5 6 


7 


8 


T emp : 


8.57 


in 

cn 


9.03 8.95 -99.99 3.48 


8.51 


9.06 


TC No: 


9 


10 


11 12 13 14 


15 


16 


T emp : 


16.90 


27.92 


21.99 26.. 59 10.03 27.90 


13.44 


23.21 


Tua 


ATua 


T 1 iqd 


Tliqd2 Tvapr Paat Taump 




8.77 


21 .53 


2.27 


2.20 2.25 -1.72 -15. 


9 




Thetab 


Htube 


Qdo Athetab Ahtube 


AuQdp 




6.528 


1 . 292E+03 9, 


.43IE+05 19. 399 5.I3SE.0: 


9. 904E+03 
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Oata Set 


Number 


■ l 1 


8ulk Oil X * 0.0 








TIME: 


16:37 


: 3 1 










TC No: 


1 




5 4 5 


6 


7 


0 


Temp : 


8.67 


9.15 


9.07 8.96 -99.99 


8.57 


8.54 


9.04 


TC No: 


9 


10 


11 12 13 


1 4 


15 


16 


T emo • 


16.03 


27.05 


22.01 25.61 10.04 


27.93 


10.48 


23.24 


Tua 


ATua 


T1 iqd 


Tliqd2 Tvapr Paat Taump 




■8.01 


21 .56 


2.28 


T 1 *1*7 _ J 


71 -15. 


9 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




6.549 


1.298E+03 8. 


433E+03 19.302 5. 


141E+02 


9.924E+03 


Data Set 


Number 


- 12 


8ulk Oil 7. - 0.0 








TIME: 


16:46 


:52 










TC No: 


1 


2 


3 4 5 


6 


7 


8 


T emp : 


10.81 


11.49 


11.44 M .32 -99.99 


10.67 


10.70 


11.45 


TC No: 


9 


10 


11 12 13 


1 4 


15 


16 


Temp : 


16.85 


27.92 


22.33 26.67 10.02 


28.00 


18.49 


23.28 


Tua 


ATua 


Ti iqd 


Tliqd2 Tvapr Paat Taump 




t 1 .06 


21 .59 


2.32 


2.26 2.28 -1. 


67 -15. 


5 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




8.771 


1.3UE+03 1. 


152E+04 19.305 5. 


IS2E+02 


9 . 946E+03 



Oata Set 


Number 


- 13 


8uik Oil 7. - 0.0 






TIME: 


16:47 


: 44 








TC No.- 


1 


2 


3 4 S 6 


7 


8 


Temp : 


10.76 


11.52 


11.44 It .21 -99.99 10.73 


10.71 


1 1 .36 


TC No: 


9 


10 


11 • 12 13 14 


15 


16 


Temp : 


16.85 


27.89 


22.01 26.65 10.01 27.97 


18.49 


23.25 


Tya 


ATua 


Tl iqd 


Tiiqd2 Tvapr Peat Taump 




1 1.03 


21 .57 


2.32 


2.24 2.27 -1.68 -15. 


,4 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




8.759 


1.318E+03 1. 


.IS4E+04 19.296 5 . 1 66E+02 


9. 969E+03 


Oata Set 


Number 


- 14 


0uik Oil 7. - 0.0 






TIME: 


16:57 


:5S 








TC No: 


1 


2 


3 4 5 6 


7 


3 


T emp : 


12.48 


13.36 


13.26 13.10 -99.99 12.29 


12.30 


13. 1 9 


TC No: 


9 


10 


11 12 13 14 


15 


16 


T emp : 


16.00 


27.91 


21.99 26.66 9.97 27.99 


18.45 


23.25 


Tya 


ATua 


Tl iqd 


TUqd2 Tvapr Paat Taump 




12.77 


21 .56 


2.23 


2.18 2.18 -1.77 -14. 


,9 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




10.576 


1 . 332E+03 1, 


. 403E+04 19.366 5 . 1 58E+02 


9. 989E+03 


Oata Set 


Number 


- 15 


8ulk Oil 7. - 0.0 






TIME: 


16:58 


:31 








TC No: 


1 




3 4 5 6 


7 


8 


T emp : 


12.40 


13.36 


13.29 13.14 -99.99 12.46 


12.46 


1 3“. 36 


TC No : 


9 


10 


11 12 13 14 


15 


16 


T emp : 


16.94 


27.32 


22.02 26.68 10.00 28.02 


18.49 


23.29 


Tya 


ATua 


T l iqd 


Tliqd2 Tvapr Paat Taump 




12.84 


21 .59 


2.24 


2.19 2.22 -1.74 -14.8 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




10.619 


1.330E+03 1.412E+04 19.369 5.171E+02 


l . 002E+04 


Oata Set 


Number 


- 16 


8ulk Gil 7. * 0.0 






TIME: 


17:05 


:00 








TC No: 


1 




3 4 5 6 


7 


9 


T emp : 


14.01 


15.25 


15. U 14.72 -99.99 14.13 


14.08 


15.06 
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TC No: 


9 10 


1 1 


12 13 14 


15 


16 


T enp : 


16.84 27.54 


22.04 26 


.69 10.02 29.02 


18.49 


23.26 


Tua 


ATua Tliqd 


T I iqd2 


Tvapr Paat Taunp 




14.52 


21.59 2.23 


2.23 


2.19 -1.74 -14.5 




Thetab 


Htube Odp 


Athetab Ahtube 


AuOdp 




12.500 


1.3G5E+03 1. 


, 679E+04 


19.370 5 . 1 70E+02 


1 . 00 1 E +04 


Data Set 


Nunber * 17 


8ulk Oil 


7. - 0.0 






TIME: 


17:05:42 










TC No: 


1 2 


3 


4 5 6 


7 


8 


Tenp : 


14.12 15.39 


15.28 14 


.77 -99.99 14. 15 


14.16 


15.10 


TC No: 


9 10 


1 1 


12 13 14 


15 


16 


Tenp : 


16.83 27.93 


22.02 26 


.69 9.98 28.02 


18.49 


23.28 


Tua 


ATua Tliqd 


Tl iqd2 


Tvapr Paat Taunp 




14. G1 


21.59 2.27 


2. 18 


2.20 -1.75 -14. 


4 




. Thetab 


Htube Odp 


Athetab 'Ahtube 


AuQdp 




12.398 


1. 3546+03 1 .5796+04 


19.376 5 . 1 54E+02 


9. 986E+03 


Oata Set 


Nunber ■ 18 


Bulk Oil 


7. - 0.0 






TIME: 


17:10:45 










TC No: 


1 2 


3 


4 5 6 


7 


8 


T enp : 


15.64 16.91 


cn 

CO 

<s> 

cn 


.36 -99.99 16.07 


16.01 


16.72 


TC No: 


9 10 


1 1 


12 13 14 


15 


16 


Tenp : 


16.68 27.63 


21.81 26 


.40 9.92 27.71 


18.31 


23.03 


Tua 


ATua Tliqd 


Tl iqd2 


Tvapr Paat Taunp 




16.24 


21.37 2.25 


2.16 


2.14 -1.80 -14. 


2 • 




Thetab 


Htube Qdp 


Athetab Ahtube 


AuQdp 




14.067 


1 . 406E+03 1. 


. 978E+04 


19.198 5. 1 47E+02 


9.881 E+03 


Oata Set 


Nunber • 19 


3ulk Oil 


% - 0.0 






TIME: 


17:11 :22 










TC No: 


1 2 


3 


4 5 6 


7 


8 


Tenp : 


15.54 15.84 


16.73 16 


.36 -99.99 16.06 


15.96 


16.80 


TC No: 


9 10 


1 1 


12 13 14 


15 


16 


Tenp : 


16.58 27.54 


21.80 26 


.41 9.94 27.72 


18.32 


23.04 


Tua 


ATua Tliqd 


Tltqd2 


Tvapr Paat Taunp 




16.21 


21.38 2.25 


2.15 


2.15 -1.80 -14. 


1 




Thetab 


Htube Odp 


Athetab Ahtube 


AuQdp 




14.036 


I.407E+03 1, 


. 975E+04 


19.207 5. 1 39E+02 


9. 870E+03 



Oata Set 


Nunber - 20 


Bulk Oil 7 . - 0.0 








TIME: 


17:16:58 










TC No: 


1 2 


3 4 5 


6 


7 


8 


Tenp : 


17.40 18.39 


18.23 17.64 -99.99 


17.89 


17.62 


18.27 


TC No: 


9 10 


11 12 13 


1 4 


15 


16 


Tenp : 


16.77 27.72 


21.88 26.50 9.98 


27.79 


18.39 


23.08 


Tua 


ATua Tliqd 


Tiiqd2 Tvapr Paat Taunp 




17.78 


21.45 2.32 


2.25 2.20 -1. 


72 -13. 


,9 




Thetab 


Htube Qdp Athetab Ahtube 


AuOdp . 




15.539 


1.51 6E+03 2, 


.35SE+04 19.206 5. 


1 48E+02 


9. 986E+03 


Data Set 


Nunber » 21 


8ulk OH Z - 0.0 








TIME: 


17:17:40 










TC No: 


1 2 


3 4 5 


6 


7 


8 


Tenp : 


17.39 18.33 


18.18 17.62 -99.99 


17.82 


17.64 


18.24 


TC No: 


9 10 


11 12 13 


1 4 


15 


16 


Tenp : 


16.77 27.74 


21.91 26.52 10.02 


27.82 


18.42 


23. 1 1 


Tua 


A Tua Tliqd 


TliqdZ Tvapr Paat Taunp 




17.75 


21.47 2.32 


2.25 2.19 -1. 


73 -13. 


. 9 




The tab 


Htube Odp Athetab Ahtube 


AuOdp 
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IS. 511 1.SZ2E+03 2.360E*04 19.234 5.1596*02 



9.9Z3E+03 



Data Set 
TIME : 

TC No: 
Tenp : 

TC No: 
Tenp : 
Tua 
18.65 
Thetab 
16.455 

Data Set 
TIME: 

TC No: 
Tenp : 

TC No: 
Tenp' : 
Tua 
18.63 
Thetab 
16.442 

Data Set 
TIME: 

TC No: 
Tenp : 

TC No: 
Tenp : 
Tua 

19.23 
Thetab 
17.029 

Data Set 
TIME: 

TC No: 

T enp : 

TC No: 
Tenp : 
Tua 

19.24 
Thetab 
17.010 

Oata Set 
TIME: 

TC No: 
Tenp : 

TC No: 

T enp : 
Tua 
19.84 
The tab 
17.705 

Data Set 
TIME: 

TC No: 



Nunber ■ 22 

17:23:15 
1 • 2 


Bulk Oil X - 0.0 

3 4 5 


6 7 


8 


18.56 


18.95 


19.94 18.14 -99.99 


19.56 18.75 


18.83 


9 


10 


11 12 13 


14 15 


16 


16.77 


27.73 


21.90 26.50 10.07 


27.81 18.40 


23. 1 1 


ATua 


T 1 lqd 


TliqdZ Tvapr Psat 


T sunp 




21 .47 


2. 28 


2.25 2.12 -1.77 


-13.5 





Htube QdQ Athetab Ahtube AuQdp 

I.725E+03 2 . 838E+04 19.274 S.138E+02 9.903E+03 



Nunber 


- 23 


8u lk 


Oil 7. - 0.0 






17:23 


: 57 










1 




3 


4 5 6 


7 


8 


18.50 


18.93 


18.38 


18.11 -99.99 19.63 


18.72 


18.81 


9 


10 


1 1 


12 13 14 


15 


16 


16.77 


27.74 


21.90 


26.52 10.05 27.83 


18.42 


23.12 


AT ua 


T 1 iqd 


TliqdZ Tvapr Psat Tsunp 




21 .48 


2.28 


2.22 


2. 12 -1.78 -13. 


,5 




Htube 


Qdp 


Athetab Ahtube 


AuQdp 




1.729E+03 2. 


. 843E+04 19.291 5. 1 26E+02 


9. 889E+03 


Nunber 


- 24 


Bulk 


Oil 7. - 0.0 




• 


17:29 


: 40 










1 


2 


3 


4 5 6 


7 


8 


19.09 


19.38 


19.33 


18.42 -99.99 21.01 


19.60 


19.18 


9 


10 


1 1 


12 13 14 


15 


16 


16.83 


27.80 


21.96 


26.58 10.09 27.87 


18.47 


23.14 


ATua 


T 1 iqd 


TliqdZ Tvapr Psat Tsunp 




21 .53 


2.29 


2.18 


2.17 -1.76 -13. 


, 1 




Htube 


Qdp 


Athetab Ahtube 


AuQdp 





1 . 969E+03 3 . 353E+04 19.320 S.1Z2E+02 9.896E+03 

Nunber - 25 8ulk Oil X 9 0.0 

17:30:27 



1 


n 


3 


4 5 


6 


7 


8 


19.10 


19.42 


19.32 


18.42 -99.99 


21 .02 


19.60 


19.16 


9 


10 


1 1 


12 13 


14 


15 


16 


16.84 


27.81 


Z1 .97 


26.60 10.11 


27.89 


18.48 


23.18 


ATua 


T 1 iqd 


Tliqd2 Tvapr Psat Tsunp 




21 .54 


2.32 


2.21 


2.19 -1. 


74 -13. 


1 




Htube 


Qdp 


Athetab Ahtube 


AuQdp 




1 . 968E+03 3. 


. 348E+04 19.317 5. 


1 07E+02 


9 . 865E+03 


Nunber 


* 26 


8ulk 


Oil 7, - 0.0 








17:35 


: 53 












I 


n 


3 


4 5 


6 


7 


8 


19.60 


19.79 


17.81 


18.74 -99.99 


22.75 


20.47 


19.48 


9 


10 


1 1 


12 13 


14 


15 


16 


16.81 


27.75 


21 .93 


26.54 10.06 


27.83 


CD 


23. 13 



ATua Tliqd TliqdZ Tvapr Psat Tsunp 

21 .49 2.21 Z. 10 2.11 -1 .84 -12.7 

Htube Qdp Athetab Ahtube AuQdp 

2 . 4495+03 4 . J35E+04 19.361 5.096E*02 9.866E+0I} 

Nunber » 27 8u Ik Oil X * 0.0 

17:36:50 

1 2 3 4 5 6 7 8 
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Tenp : 19.62 19.91 19.82 18.76 -99.99 22.76 20.54 19.55 

TC No: 9 10 11 12 15 14 15 16 

Tenp : 16.37 27.80 21.98 26.58 10.18 27.88 18.51 25.18 

Twa ATwa Tliqd Tliqd2 Tvapr Psat Tsunp 

19.37 21.55 2.55 2.19 2.20 -1.75 -12.6 

Thetab Htube Odp Athetab Ahtube AuQdp 

17.659 2 . 456E+05 4.332E+04 19.526 5.M0E+02 9.875E+03 



Data Set Nunber » 

TIME: 17:40:06 

TC No: 1 

Tenp : 

TC No : 9 

Tenp : 16.89 
Tuia ATwa 
20.50 21.58 

Thetab Htube Qdp 

18.277 3.312E+05 6.054E+04 



7 8 

20.00 

15 16 

18.55 23.22 

Tvapr Psat Tsunp 
2.14 -1,75 -12.3 

Athetab Ahtube AuOdp 
19.359 5 . 083E+02 9.839E+03 



28 8ulk Oil % - 0.0 

3 4 5 6 

20.25 20.32 20.40 19.23-99.99 24.53 21.18 

10 11 12 13 14 

27.83 21.99 26.60 10.17 27.91 

Tliqd Tiiqd2 
2.32 2.28 



Oata Set Nunber » 29 Bulk Oil X » 0.0 

TIME: 17:40:48 

TC No: 1 2 3 4 5 6 7 8 

Tenp : 20.23 20.30 20.38 19.22-99.99 24.51 21.17 20.00 

TC No : 9 10 11 12 13 14 15 16 

Tenp : 16.85. 27.74 21.93 26.51 10.13 27.82 18.47 23.15 

Twa ATwa Tliqd Tliqd2 Tvapr Peat Taunp 

20.48 21.51 2.33 2.25 2.16 -1.74 -12.2 

Thetab Htube Odp Athetab Ahtube AuOdp 

18.254 3.306E+03 6.034E+04 19.283 5.096E+02 9.826E+03 



Oata Set 


Nunber * 30 


Bulk Oil 


X - 0.0 








TIME: 


17:41 : 22 • 












TC No: 


1 2 


3 


4 5 6 




7 


8 


Tenp : 


20.23 20.31 


20.38 19.22 -99.99 24.51 


21 


.17 


20.00 


TC No: 


9 10 


1 1 


12 13 14 




15 


16 


Tenp : 


16.86 27.77 


21.95 26 


.55 10.14 27.85 


18 


t.S0 


23.16 


Twa 


ATwa Tliqd 


Tliqd2 


Tvapr Psat Tsunp 






20.48 


21.53 2.30 


2.27 


2.13 -1.76 -12. 


.2 






Thetab 


Htube Qdp 


Athetab Ahtube 


AuOdp 




18.275 


3.310E+03 6, 


+ 

UJ 

cn 

s> 


19.323 5. 1 00E+02 


9. 


8545+03 



NOTE: 29 data rune were stored in file DAT0406I42 



Date : 6 Apr 1992 



NOTE: Progran nane : QRP72 
Disk nunber * 00 
New file nane: OAT0406D42 
TC is defective at location 5 
No defective AUX TCs e*ist 
Tube Nunber: 4 



Oata Set 
TIME: 

TC No: 


Nunber » 1 

17:48:16 
1 2 


8ulk Oil x - 0.0 

3 4 5 


6 7 8 


Tenp : 


19.64 


19.35 


19.85 18.76 -99.99 


22.85 20.54 19.54 


TC Mo: 


9 


10 


11 12 13 


14 15 16 


T enp : 


16.80 


27.70 


21.39 26.48 10.07 


27.79 18.44 23.12 


Twa 


ATwa 


T 1 iad 


Tliqd2 Tvapr Psat 


T sunp 


1 9.89 


21.47 


2.31 


2.23 2.M -1.78 


-1 1 .7 
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Thetab Htube Qdp Athetab Ahtube AuQdp 

17.700 2.48IE+03 4 , 392E+04 19.279 5.110E+02 9.850E+03 



Oata Set 


Number 


- 2 


8ulk Oil X - 0.0 






TIME: 


17:43 


: 07 








TC No: 


1 


2 


3 4 5 6 


7 


8 


Temp : 


19.64 


19.38 


19.88 18.78 -99.99 22.85 


20. S2 


19.56 


TC No: 


9 


10 


II 12 13 14 


15 


16 


Temp : 


16.79 


27.66 


21.87 26.45 10.11 27.75 


18.42 


23.07 


Twa 


ATwa 


T1 tqd 


Tliqd2 Tvapr Paat Taump 




13.80 


21.45 


2.32 


2.28 2.14 -1.75 -11, 


,7 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




I7.S85 


2.476E+03 4.378E+04 19.231 S.II2E+02 


9 . 830E+03 


Oata Set 


Number 


- 3 


8ulk Oil 7. - 0.0 






TIME: 


17:54:28 








TC No: 


1 


2 


3 4 5 6 


7 


8 


Temp : 


18.55 


19.04 


18.96 18.18 -99.99 19.76 


18.77 


18.87 


TC No: 


9 * 


10 


11 12 13 14 


15 


16 


Temp : 


16.70 


27.60 


21.79 26.38 9.99 27.68 


18.33 


23.09 


Twa 


ATwa 


T1 iqd 


Tiiqd2 Tvapr Paat Taump 




18.70 


21 .38 


2.23 


2.19 2.09 -1.82 -11. 


.5 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




16.552 


1.7SSE+03 2, 


. 905E+04 19.226 S.12IE+02 


9. 845E+03 


Oata Set 


Number 


* 4 


Bulk Oil X - 0.0 






TIME: 


17:55 


: 10 








TC No: 


1 


2 


3 4 5 6 


7 


8 


Temp : 


19.61 


19.05 


18.96 18.17 -99.99 19.76 


18.85 


18.89 


TC No: 


9 


10 


11 12 13 14 


15 


16 


Temp : 


16.71 


27.63 


21.80 26.41 9.99 27.71 


18.35 


23.12 


Twa 


A Twa 


T l iqd 


Tliqd2 Tvapr Paat Taump 




18.73 


21 .40 


2.26 


2.20 2.12 -1.79 -11. 


.5 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




IS.S5I 


1 . 756E+03 2, 


. 906E+04 19.222 5.122E+02 


9.345E+03 


Oata Set 


Number 


- 5 


8ulk Oil X - 0.0 






TIME: 


18:01 


: 43 








TC No: 


1 


2 


3 4 5 6 


7 


8 


Temp : 


17.47 


18.37 


18.30 17.72 -99.99 17.94 


17.61 


18.35 


TC No: 


9 


10 


11 12 13 14 


• 15 


16 


Temp : 


16.70 


27.60 


21.78 26.39 9.96 27.67 


18.32 


23.07 


Twa 


ATwa 


T l iqd 


Tliqd2 Tvapr Paat Taump 




17.82 


21 .37 


2.23 


2.18 2.09 -1.83 -11. 


.4 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




IS. 678 


1 . 555E+03 2, 


. 439E+04 19.225 S.I40E+02 


9. 882E+03 


Oata Set 


Number 


6 


Sulk Oil X - 0.0 






TIME: 


18:02 


: 25 








TC No: 


1 


2 


3 4 5 6 


7 


8 


T emp : 


17.49 


18.30 


18.25 17.75 -99.99 17.96 


17.62 


18.23 


TC No: 


9 


10 


II 12 13 14 


IS 


16 


Temp : 


16.75 


27.71 


21.86 26.50 9.99 27.79 


18.39 


23.18 


Twa 


ATwa 


Ti iqd 


Tliqd2 Tvapr Paat Taump 




17.80 


21 .45 


2.25 


2.22 2.10 -1.80 -11. 


. 4 




Thetab 


Htube 


Qdp Athetab Ahtube 


AuQdp 




15.627 


1.56ZE+03 2, 


, 4 42 E a 04 19.295 5 . 1 53E+0Z 


9 . 938E+03 



Oata Set Number ’ 7 Bulk OtI % * 0.0 

TIME: 18:07:14 
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TC No: 


1 


-> 


3 4 5 


6 7 


8 


Temp : 


15.29 


16.40 


16.31 15.94 -99.99 


15.69 15.53 


16.35 


TC No: 


9 


10 


11 12 13 


14 15 


16 


T emp : 


15.74 


27.72 


21.84 26.50 9.89 


27.90 13.38 


23.19 


T ua 


ATua 


T l iqd 


Tliqd2 Tvapr Psat 


T sump 




15.92 


21.44 


U . 4 - 


2.13 2.14 -1 .31 


-11.4 





Thetab Htube Qdp Athetab Ahtube AuQdp 

1 3 . 656 1.402E+03 1.91SE+04 19.381 S.1S3E+02 9.937E+03 

Data Set Number - 8 8ulk Oil X * 0.0 



TIME: 


18:07:50 










TC No: 


1 2 


3 4 5 


6 


7 


8 


T emp : 


15.33 16.51 


16.43 15.99 -99.99 


15.59 


15.55 


16.40 


TC No: 


9 10 


11 12 13 


1 4 


15 


16 


Temp : 


16.75 27.73 


21.37 26.50 9.94 


27.31 


18.39 


23.20 


T ua 


ATua Tliqd 


Tliqd2 Tvapr Psat 


T sump 




15.86 


21.46 2.23 


2.16 2.14 -1.30 


-1 1 . 


,4 




Thetab 


Htube Qdp Athetab Ahtube 


AuQdp 




13.690 


1.397E+03 1, 


•9I2E+01 19.290 S . 1 39E+02 


9. 91 2E+03 


Oata Set 


Number - 9 


8ulk Oil X - 0.0 








TIME: 


18:12:30 










TC No: 


1 2 


3 4 5 


6 


7 


8 


Temp : 


12.40 13.42 


15.33 13.03 -99.99 


12.45 


12.33 


13.24 


TC No: 


9 10 


11 12 13 


14 


IS 


16 


Temp : 


17.23 28.71 


22.91 27.44 10.37 


28.81 


18.97 


24.02 


Tua 


ATua Tliqd 


Tliqd2 Tvapr Psat 


T sump 




12.80 


22.24 2.28 


2.23 2.22 -1.73 


-1 1 . 


,4 




Thetab 


Htube Qdp Athetab Ahtube 


AuQdp 




10.564 


1 . 337E+03 1. 


.412E+04 20.002 5 . 1 98E+02 


1 . 040E+04 


Oata Set 


Number * 10 


8ulk Oil X - 0.0 








TIME: 


18:13:45 










TC No: 


1 2 


3 4 5 


6 


7 


8 


Temp : 


12.22 13.34 


13.25 12.95 -99.99 


12.41 


12.41 


13.25 


TC No: 


9 10 


11 12 13 


14 


IS 


16 


Temp : 


16.82 27.34 


22.25 26.62 10.12 


27.93 


18.43 


23.31 



Tua ATua Tliqd Tliqd2 Tvapr Psat Tsump 

12.75 21.60 2.27 2.18 2.20 -1.75 -11.4 

Thetab Htube Qdp Athetab Ahtube AuQdp 

10.536 1.340E+03 1.412E+04 19.392 5.143E+02 9.973E+03 



Oata Set 
TIME: 

TC No: 


Number * 1 1 

18:14:39 
1 2 


8uik Oil X - 0.0 

3 4 5 


6 7 


8 


Temp : 


12.33 


13.31 


13.25 12.95 -99.99 


12.39 12.36 


13.24 


TC No: 


9 


10 


11 12 13 


14 IS 


16 


T emp : 


16.97 


27.94 


22.33 26.72 10.10 


28.04 13.53 


23.37 


T wa 


A Twa 


Tliqd 


Tiiqd2 Tvapr Psat 


T sump 




12.75 


21 .67 


2.23 


2.20 2.19 -1.76 


-11.4 





Thetab Htube Qdp Athetab Ahtube AuQdp 

10.547 1.340E+03 1.413E+04 19.469 5.I61E+02 1 .005E+04 

Oata Set Number - 12 8ulk Oil 7. » 0.0 

TIME: 18:15:27 



TC No: 


1 


n 


3 4 5 6 7 


9 


T emp : 


12.28 


13.35 


13.22 13.00 -99.99 12.49 12.49 


13.19 


TC No: 


9 


10 


11 12 13 14 15 


16 


T emp : 


16.63 


27.45 


21.92 26.21 9.39 27.50 18.16 


22.96 


Twa 


A T wa 


Tliqd 


Tliqd2 Tvapr Psat Tsump 
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12.79 21.27 2.29 2.2-1 2.23 -1.72 -11.4 

Thetab Htube Qdp Athetab Ahtube AuQdp 

10.545 I.334E+03 1.406E+04 19.029 S.060E+02 9.628E+03 



Data Set 


Number a 1 3 


8u 1 k 


Oil Z • 0.0 








TIME: 


18:20:06 












TC No: 


1 2 


3 


4 5 


6 


7 


8 


Temp : 


9.58 10.38 


10.34 


10.11 -99.99 


9.65 


9.65 


10.23 


TC No: 


9 10 


1 1 


12 13 


1 4 


15 


16 


Temp : 


16.51 27.11 


21.70 


25.94 9.92 


27.19 


18.05 


22.71 


Twa 


ATwa Tliqd 


Tliqd2 Tvapr Paat 


T sump 




9.93 


21.08 2.25 


2.19 2.21 -1.75 


-1 1 .5 






Thet ab 


Htube Qdp 


Athetab Ahtube 


AuQdp 




7.718 


1 . 297E+03 1. 


.00IE+04 13.862 S.090E+02 


9.60IE+03 


Oata Set 


Number * 1 4 


8u lk 


Oil Z - 0.0 








TIME: 


18:20:40 












TC No: 


1 2 


3 


4 5 


6 


7 


8 


T emp : 


9.64 10.30 


10.21 


9.96 -99.99 


9.60 


9.57 


10.20 


TC No: 


9 10 


1 1 


12 13 


1 4 


15 


16 


Temp : 


16.49 27.08 


21 .68 


25.90 9.91 


27.16 


18.04 


22.69 



Twa ATwa Tliqd Tllqd2 Tvapr Paat Taump 

9.87 21.06 2.22 2.17 2.16 -1.79 -11.5 

Thetab Htube Qdp Athetab Ahtube AuQdp 

7.68S 1 .3036+03 1.001E+04 18.874 S.074E+02 9.S77E+03 



Oata Set 
TIME: 

TC No: 


Number 

18:25 

1 


- IS 
: 06 
2 


8ulk Oil X • 0.0 

3 4 5 


6 


7 


8 


Temp : 


6.06 


6.30 


6 . 25 6.16 -99.99 


6.00 _ 5 . 99 


6.27 


TC No: 


9 


10 


11 12 13 


1 4 


15 


16 


Temp : 


16.55 


27.13 


21.72 25.96 9.99 


27.21 


13.11 


22.75 


Twa 

6.12 


ATwa 

21.11 


T liqd 
2.24 


Tliqd2 Tvapr Paat 
2.20 2.23 -1.74 


T aump 
-1 1 .6 





Thetab Htube Qdp Athetab Ahtube AuQdp 



3.893 1 . 267E+03 4.931E+03 18.889 5.062E+02 9.S6IE+03 



Oata Set 


Number 


- 16 


8ulk Oil Z - 0.0 






TIME: 


18:25 


: 4 1 








TC No: 


1 


2 


3 4 5 


6 7 


8 


Temp : 


6.02 


6.37 


6.33 6.22 -99.99 


5.98 5.99 


6.32 


TC No: 


9 


10 


11 12 13 


14 15 


16 


T emp : 


16.58 


27.16 


21.75 25.99 10.01 


27.24 18.15 


22.77 


Twa 


ATwa 


Tliqd 


Tliqd2 Tvapr Paat 


T aump 




6.15 


21.14 


2.28 


2.23 2.27 -1.69 


-11.6 




Thet ab 


Htube 


Qdp Athetab Ahtube AuQdp 




3.881 


I.27IE+03 4. 


933E+03 18.876 S.069E+02 9.567E+03 



Oata Set 


Number 


- 17 


8ulk Oil 7. - 0.0 






TIME: 


18:30 


: 07 








TC No: 


1 


n 


3 4 5 


6 7 


8 


Temp : 


cn 


4.65 


4.63 4.56 -99.99 


4.42 4.43 


4.64 


TC No: 


9 


10 


11 12 13 


14 15 


16 


Temp : 


16.75 


27.50 


22.00 26.31 10.14 


27.59 18.34 


23.05 


Twa 


ATwa 


T l lqd 


TLiqd2 Tvapr Paat 


T aump 




4.53 


21 .39 


2.30 


2.24 2.27 -1.69 


-11.7 




Thetab 


Htube 


Qdp Athetab Ahtube AuQdp 




2.255 


1.258E+03 2. 


, 359E+03 19.121 5.090E+02 9.732E 


*03 



Oat a Set Number - 18 8u lk Oil *1 51 0.0 
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TIME: 

TC No: 
Temp : 

TC No: 
Temp : 
Tua 
4.54 
Thetab 
2.262 

Oata Set 
TIME: 

TC No: 
Temp : 

TC No: 
Temp : 
Tua 
3.40 
Thetab 
1.176 

Oata Set 
TIME: 

TC No: 
Temp : 

TC No.- 
Temp : 
Tua 
3.43 
Thetab 
I .165 

Oata Set 
TIME: 

TC No: 
Temp : 

TC No: 
Temp : 
Tua 
2.67 
Thetab 
.444 

Oata Set 
TIME: 

TC No.- 
Temp : 

TC No: 
Temp : 
Tua 
2 . 70 
Thetab 
.459 
NOTE: 22 



18:30:52 



1 


n 


3 4 5 


6 7 


3 


4 . 48 


4.65 


4.64 4.58 -99.99 


4.44 4.44 


4.64 


9 


10 


11 12 13 


14 15 


16 


16.75 


27.48 


21.98 26.30 10.11 


27.57 18,33 


23.05 


A Tua 


T 1 iqd 


Tliqd2 Tvapr Psat 


T sump 




21.38 


2.30 


2.24 2.28 -1.68 


-1 1 .7 





Htube Qdp Athetab Ahtube AuQdp 

I.259E+03 2.847E+03 19.104 S.067E+02 9.681E+03 

Number • 19 Bulk Oil 7 » 0.0 

18:35:00 



1 


*■* 

L. 


3 4 5 


6 


7 


8 


3.39 


3.46 


3.46 3.42 -99.99 


3.33 


3.34 


3.46 


9 


10 


11 12 13 


14 


15 


16 


16.68 


27.38 


21.91 26.19 10.12 


27.45 


18.27 


22.93 


ATua 


T1 iqd 


Tliqd2 Tvapr Psat 


Tsump 




21 .30 


2.25 


2.20 2.22 -1.74 


-1 1 . 


8 




Htube 


Qdp Athetab Ahtube 


AuQdp 




1 . 306E+03 1, 


. S36E+03 19.079 S.060E+02 


9.654E+03 


Number 


- 20 


8ulk Oil X - 0.0 








18:35 


: 42 










1 


2 


3 4 S 


6 


7 


8 


3.42 


3.48 


3.49 3.46 -99.99 


3.37 


3.38 


3.50 


9 


10 


11 12 13 


14 


15 


16 


16.70 


27.38 


21.91 26.19 10.09 


27.46 


18.27 


22.95 


ATua 


Tl iqd 


Tliqd2 Tvapr Psat 


Tsump 




21 .30 


2.30 


2.23 2.27 -1.69 


-11. 


8 





Htube 


Qdp Athetab Ahtube AuQdp 




1.31 9E*03 1. 


. S37E+03 19.034 S.083E+02 9.67SE+03 


Number 


* 2 1 


Sulk Oil 7. - 0.0 






18:40 


: 25 








1 


2 


3 4 5 


6 7 


8 


2.68 


2.68 


2.69 2.67 -99.99 


2.64 2.65 


2.71 


9 


10 


11 . 12 13 


14 15 


16 


16.66 


27.28 


21.83 26.10 10.16 


27.37 18.22 


22.88 


ATua 


T 1 iqd 


Tliqd2 Tvapr Psat 


T sump 




21 .25 


2;23 


2.21 2.23 -1.73 


-11.9 





Htube Qdp Athetab Ahtube AuQdp 

I . 40 4E+03 6 . 232E+02 19.020 S.044E+02 9.S93E+03 



Number 

18:4! 

1 


- 22 
: 07 


Bulk Oil 7. * 0.0 

3 4 5 


6 7 


3 


2.71 


2.71 


2.73 2.70 -99.99 


2.67 2.67 


2.73 


9 


10 


11 12 13 


14 15 


16 


16.67 


27.30 


21.84 26.11 10.07 


27.38 18.22 


22.90 


ATua 


T 1 iqd 


Tliqd2 Tvapr Psat 


T sump 




21.25 


2.29 


2.22 2.26 -1.70 


-11.9 





Htube Qdp Athetab Ahtube AuQdp 

1.41 5E+03 6 . 2 1 6E+02 18.985 S.03IE+02 9.SSIE+03 

data runs were stored in file OAT040SQ42 
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APPENDIX D. PROGRAM SETUP72 



Program SETUP72 is listed on the following pages. The program enables the user 
to: 

1. Monitor coolant sump temperature. 

2. Monitor liquid pool average temperature. 

3. Monitor all thermocouple channel output temperatures. 

4. Monitor voltage, current, and resulting power supplied to the upper tube heater 
(main heater) as well as the resistance of the heater. 

5. Monitor the voltage, current, and resulting power supplied to the lower tube heater 
(auxiliary heater) as well as the resistance of the heater. 

Program SETUP72 is written in Hewlett-Packard Basic 5.0 for the HP 9300 series 
computer. 
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; i 
2 i 
4 i 
c t 
! 0 
:0 
j t 



PROGRAM: SETUP 
n^T£ : AUGUST 3 ! 99 1 

PROGRAMMER: LT OE.AN SUGIYAMA 
MODIFIED BY LAMNIE LAKE JAM 30, 1992 
COM /Cc/ C(7) 

DATA 3 . ! 0095931 ,25727 . 34359 .-757345 . 9295 .79925595 . 9 ! ,-9247495589 ,5 . 976595^ 
.551925+13 
DATA 3.94075E+U 
READ C( • ? 

ON KEY 1 .15 SOTO 27 
PRINTER IS ! 

PRINT 



PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
• PRINT 
PRINT 
PRINT 
PRINT 
BEEP 
INPUT 
IF Ido 
IF Ido 
IF Ido 
IF Ido 



USING " 4 X , " " 
USING * 5 X , " " 
USING "SX."" 
USING "EXC- 
USING "6X , " " 
USING "5X , " " 
USING EXC- 
USING "4X , " " 

Ido 

;S THEN Ido- 
-0 THEN 50 
! THEN 155 
-2 THEN 1.73 



SELECT OPTION""" 
0-MGNITOR SUMP""" 

! "MONITOR LIQUID""* 

2- CHECK THERMOCOUPLES" 

3- CHECK MAIN HEATER""" 

4- CHECK AUX HEATERS""" 

5- EXIT PROGRAM""" 

NOTE: KEY I - ESCAPE"" 



IF Ido-3 THEN 195 
IF Ido-4 THEN 195 
IF Ido-5 THEN 231 
PRINT 



PRINT 

PRINT "SUMP TEMPERATURE OEG C " 
PRINT 

OUTPUT 709; "AR AF ! 9 All 9 URS" 
OUTPUT 709; "AS SA" 

Sun-0 

FOR J-l TO 5 
ENTER 709 ;£ 

Sun-Sun+E 
NEXT J 
Eave-Sun/5 
Tonp-FNTvsv( Save ) 

PRINT USING "4X ,MCO , 00 " ; T-rp 
BEEP 
PRINT 
UA I T 5 
GOTO 59 

PRINT 

PRINT "LIQUID TEMPERATURE OEG 0" 
PRINT 

OUTPUT 799 s" AR AF 1 G AL17 URS" 
Sun-0 

FOR I -I TO 2 
OUTPUT 709;" AS SA " 

’S3 ENTER **09 ; E 
! S 4 Sun-Sun A E 
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Men r 

Eava-Su.-s/O 
T 2 r*.p * FN T v s v ( E a v a ) 

PRINT USING " 4 X MOO . 00 " ; Tar.p 

9EEP 

WAIT 5 

GOTO !S5 



PRINT 

PRINT "CHANNEL TEMPERATURE OEG C" 

OUTPUT 709 ; ”AR AF00 AL!9 UR 5" 

FOR 1-1 TO 23 
OUTPUT 739: "AS SA" 

Sun-G 

FOR J - ! TO 5 
ENTER 739 ;E 
Sun-Sun+E 
NEXT J 
Eave-Sun/5 
Tenp-FNTvav< Eava ) 

PRINT TA0(3):IiTA0( ISJjTenp 
NEXT I 
9EEP 
WAIT S 
GOTO 173 
f 

PRINT 

OUTPUT 739; “AR AF23 AL22 URS" 

FOR >! TO 3 
OUTPUT 739; “AS SA" 

Saw-O 

FOR J-l TO 5 
ENTER 739 ;E 
Sun-Sum+E 
NEXT J 

IF 1-1 THEN Uoit-Sum/5 
IF 1-2 AND Ido-3 THEN 

PRINT "MAKE SURE U OUTAGE 90X IS SET TO MAIN HEATERS" 
Amp-Sun/S 
END IF 

IF 1-3 AND Ido-4 THEN 

PRINT “MAKE SURE 'JOLT AGE 90X IS SET TO AUX HEATERS" 
Anp-Suw/S 
ENO IF 
NEXT I 

Arp«A9S( Anp*l .9192 ? 

Uc 1 t-A0S( Uo i t • 25 ) 

Power-Uo 1 1 •Amp 
Ras is t-nca-Uol t /Arcp 
ZZi PRINT 
225 9EEP 

225 PRINT "UOLTAGE'U) CURRENT ( A ) RES ISTENCE( chns ) POUER(U)" 
227 PRINT 

229 PRINT USING - IX ,S< M0000. 00 ,4X )" :Uolt .Anp .Resistance ,Pouer 
229 UAIT 5 

233 GOTO 195 

231 9EEP 

232 P°INT 

333 PRINT "THAT’S ALL FOLKS'" 

234 ENO 
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ut u> r* o) 01 e> — rj c> 
i ) i') hi hi r<) tj ^ in 
ririririM M ri r*i r i 



ncc FNTvsv'U > 

COM / Cc / C \ 7 ) 

T-C<0? 

FOR !-! TO 7 
T-T+C( I >*U"I 
NEXT I 

T-T+9.62S897E-2 A T*{ 3.76 ! ! 99E-3-T*S . 0ES9C59E-5 ) 

RETURN T 

FNEND 
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APPENDIX E. PROGRAM DRP72 



The data acquisition and reduction program DRP72 is listed on the following pages. 
Program DRP72 is written in Hewlett-Packard Basic 5.0 for the HP 9300 series com- 
puter. 
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!0 1 f T’_E MflMC ; 2PP72 



CO ' DATE: MARCH 5 ; ca~* 



C3 * RE'. 1 1 9 ED VERSION GF DP.P71 FOR TWO TUBE DATA 
‘IS i REUISEO BY Larrtie Like 
SS COM /Ido/ Ido 

S0 PRINTER IS ! 

70 CALL Select 

53 INPUT 11 'JAN T TO SELECT ANOTHER OPTION ( I »Y ,0-N )?" . I sei 

90 if THEM GOTO 79 

130 BEEP 

! I 0 9EEP 

!C0 PRINTER IS ! 

130 PRINT "DATA COLLECTION/REPROCESSING COMPLETED" 

HQ END 

150 SUB Main 

150 COM /Ido/ Ido 

170 COM /Cc/ C( 7) .leal 

190 COM /Wil/ 02 ,0i ,Do ,L ,Lu ,Kcu 

190 OIM Emf ( 20> ,T< 20> ,Ola( 13) ,02a( 13 ) ,Oia< 1 3 > ,0oa< !3\La( 13 ) ,lua< 13? ,Kcua< 13) , 
Et ( ! 9 ) ,TnS( 4 ) C 1 S 3 

200 DATA 0.10095091 ,25727.94369,-757345.9295,79025595.9! 

210 DATA -924748S599 ,5 . 975995+1 1 ,-2 .55 1 92E+! 3 ,3. 94079EH 4 
220 READ CC • ) 

230 ! DATA " Smooth" , "High F lux " , - Thermoexei-E" , " Thermoexel-HE" 

240 DATA Smooth, High Flux , Turbo-0 .High Flux Mod,Turho-9 Mod 
250 READ Tn3( * > 

259 PRINTER IS 70 1 
279 SEEP 

290 IF Idp-4 THEN PRINTER IS 1 
290 IF I dp-4 THEN GOTO 3030 



300 i INPUT "ENTER MONTH. DATE AND TIME ( MM : 00 : HH : MM : SS 
310-1 OUTPUT 709; "TO" ?QateS 
320 i OUTPUT 709; "TO" 

330 I ENTER 709;DateS 
340 PRINT 



353 ! PRINT 
350 PRINT 
370 PRINT 
3 S3 PRINT 

390 BEEP 
400 INPUT 
410 PRINT 
420 SEEP 



Month, Date and Time :Date$ 
Date : " ,DATE$( TIMEDATE > 

USING " l 0X , * "NOTE! Program name : 0RP72""" 

"ENTER DISK NUMBER * ,0n 

USING "1SX ."-Disk number - ‘•".ZZ'tOn 



430 INPUT "ENTER INPUT MCOE ( 0-3054A , 1 -F ILE ) " , Im 
440 BEEP 



.OateS 



450 INPUT "1 0R 2 TUBE OPERATION ( ENTER l OR 2>",Humntu- 

acQ dccd 

470 INPUT "SELECT HEATING MODE (0-ELEC; ! -WATER Ihn 
4 50 SEEP 

490 INPUT "ENTER THERMOCOUPLE TYPE ( 0-NEU .! -OLD leal 
500 IF In-0 THEN 
510 SEEP 

520 INPUT "GIVE A NAME FOR THE RAW DATA FILE" ,02_* i leS 

530 PRINT USING "!SX,""New file name: " “ , ! 4A " : 02^? i l e3 

540 S : I -20 

550 CREATE 90AT 02_f i IeS .Size! 

550 ASSIGN OF i L eC TO 02_r*iLeS 
E70 i 

590 * DUMMY FILE UNTIL Mrun KNOWN 
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533 0! f 1 1 “'3* ' DUMMY “ 

520 CREATE SCAT 0 1 _? i i e$ , 3 i :e ! 

512 ASSIGN 2Fi la I TO D 1 ; L-S 
922 OUTPUT ?F 1 1 s ! ; Ga t -S 
533 IF Ihn-0 THEN 
5^0 SEEP 

550 INPUT “ENTER NUMBER OF OEFECTIUE TCS ( 3 -DEFAULT ) " , Idtc 
559 IF Idtc *3 THEM 
573 Ldtc ! *0 
520 Ldtc2-0 

530 PRINT USING "!BX, -- No defective TCs exist"-” 

703 END IF 

710 IF Idtc-! THEN 

720 SEEP 

720 INPUT “ENTER OEFECTIUE TC LOCATION ( 1 -9 > ” ,Ldtc 1 

740 PRINT USING “1SX,'*“TC is defective at location “ " ,00“ ;Ldtc ! 

750 Ldtc2-0 
7E0 END IF 
773 IF Idtc-2 THEN 
790 SEEP 

730 INPUT “ENTER OEFECTIUE TC LOCATIONS ( 1 -9 ) “ ,Ldt c ! ,Ldtc2 

900 PRINT USING “ISX,““TC are defective at locations * - .00 ,4X .00" ?Ldtc 1 ,Ldtc2 

910 ENO IF 

920 IF Idtc'* 2 THEN 

930 BEEP 

940 PRINTER IS I 

950 SEEP 

950 PRINT “INVALID ENTRY- 
970 PRINTER IS 70! 

990 GOTO 540 
590 ENO IF 
900 ENO IF 

919 OUTPUT 9Filel iLdtcl ,Ldtc2 

920 IF Hunntu-1 THEN GOTO ! ! 50 

930 INPUT ’ENTER NUMBER OF OEFECTIUE AUX TCS ( 0-QEFAULT )" ,A idtc 
940 IF Aidtc-0 THEN 
550 Aldtc ! -0 
550 Aldtc2-0 

970 PRINT USING “ISX,““No defective AUX TCs exist- - " 

990 ENO IF 

990 IF A idtc- 1 THEN 

1000 SEEP 

1010 INPUT “ENTER OEFECTIUE TC LOCATION ( 9- 1 6 > ” ,A Idtc I 

1020 PRINT USING -|5X,- - TC is defective at location - * ,00" : Aldtc ! 

1030 Aldtc2-0 
1040 ENO IF 
1350 IF Aidtc-2 THEN 
I0S0 SEEP 

1970 INPUT “ENTER OEFECTIUE TC LOCATIONS ( 5-1 5 )” .Aldtcl ,Aldtc2 

1080 PRINT USING ’’15X.--TC are defective at locations -- ,00 ,4X ,00" jAldtc 1 .Aldtc 

n 

1990 ENO IF 
! 100 IF Aidtc; 2 THEN 
1110 9EEP 
! ! 29 PRINTER IS I 
1 1 30 SEEP 

! ! 40 PRINT “ INVALID ENTRY” 

! 1 50 PP INTER IS 701 
1 1 59 GOTO 540 
M70 END IF 
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! ! 3G OUTPUT 0F i i s ! ; .A i d t c ! Aldtc2 

| ! an i i»|a] 2 p t :cn 

1 200 ci_ zc 

! ? t .0 9EEP 

! 220 INPUT "GIUE THE NAME OF THE EXISTING OATA c ILE " ,02_f : l s$ 

1230 PRINT USING " ! SX , " “Old file name: • " . ! 4A" :D2_f i ia* 

1210 ASSIGN SFileZ TO 02_fite3 
! 250 ENTER *5Fi 1-2 :Nrun 
12S3 ENTER 0F i 1-2 ; Oo l d3 

! 279 PRINT USING - ! SX , “ " This data set taken on : " - , 14A" ? Do Id* 

1 290 ENTER OFi le2 ;Ldtc ! ,Ldtc2 ,A1 tds 1 ,A1 tdc2 
1230 IF Ldtc 1 >0 OR Ldtc2:0 THEN 

1300 PRINT USING * 1 SX , “ " Thernocoup ies were defective at locat ions : - * ,2( 30 ,4;< )** ; 
Ldtc! ,Ldts2 
1310 END IF 

1323 IF Humntu- 1 THEN GOTO 1353 
1330 IF Aldtc! >0 OR Aldtc2>0 THEN 

1313 PRINT USING " 1 GX , " ** AUX Thermocouples were defective at locat ions : - " ,2< 30 ,4 

X >" ; Aldtc I ,Aldtc2 

1350 ENO IF 

1 3S0 ENTER OFi Is2 i It t 

1370 ENO IF 

1390 Idic-0 

1393 .IF Ldtc 1 >3 THEN Idtc-Idtc*! 

1400 IF Ldtc2>0 THEN Idtc-Idtc+1 
1410 IF Humntu-1 THEN GOTO 1440 
1423 IF Aldtc 1 >0 THEN Axdtc-AtdtcH 
1433 IF Aldtc2>0 THEN A idt c*A idtc+l 
14401 IF lm»0 ANO Ihn-I THEN 1S95 
1450 SEEP 

14S0 INPUT "'JANT TO CREATE A PLOT FILE? ( 0-N , 1 -Y > " , Ip lo t 
1470 IF Ipiot-1 THEN 
1490 9EEP 

1490 INPUT "GIVE NAME FOR PLOT FILE** ,?J t le* 

1500 CREATE 90AT P_file*,4 
1513 ASSIGN 3P lot TO P.PileS 
1523 ENO IF 
1533 IF Ihm-t THEN 
1540 SEEP 

1550 INPUT -UANT TO CREATE Uo FILE? (3-N.1-Yr.Iuf 
1 5S0 IF Iuf-1 THEN . 

1570 BEEP 

1580 INPUT -ENTER Uo FILE NAME” ,Uf i le* 

1590 CREATE BOAT UfileS ,4 
1S00 ASSIGN OUftle TO Ufile* 

IS 10 ENO IF 

1 S23 SEEP 

1530 INPUT "UANT TO CREATE Re FILE? ( 0-N r 1 -V > - . Ire 
1540 IF Ire-1 THEN 
1550 k 9EEP 

1SS0 ‘input -ENTER Re FILE NAME" .Ref i la* 

1670 CREATE BOAT Ref lie*. 10 

1S30 ASSIGN ORefile TO Refile* 

! E90 ENO IF 

1730 ENO IF 

1710 °RINTER IS 1 
1723 IF In-0 THEN 
1733 SEEP 

1^40 PRINT USING -4X f ’* -Select tube number""" 

1750 IF Ihrt-0 THEN 
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1753 


PRINT 


US IMG " SX , " 0 


smooth 4 inch Ref 


1 770 


PRINT 


USING " 5 X , “ “ ! 


Smooth 4 inch Cu (Presa/S 


i 730 


PRINT 


USING "SX , " "2 


Soft Solder 4 inch Cu “ “ “ 


! 790 


PRINT 


USING “SX," "3 


Soft Solder 4 inch HIGH F! 


1903 


PRINT 


USING “SX,-" 1 


Uieland Hard 9 inch"*" 


1910 


PR IMT 


USING “ 5 X , “ " 5 


HIGH FLUX 9 inch"*'* 


1920 


PRINT 


USING “G:<,““S 


GEUA-K 40 Fins/ in 


1930 


PRINT 


USING "SX,""7 


GEUA-K 2S Fins / in**"" 


1940 


PRINT 


USING "SX , " “ 9 


GEUA-T 19 Fms/in““" 


! 950 


PRINT 


USING “SX , " "9 


GEUA-T OR GE'JA-TY 26 Fine 


1960 


PRINT 


USING "SX , " “ ! 0 


THEPMGEXCSL-E* " " 


1970 


PRINT 


USING "SX ,"*! 1 


THERMQEXCEL-HE” " " 


1933 


PRINT 


USING -sx."-i: 


TURBO- 8 " " " 


1990 


PRINT 


USING *SX . ** 1 3 


GEUA-K 19 F i n s / i n “ “ “ 


1933 


ELSE 






1 910 


PRINT 


USING "SX ,""0 


Smooth tube"'”* 


1923 


PRINT 


USING "SX .""1 


High Flux""" 


1930 


PRINT 


USING “SX ,**2 


Turbo- 8 " " " 


1943 


PRINT 


USING "SX , * "3 


High Flux Mod""" 


1950 


PRINT 


USING "SX ,"“4 


Turbo-9 Mod" 


1960 


ENO IF 




1373 


INPUT 


I t t 





1390 OUTPUT SFllel j 1 1 1 

1990 ENO IF 

2000 PRINTER IS 701 

2010 IF Itt<10 THEM PRINT USING "lex^-Tube Number: * * ,0* ? 1 1 1 

2020 IF 1 1 i>9 THEM PRINT USING "16X, -- Tube Number: * * .OO* ; 1 1 1 

2030 IF Ihm-1 THEN PRINT USING " 1 GX , " "Tube Type: - " , 1 SA" ; Tn* < 1 1 1 > 

2040 BEEP 

2050 INPUT “ENTER OUTPUT VERSION ( 0-LONG ,! -SHORT ,2*N0NE > * , Iov 
20S0 BEEP 

2070 INPUT “SELECT ( 0-LIQ , 1 »VAP ,2*< LIQK'AP >/2 J " , I Iqv 
20901 

20901 DIMENSIONS FO TESTED TUBES 
21001 ELECTRIC HEATED MOOE 

2110 * QI*Oiameter at thermocouple positions 

2120 DATA .01 1 1 125, .01 1 1 1 25 , .01 1 1 1 25 , .0129540 ,.012446 , .01 29540 , .0100965 
2130 DATA .0100965 , .0! 157 . . 01 1 57 , .01 1 57 , .01 1 57 , .011 57 , .0100965 
2140 READ 01 a( ♦ ) 

2150 Ol-Ola(Itt) 

2160* 

21701 02*0iameter of test section to the base of fins 

2190 OATA .015975,. 01 5975 ,.01 5975,. 01 5924 , . 0 1 5975 , . 0 1 5924 ,.01270 

2190 OATA .01 27, .01 39, .0139,. 0139, .0139, .0138, .01 27 

2200 READ 02a< * > 

2210 * 

2220 * Oi*Inside diameter cf unenhanced ends 

2230 OATA* .0127 . .01 27 , .01 27 . .01 32 . .0127 . .01 32 , .01 I I 125 , .01 1 1125 
2240 OATA .01 19 , .01 19 ,.01 1 9 , .01 19 , .01 19 , .01 1 1 125 
2250 READ Oia< * ) 

22501 

2270* Do-Outside diameter of unenhanced ends 

2230 DATA . 01 5975 ,. 0 1 59^5 , .0 1 5975 ,. 01 592 4 , .015975 ,.015924 , . 0 1 270 , . 0 ! 270 
2230 DATA .31331 ..31331 ,.01331 ,.01331 ,.0159, .3127 
2300 READ Dca( * ) 

2310! 

2320' L -Lar.g t h of enhanced surface 

2 330 DATA . 101S , . 1016, .1015, .1316 , .2332 , .2352 , .2332 . .2332 , . 3032 , .2332 ,.2332 . .20 
32 ,.2032 , .2332 
2340 pc an i_ 3 / * > 
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2350 i 



23501 t_u “Length of unenharced surface at the ends 

2370 DATA .0254 , .0254 , .0254 , .0254 , .0752 , .0752 , .0752 . .0752 , .0752 . .07S2 , , 0752 , .07 
52 .0752 .0762 
2390 A SAG Lua( * ? 

2390! 

24001 Kcu-Thermai Conductivity of tubs 

2410 DATA 399 ,344 ,344 ,45 ,344 ,45 ,344 ,344 ,399 ,399 .399 ,399 ,399 ,344 
2420 READ KcuatO 
2430 IF Ihn-1 THEM 



2440i 

2 4501 

2450' 

2470 

2490 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2550 

2570 

2590 

2599 

2500 

2510 

2620 

2630 

2540 

2550 

2SS0 

2570 

2S90 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2790 

I tv 



Data statements for water heating mode 

DATA 0.015975 ,0.015875,3.9169 ,9 . 3 1 39 ,0 . 0! 69 ,0 ,0 ,0 ,3 ,0 
READ 02a< * ) 

DATA 0.0! 27 ,9.9127 ,9.01 45 ,0.0127 ,0.01 45 ,0 ,0 ,9 .0 ,0 
READ Oia< ♦ > 

DATA 0.015975 ,0.015975,3.0169 ,0.015875 ,0 . 9 1 69 ,0 ,0 ,9 ,0,0 
REAO 0oa< ♦ ) 

OflTA 0.3048 ,0.3048 ,9.3048 ,0.3048 ,0.3048 ,0 ,0 ,0 ,0,3 
REAO ta( ♦ ) 

DATA 0.0254 ,0.0254 .0.02S4 ,0.0254 ,0.0254 ,0 ,0 ,0 ,0 ,9 
REAO Lua< ♦ ) 

OATA 399,45,399.45,399,0,3,0,3,0 
REAO Kcua<*> 

ENO IF 
02-Q2a( Itt > 

Oi *0ta< Itt > 

0c-0oa< Itt) 

L-La< Itt ) 

Lu-Lua< Itt ) 

Kcu-Kcua< Itt) 

Xn-.9 

Fr*.3 

IF Itt-0 THEN Cf-1 . 72E X 9 

IF 1 1 1 > 0 THEN Cf *3 . 7037E+1 0 

A-P I ♦ ( Da ' 2-Q i " 2 )/ 4 

P*PI *00 

IF Ihm-1 THEN 

9EEP 

INPUT "TU8E INITIATION MCOE. ( l -HOT WATER ,2-STEAM ,3-COLO WATER )" , 1 1 tm 
IF 1 1 tm- 1 THEM PRINT USING '* 1 SX , " "Tube Initiate: Hot Water 1 '"" 

IF Itim-2 THEN PRINT USING "16X,""Tube Initiate: Steam""" 

IF I tin-3 THEN PRINT USING " l SX , * ‘'Tube Initiate: Cold Water""" 

INPUT "TEMP/VEL MOOE : ( 0-T-CONST .U-OEC: 1 -T-OEC .V-CONST ? 2-T-INC ,U-CONST P , 



2790 IF I tv-0 THEN PRINT USING " 1 SX , " "Temp/Uel Mode: Cons tant /Oecreas mg 

2900 IF Itv-t THEN PRINT USING " l SX ,”T emp/Uel Mode: Oecreae mg/Constant " " " 

2910 IF I tv-2 THEM PRINT USING " 1 6X , " "Temp/Uel Mode: Increas i.ng/Cons tant " " " 

2920 INPUT "WANT TO RUN WILSON PLOT? < 1 -Y ,0-N ) " , Iu 1 1 

2930 IF Ihn-l ANO Iuil-0 THEM 

2940 IF Itt-0 THEN C1-.032 

2950 IF Itt-t OR Itt-3 THEN C1-.059 
2950 IF Itt-2 OR Itt-4 THEN C:-.062 
BEEP 

2990 INPUT '‘ENTER Cl 1 QEF : WH= . 032 ,HF- . 059 , T9- . 062 ) " ,C i 

2990 PRINT USING - ! SX , "" S i sdsr-Tat e 

2900 PRINT USING - 1 5X , ‘ " Constant - " " , l . 40 " ; C i 
2910 ENO IF 

2920 ENO IF 
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2939 IF [ hn= ! AND I.r-! AN 0 luil-l THEN 
29=10 IF Itt-G THEN Ci». 032 
2950 IF ltt-1 OR Iit-3 THEM Ci-.0S9 
2950 IF Itt-2 OR Itt»4 THEM C:-.062 
2970 ASSIGN 9F i I s2 TO ♦ 

2990 CALL UilsenCCf ,Ci > 

2990 ASSIGN 0Fil-2 TO 02_fii2S 

3000 ENTER ®Fil-2;Nrun f 0cid3 ( Ldis! ,Ldtc2 . £ 1 1 

3010 EMO IF 

3020 Naub -0 

3030 IF Idp-4 THEM Ihw-1 
3040 IF Ihm- ! THEN Naub-9 
3050 Ntc-S 

3060 IF Ihr-0 THEN Ntc-20 
3070 J - ! 

3090 S*-0 

3090 Sy-0 
3100 Sx2-0 

3110 Sxy-0 
3120 Rapeai: ! 

3130 IF 'In-0 THEN 
3U0 Ot Id-2.22 
3150 I dc*2 

3150 ON KEY 1,15 RECOVER 3120 
3170 PRINTER IS 1 

3190 PRINT USING "4X "SELECT OPTION*"" 

3190 PRINT USING "GX "0-TAKE DATA""" 

3200 IF Ihn-0 THEN PRINT USING "SX, ""1-SET HEAT FLUX""" 

3210 IF Ihm-1 THEM PRINT USING "SX,""1-SET WATER FLOW RATE"" 

3220 PRINT USING "SX,""2-SET Taat*"" 

3230 PRINT USING "SX , " "3-SET AUX HEAT FLUX""" 

3240 PRINT USING "4X,""NOTE: KEY 1 - ESCAPE""" 

3250 BEEP 

3250 INPUT Ido 

3270 IF I dc> 3 THEM Ido-3 

3290 IF Ido-0 THEN 5290 

3290! 

33001 LOOP TO SET HEAT FLUX OR FLOWMETER SETTIMG 

3310 IF Ido-1 THEN 

3320 IF Ihm-0 THEM 

3330 OUTPUT 709; "AR AF20 AL21 UPS" 

3340 GEEP 

3350 INPUT "ENTER DESIRED Qdp" ,Oqdp 

3350 PRINT USING "4X , " "DESIRED Qdp ACTUAL Qdp""" 

3370 Err- 1000 

3390 FOR 1-1 TO 2 

3390 OUTPUT 709 : "AS SA" 

3400 Sun-0 
3410 FOR it- 1 TO 5 
3420 ENTER 7G9:E 
3430 Sum-S un +5 
3440 NEXT Ji 

3450 IF 1-1 THEN Volt -San/ 5 
3460 IF 1-2 THEN Ano-Sun/5 
3470 NEXT I 

34 90 Anp-ABS ( Arc ♦ I .3132 ? 

3490 Uolt-A0S(Uoit-25 ' 

3590 Aqdc-Uc 1 t •Arp /(PI*D2*L 
3510 IF ASS 1 Aqdc -Qqdo > E'"” r HEL‘ 

3520 IF Aodp.Q"d" 7 HEM 
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3530 8£EP 4000 ,.3 
3540 BEEP 40G0 .3 
3550 BEEP 4000,.: 

3550 ELSE 
3570 BEEP 350 .3 
3590 SEEP 350,. 3 
3590 BEEP 350 \.2 
3500 ENO IF 

3510 PRINT USING " 4X . 3DE ,3X ,NZ . 30E " ;Oqdq ,Aqdp 

3530 WAIT 2 

3530 GOTO 3390 

3540 ELSE 

3550 BEEP 

3550 PRINT USING " 4X ,MZ . 3CE ,3X ,MZ . 3DE" ;Oqdp ,Aqdp 

3570 Err-500 

3530 • WAIT 3 

3550 SOTO 3390 

3700 END IF 

37!0 ELSE 

3720 BEEP 

3730 INPUT "ENTER FLOWMETER SETTING" ,Fma 

3740 GOTO 3190 

3750 END IF 

3750 END IF 

37701 

3790! LOOP TO SET Taat 
3790 IF Ido-2 THEN 
3900 IF Ikdt-I THEN 3850 
38 ! 0 BEEP 

3920 INPUT "ENTER DESIRED Tsaf“,0tld 

3930! PRINT USING "4X,"" OTaat ATaat Rat- Tv Rats""" 

3540 Ikdt-I 
3950 0 1 d 1 -0 

3950 01 d2-0 

3970 Nn- ! 

3990 Nra-Nn MOO IS 
3590 Nn-Nn+t 
3900 IF Np a — ! THEN 

3910 IF Ihn-0 THEN PRINT USING "4X,"" Taat TId! Tld2 Tv 

3920 IF Ihn-I THEN PRINT USING "4X,"" Taat TId! Tld2 Tv Taunp 

t Tqile Tout"" 

3930 END IF 

3940 IF Ihn-0 THEN OUTPUT 709:"AR AFlS ALI9 UR 5" 

3950 IF Ihn-I THEN OUTPUT 739; "AR AF0 AL5 UR 5 " 

39S0 FOR I-l TO 5 

3970 IF Ihn-0 AND I>4 THEN 4240 

3990 5un-0 

3999 OUTPUT 709; "AS SA" 

4000 FOR Ji-I TO 20 

4010 ENTER 709; E i iq 
4030 Sun-Sun+EIiq 
4030 NEXT J i 

4040 El iq-Sun/20 
4050 T 1 d- c NTva v( £1 iq ) 

4050 IF I - 1 THEN TIdl -TId 

4070 IF 1-2 THEN Tld3 s,r ld 

4090 IF 1-5 THEN T.'-Tld 

4090 IF 1-4 THEN Tsunq-TId 

4100 IF 1-5 THEM Tmlst-Tld 



T aunQ 
T tnle 
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4 ! ! 3 IF 1-5 tucm Tcut- r L d 

a i 20 next i 

4 ! 39 IF Ihm-1 THEM 

4 1-10 OUTPUT 709; \AR AF00 AL00 VP.S" 

4 150 OUTPUT 709; "AS SA" 

4 I EG 3um-0 

4170 FOR Kk-1 TO 70 

4190 ENTER 709:5 

4 I 90 Sun-Sum+E 

4200 NEXT Kk 

4219 Emf'7 >-A9S( Sun/29 ? 

4229 Tp iie-Emf< 7)/ 3. 9GE- 4 
1220 ENO IP 

4240 Atld-< Tidl J ’TLd2 >*.S 

1250 IF A8S( At 1 d-Ot Id )> . 2 THEN 

4250 IF At I d>Ot Id THEN 

4270 2GEP 4000,. 2 

4290 BEEP 4000,. 2 

4290 9EEP 4000,. 2 

4300 ELSE 

4310 BEEP 250,. 2 

4329 BEEP 250,. 2 

4330 BEEP 250 1.2 
4340 ENO IF 

4350 Errl-Atid-Oldl 
4360 Oldl-Atid 
4370 Err2*Tv-Qid2 

4399 OId2»Tv 

4390 IF Tidl > 1 00 . THEN 4440 

4400 IF Ihn-0 THEM PRINT USING "4X ,5( MOOO .00 ,2X V ;Qt id ,Tidl ,TId2 ,Tv ,T sump 

4410 IF Ihm-1 AMO Idp-0 THEN POINT USING "4X ,7( MOO. 00 ,2X )" ;0t Id ,TIdl ,Tld2 ,Tv ,Ts 
unp ,T m let , Tp t le 

4420 IF Ihm-1 ANO Idp-4 THEN PRINT USING “ 4X ,5( MOO .00 ,2X ) ,3( M30.00 ,2X ;0t id , Ti 

d ! , T Ld2 , Tv , T sump , Tin 1 et ,Tpt!s,Tout 

4420 UAIT 2 

4440 GOTO 2990 

4450 ELSE 

44G0 IF A0S< At Id-Ot I d > > . 1 THEN 

4470 IF At id>Ot id THEN 

4420 BEEP 3000,. 2 

4490 9EEP 3000,. 2 

4500 ELSE 

4510 9EEP 900,. 2 

4520 BEEP 900 ,.2 

4530 END IF 

4540 Errl-At Id-OLdl 

4550 Old 1 -At Id 

4550 Err2»Tv-0!d2 

4570 0ld2-Tv 

1530 IF Ihm-0 THEM PRINT USING " 4X ,S( MOOO .00 ,2X ? " :0 1 Id ,T Id! ,T Ld2 , Tv ,Tsump 
4590 IF Ihm-1 THEM PRINT USING "4X ,5( NOD .00 ,2X > ,3( M30 .00 ,lX?";Otld,Tldl ,TId2 ,Tv 
, T sump , Tin let ,TpiL- .Tout 
1500 UAIT 2 
4610 GOTO 3990 
1S20 ELSE 
1620 SEEP 

1540 Err ! = A t L d-0 Id 1 

4650 0 1 d 1 -At Id 

4560 *Er’"2 - f .'-Q L d2 
1670 0 1 d2-T .• 
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4630 IF I hr 3 *# THEN PR I .NT USING " 4X ,5( MOOD .DO , 2X V* ; Q t id , T 1 d ! , T Id2 ,Tv , T *unp 

4 c 00 tc tu£n oqt^jt USING -ay g'MQO.QC ZX)’ 1 ;Otid Tld! TldZ Tv Tsunp Tirslst 

, Tp: is , Tcut 

4‘09 WAIT Z 

ir 1 0 GOTO 3999 

*1739 END IF 

4730 END IF 

4740 END IF 

4759 1 

4750! LOOP TO SET AUX HEAT FLUX 

4779 IF Ido-3 THEN 

4790 IF Ihn-0 THEN 

4790 PF 3 . 1 N T " SET VOLT SOX TO AUX " 

4900 OUTPUT 709: "AR AF20 AL22 VR5 M 
49! 9 SEEP 

49Z9 INPUT “ENTER OESIREO AuQdp" .Ouxqdp 

4930 PRINT USING "2X "DESIRED AuxQdp ACTUAL AuxQdp 

4040 Err- 1000 

4950 FOR I- I TO 3 

4850 OUTPUT 709? "AS SA" 

4970 Sun-0 
4990 FOR Jt-1 TO 5 

4990 ENTER 709 ?E 
4900 Sun-Sun+E 
49 ! 9 NEXT Ji 

4920 IF I-l THEN Uolt-Sun/5 
4930 IF 1-3 THEN Amp-Sun/S 
4940 NEXT I 

4950 Anp-A9S( Ar?p* 1 .9192 > 

4950 Vo I t-A9S< Vo 1 1 *25 ) 

4970 Auxqdp-Vol t *Anp/( PI •OZ'L > 

4990 IF ASS< Auxqdp-Ouxqdp >>Err THEM 

4990 IF Auxqdp>Ouxqdp THEN 

5000 9EEP 4000,. 2 

5010 SEEP 4000, .2 

5020 BEEP 4000,. 2 

5030 ELSE 

5040 9EEP 250,. 2 

5050 9EEP 250,. 2 

5050 SEEP 250 . . 2 

5070 ENO IF 

5090 PRINT USING " 4X , M2 . 30E ,2X ,MZ . 30E" ? Ouxqdp .Auxqdp 

5090 WAIT 2 

5100 GOTO 4850 

5110 ELSE 

5120 SEEP 

5130 PRINT USING “ 4X . NZ . 30E ,2X ,MZ . 30E" ; Ouxqdp , Auxqdp 

5140 Err-500 

5150 UA I T 2 

5150 GOTO 4850 

5170 END IF 

5190 GOTO 3190 

5190 ENO IF 

5200 END IF 

5210' ERROR TRAP FOR Ido OUT OF BOUNDS 

5220 IF Ido.’3 THEN 

5230 SEEP 

5240 GOTO 3!°0 

5250 ENO IF 

52E0 1 



91 



9270' TAKE DATA IF In-0 LOOP 
9230 IF I k a 1 - 1 THEN 9320 
9290 9EEP 

9300 INPUT M ENTER SULK OIL 7." ,Bcp 
93 i 0 Ike I* I 

9323 IF Ihn-3 THEM OUTPUT 709; “AR AF00 AL ! 9 UPS" 

9330 IF Ihn* ! THEN OUTPUT 739 ; "AR AF3 AL9 URS" 

9340 IF Ihn«3 THEN Ntc-29 
9390 FOR I-l TO Ntc 
9350 OUTPUT 70S; "AS SA 1 * 

9370 San-3 
5390 FOR Ji-I TO 20 
5390 ENTER 70S;E 
5400 Sun-Sun+E 

5410 IF I*( I 7-Msub ) OR I-(I9-Nsub> THEN Et<Ji-t!-E 
5420 NEXT Ji 
5430 Kdl*0 

5440 IF I -( 1 7-Nsub > OR I»'l 8 -N 3 ub) THEN 
S450 Eave-Sun/20 
5460 Sun-0. 

5470 FOR Jk-0 TO 19 

9490 IF ABS<Et(Jk>-EaveX5.0E-S THEN 

5490 Sun»San+Et ( Jk > 

5900 ELSE 
5510 Kdl-Kdi+t 
SS20 ENO IF 
5530 NEXT Jk 

5540 IF I-( 1 7-Nsub > OR I-MB-Nsub) THEN PRINT USING "4X/"KdI - '’.OD-jKdl 
5550 IF Kdi > 1 0 THEN 
SS60 BEEP 
5S70 SEEP 

5S90 PRINT USING "AX," "Too nuch scattering in data - repeat data set"** 
5590 GOTO 3170 
5600 ENO IF 
55 1 0 ENO IF 

5620 Enf ( I )-Sun/ ( 20-Kdl ? 

5630 NEXT I 

5640 IF Ihn-I THEN 

56S0 OUTPUT 709; "AR AF00 AL00 UR 5 ’* 

5660 OUTPUT 709; "AS SA" 

5670 Sun-0 

5680 FOR Kk-I TO 20 

S590 ENTER 709 ;E 

5700 Sun-Sun+E 

5710* NEXT Kk 

5720 En?(7)-ABS(Sun>/20 

5730 ENO IF 

S710 IF Ihn-0 THEN 

5750 Coun- 0 . 

5750 OUTPUT 709; ”AR AF20 AL22 URS" 

5770 FOR I-l TO 3 
S7S0 OUTPUT 709; “AS SA" 

5790 Su.n-0 

5900 FOR Ji-! TO 5 

5910 ENTER 709 ;E 

5920 Sun-Sun+E 

5930 NEXT Ji 

5940 [F Ccun-0. THEM 

cac^ re r a i THEN Ur-San. 5 

59E0 IF 1-2 THEN Ir-Sun/5 
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S3~0 ELSE 

3332 IF I =* ! THEN A.vr-Sun/5 
5923 IF 1-3 THEN Air-Sun/S 
5222 END IF 
59! 3 NEXT I 

39Z2 IF H’-Ttn tu* f THEN 22 TO EG! 3 

5933 IF CcLin-3 THEM 

59d0 PRINT "SHIFT 'JOLT 9GX TO A'JX " 

5253 END IF 

5250 IF Caun-1 THEM GOTO 5000 

5973 INPUT "TAKE AUX READINGS! ! -YES )' 7 * ,Ccsn 

5990 Ccun-Coun+Ccan 

5290 GOTO 5770 

E000 END IF 

5013 ELSE 

5020 IF Ihn-0 THEN ENTER 9FU-2 ;9op ,ToldS ,EnF( •> ,Vr , Ir ,Avr , A 

5030 IF Ihm-t THEN ENTER 0Fi 1-2 ;9op .ToidS ,En? < • ) ,Fns 

5040 ENO IF 
5050! 

5050! CONVERT snf’S TO TEMP , VOLT .CURRENT 

5070 Tua-0 

5090 Atwa-0 

5090 FOR I - 1 TO Rite 

6100 IF Idtc, N 0 THEM 

Sl!Q IF I-Ldtcl OR I-LdtcZ THEN 

6130 T ( I 5—99.39 

5130 SOTO 5300 

6140 EMO IF 

5150 ENO IF 

6150 IF Humntu-l THEN GOTO 5230 

5 l 70 IF Aidtc>0 THEN 

6130 IF I-ALdtc! OR I-Aldtc2 THEN 

5190 T( I )--92 . 99 

6200 GOTO 5300 

6210 ENO IF 

6220 ENO IF 

5230 IF Itt'4 AND Ihn-0 THEN 

6240 IF I>4 ANO I <9 THEN 

6250 T< 1 5—99.99 

5250 GOTO 5300 

5270 EMO IF 

6230 ENO IF 

5290 T ( I >FNTvsv(Enf( I 5 J 

5330 NEXT I 

53 ! 0 IF 1 1 1 < 4 THEN 

5330 FOR I-l TO 4 

5330 IF I-Ldtcl OR I-Ldts2 THEM 

6340 Twa-Tua 

E350 ELSE 

5350 Tua-Tua+T' I ) 

5370 ENO IF 
5390 NEXT I 
5390 Tua-Tua/' 4- Idtc ) 

5400 ELSE 

E4J0 IF Ihn-I THEM 5500 

5420 FOR I-l TO 9 

5430 IF I-Ldtc! OR I=Ldtc2 THEN 

6440 Tvia-T-ja 

54E0 ELSE 

5463 Tfja* Tuja * T { I ) 
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Cji70 END IF 

5490 NEXT I 

5432 Ty-Tua.M 9-Idtc ■' 

5EQ2 IF Hu:nnt(j-l THEN GOTO S530 

5510 FOR I =*9 TO 15 

EGGO IF !»ALdtcl OR 1 3 A l d t 2 THEN 

5539 Atua a Atvia 

5540 ELSE 

5553 Aiwa a A$wa*TI I ) 

5550 ENG IF 

5S79 NEXT l 

5590 Atw-Atua/«9-Aidtc > 

5590 ENO IF 

5500 Tld-T( ! 7-Msub ) 

5510 T ld2=* T ( I9-Naub > 

5530 Tida-' Tld+Tld2 >*.5 
5530 Tv-T< ! 3-Nsub ) 

5540 IF Itt\3 AND !hn-0 THEN 

5550 Tld2—99.99 

5550 T v » ( T ( 10 > + T( 1 1 ? S/2 

5579 ENO IF 

5590 Tsump-T( 20-Nsub ) 

5590 IF Ihn-0 THEN 5730 
6700 Tinist-T( 13-Nsub > 

5710 Tout-H U-Naub J 
5720 IF Xhn-0 THEN 
5730 A«p-A0S( Ir*1 .9192 > 

5740 Voit-A0S<Ur>*25 
5750 Q-Uoit«Anp 

5750 IF Humntu-J THEN GOTO 5900 
5770 Aua«p-A9S( Air«1 .9192 > 

5790 Auvo I t»A9S( Avr >*25 
5790 Auq-Auvo 1 t *Auanp 
5900 ENO IF 

5910 IF Itt-0 AND Ih«-0 THEM 
5920 Kcu-FNKcuITuS 
5930 ELSE 
5940 Kcu-Kcua< Itt } 

5950 ENO IF 
59501 

58701 FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
5999 IF I hn*0 THEN Tu*Tw-Q*LOG( 02/01 >/( 2*PI*Keu*L ) 

5990 IF Hwtnnta-1 THEN GOTO 5910 

5900 IF Ihn-0 THEN Atw-Atw-Auq*L06( 02/0! )/( 2*PI *Kcu*L ) 

5910 IF Ilqv-0 THEN Tsat-Tida 

5920 IF I Iqv* 1 THEN T sat -T v 

5930 IF Iiqv-2 THEN Tsat -< T! da+Tv S • .5 

5940 IF Ihm-1 THEN 

S3E0 Tavg-Tmist 

5950 Grad-37. 9953^ . !04293*Tavg 

5970 Tdrop-A8SI Enf ( 7 1 )*1 .E-S/< 10*Grad > 

5990 Tavgc-T ir. I at -Tdrcp* . 5 

5990 IF A9S( Tavg-Tavgc >: .0! THEN 

7000 Tavg»( Tavg-Tavgc >• .5 

7010 GOTO 5959 

7920 ENO IF 

7939' 

70401 COMPUTE WATER ponpcpyrcc 
7950 IF Ihn-1 ru EM 
"050 Ku-FNKvi'. F avg • 
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7973 f1uu3 3 f*Mf1uuj'. T 3 v g ) 

7090 Cpu*FNCpui T > 

7099 Pr*4*FNPru< Tavg 
7190 Rhcu*FNRhou( Tavg ) 

7 ! ! 0 Tw i *T avg 
7!20l 

7 ! 10 1 Ccnputs MOOT 

7U0 Mdc t *3 . 35575*3"^ ns ♦ f 3 . 5 l SSSE-^-Fns • ( 9 . SZOQSE-S-Fns ♦ ( 1 .23S99E-7-Fna*4.3l99 
E - ! 9 ) ) ) 

71 E0 1 Mdot»Mdot*( l .03E5-T inlet •< 1 . 9ES44E-3-T in 1st *5 . 252E-S > >/ I .0037 
71 E0 Kdt*0 

7170 Q-Mdot ♦Cpu^Tdrop 

7190 L.nt d* T dr op /LOG < (Tini-t-Taat )/(T in l et -Tdrop-Taat ) ) 

7190 Uo-9/(PI •Qc*L*Lmtd ) 

7200 Ry*0o ♦L0 G( Qq/Q i >/'2.«KeuJ 

7210 Tw-Tsat*Fr*Lntd 

7220 Uui-fldot / ( Rhou*PI *0 1 *2/4 ) 

7230 Raw-Rhou^Uu^Oi/Muua 

7240 Hi«Ci*Kui/Qi*Reu\9*Pru*< 1/3. >♦ ( Muwa/FMMuu< Twi >r .14 
72E0 Tuic*Tavg-Q/( PI*Go*L*Ht > 

7250 IF A9S( Tui-Tuic >> .01 THEM 

7270 Twi*< Twi^Twtc )• .5 

7290 SOTO 7240 

7290 ENO IF 

7300 Tut»( Tui+Tuic >*.5 

7310 Ho-l/( 1 /Uo-Qo/'. Oi ♦Hi >-Ru ) 

7320 ENO IF 

7330 ENO IF 

7340 IF Ihn-l THEN 

73S0 Thefcab*Q/( Hc*PI ♦C?c»L ) 

7380 Tu-Taat+Thetab 
7370 ELSE 

7390 Thst3b*Tw-Taa t 

739! IF Hunntu-1 THEN SOTO 7400 

7390 Athetab*Atu/-Taat 

7400 END IF 

7410 IF Thetab/0 THEN 

7420 BEEP 

7430 INPUT "TUALL.TSAT (0-CONTIMUE, !«ENO>-,Iev 
7440 IF Iev-0 THEN GOTO 3130 
7450 IF lev-! THEN 9930 
7460 ENO IF 

7470 IF Hwmntu-l THEN SOTO 7540 
7490 IF A thetab<0 THEN 
7490 9EEP 

7500 INPUT " AUX TUALL<TSAT (0-CONTINUE. 1 -ENO ? * . A i -v 

75 !0 IF Aiev-0 THEM SOTO 3 130 

7520 IF Aiev-I THEM 9930 

7530 ENO IF 

75401 

75501 COMPUTE VARIOUS °ROPEPTIES 
75EG T f i in*»( Tu+T 3at )*.S 
7570 Rho- c MRho( TP i In > 

7590 Mu«FNMu< T f i Lm > 

7590 V-FMK< T f i In ) 

7500 Cp-FNCp( Tf i In 1 
75 ! 3 3eta-FM0eta<TFiln> 

""529 u f g^CMH ?g< T- a t ) 

"539 *1 1 =Mu. Rho 

tc 4 0 ^l n ha 9>/ /l Ph?*C" ■ 
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7553 Pr*Nt. Alpha 
7553 Psat-FNPsatl Taat J 
7570 1 

7530 1 COMPUTE VARIOUS PROPERTIES FOR AUX TU8E 
7590 IF Hunntu-1 r HEM GOTO 7790 
7700 At f iln-C Atu+Tsat >*.5 
7710 Arho-FNRhc< At? ilrt ) 

7720 A-m-FNrtu' Atf tin? 

7733 Ak-FNK< At ? tin ) 

7713 Acp-FNCp< At ? tin) 

7790 Abeta-FMSet a( At f : In > 

7760 Ani-Asnu/Arho 

7770 Aa lpha-Ak/ ( Arho*Acp > 

7790 Apr* Am /Aa ipha 
77901 

79001 COMPUTE NATURAL-CONVECTIVE HEAT -TRANSFER COEFFICIENT 
7910! FOR UNENHANCED ENO( S ) 

7920 Hbar-130 

7930 Fe-( Hbar*F / ( Kcu *A > > ' .5*Lu 
7943 Tanh-FNTanh( Fa > 

7953 Theia*Theiab*Tanh/Fe 

7960 Xx-< 9.9! *0e t a^The tab *Do'3*Tanh/(Fa*Ni ‘Alpha ) r. 166667 
7970 Yy*( !♦' . 559/Pr > * ' 9/ ! S ) ) m (9/27) 

7990 Hbarc-K/Do*< .S*.397*Xx/Yy )'*2 

7990 IF ASS< ( Hbar-Hbarc >/Hbarc >> .001 THEM 

7900 Hbar-< Hbar+Hbarc )* .5 

7910 SDTO 7930 

7920 END IF 

79301 

7940' COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
7950 J FOR UNENHANCED END'S) FOR AUX TUBE 
7950 IF Hunntu-1 THEN GOTO 9090 
7970 Ahbar- 1 93 

7990 Fe-( Ahbar ♦P.M Kcu*A ) >" .S*Lu 
7990 Atanh-FNTanh( F* > 

9000 At heta “At he tab ♦At anh/ Fa 

9013 Axx-( 9.9! *Abet a*Aihetab*0o *3*Atanh/( Fe^Ani ♦Aalpha ) ) A .166657 
9020 Ayy-( lf( .559/Apr)*'(9/!6))"(9/27) 

9030 Ahbarc-Ak/Do*( . 5+ . 397*Axx/Ayy )"2 

9040 IF ABS( ( Ahbar-Ahbarc >/Ahbarc )> .00! THEN 

9350 Ahbar*( Ahbar+Ahbarc )♦ . 5 

9060 GOTO 7980 

9070 END IF 

9390 1 

9093* COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENO< S > 

9 ! 00 Q l Hbar*P ♦Kcu*A ) " . 5‘ The tab ♦ Tanh 
9110 Qc*0-2 *0 1 

9120 As= a I *02 *L 
9130 ' 

9140 'COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENO< S > OF AUX TU8E 

9150 IF Hunntu-1 THEN GOTO 9190 

9150 Aq 1 *< Ahb ar*P • K c u • A ) " . 5* At he tab ♦At anh 

9170 Aqc-Auq-2*Aql 

9 l 90 i 

9130* COMPUTE ACTUAL U EAT FLUX AMO SOILING COEFFICIENT 

3230 Odp-Oc/As 

3213 Htube-Odp/Th-tao 

3223 Cat’' Cp * That at .• u * g • • ( Qdo." Mu*Hfq.'*' .311/(3.9’ *°hc : “ .5) ' 1 . 3. J*Pr I 
3220' 

3210* CQMP'J r E ACTUAL U EAT FLUX AMO SOILING COEFFICIENT C 0R AUX TUBE 
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3 IE 3 IF Huimrstu 31 THEN GOTO 9139 

9259 .^ui5d3“*^' , *c/As 

9179 Ahtuhe a Auqdp/A thetab 

9190 AcaP=*< Acp 'Athetab .'Hfg )/( Auqdp/( Amu'Hfg ) * ( . 9 1 4/ ' 9 . 9 1 *Arho r .5 >" ( 1/3. >'Apr" 1 
. 7 ) 

9199' 

9309 1 RECORD TINE OF DATA TAKING 
9319 IF 1**0 THEN 
9319 OUTPUT 739 ; “ TO" 

9339 ENTER 799;ToldS 

9340 END IF 
93E0 ! 

9360! OUTPUT DATA TO PRINTER 
9379 PRINTER IS 701 
9390 IF Iov-9 THEN 
9390 PRINT 

9400 'PRINT USING * 1 0X , " "Data Set Number - • " ,000 ,IX , • *9ulk Oil X - - ,00 . 0 ,SX , 1 

4A- ; J ,9a p ,Tald$ 

9410 PRINT USING * 1 0X , • "Data Set Number * -* ,000 ,2X ,• "Bulk Oil 7. - * - ,00 . 0 * ; J ,B 
op 

9423 PRINT” TINE: - , TIMES! TIMEOATE) 

9430 IF Ihm-0 THEN 

9440 PRINT USING -10X,-"TC No: 1 2 3 4 S G 7 

B mmm 

9 450 PRINT USING - 1 0X ."“Tenp ,9( IX , MOO. 00 >":T( 1 ) ,T(2 ) ,T(3> .T(4) ,T( 5 ) , T< 6 > . T( 
7 ) T ( 9 ) 

9451 IF Humntu-1 THEN GOTO 9490 

9450 PRINT USING * 1 9X ,- - TC Mo: 9 10 M 12 13 14 15 

15--- 

9470 PRINT USING - 10X , -"Temp : - " ,9( IX , MOO. 00 >* ;T< 3 ) ,T( 1 0 ) ,T( l 1 ) ,T( 1 2 > ,T( l 3 ) ,T< 1 
4 ) ,T< 15) ,T( IS > 

9499 PRINT USING • 1 0X , " " Tua ATwa Tiiqd Tiiqd2 Tvapr Paat T 3 ump' , - M 
94 90 PRINT USING " ! 0X ,3( MOO . 00 , 1 X ) , I X , MOO. 00 , IX ,2( IX , MOO. 00 ) ,2X .MOO.D" ; Tw ,Atw ,T 
Id ,Tld2 ,Tv .Paat .Taump 

9500 PRINT USING - 1 0X . - " Thetab Htube Qdp Athetab Ahtube AuQd 

P 

9519 PRINT USING “ 1 0X .M00.30 , IX ,M2 .30E . IX ,M2.30E . IX .MOO. 30 , IX .MZ.30E . IX ,MZ .505" 
; Thetab , Htube ,0dp .Athetab .Ahtube .Auqdp 

9520 ELSE 

9539 PRINT USING "19X,'*- Fma Uw Taat Ttnl Tdrap Thetab q Uo 

Ho--- 

9540 PRINT USING “ ! 0X . 4 ( 20 .00 . 1 X ) , Z . 30 . 1 X .00 . 00 . 1 X , 3( MZ . 30E , 1 X ) - ; Fma ,Uu , Taa t ,Tl 
niet ,Tdrop .Thetab ,Qdp ,'Jc ,Ho 

9550 END IF 
9550 ENO IF 
9570 IF I o v * 1 THEN* 

9590 IF J*t THEN 

3590 PRINT 

9500 IF Ihm-0 THEN 



9610 

9620 


PRINT 

ELSE 


USING 


- 1 9X , 


PUN No 


OilS 


Taat 


Htube 


Qdp 


Thetab 


9530 


PRINT 


USING 


- 1 0X , 


, - - FMS 


oil: 


T SAT 


HTUBE 


OOP 


THETAB 



9540 ENO IF 
9550 ENO IF 
9550 IF Ehm-0 THEN 

9S~9 PRINT USING " 1 2X .30 ,4X ,00 ,2X . MOO. 00 ,3M X ,MZ . ODE ) " ; J ,9ap . Taat .Htube .Qdp . The 

\ 3 K 

9530 ELSE 

3599 PRINT USING " 1 2X .30 , IX .00 ,2X , MOD. 00 ,3( ! X ,M2 . 30E Fma ,9cc . Taat .Htube ,Qdp , T 
He t 3 t 
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9723 ENQ IF 
£710 ei\jn re 

a?:a if in-o then 
3732 SEEP 

9749 INPUT "OK TO STORE THIS GATA SET ( ! - Y , 0 -N >? “ ,0k 
37 c ,0 ckjq re 

37E2 IF Ok - ! OR Em - 1 THEN J-JH 
3770 IF OK- I ANO In-0 THEN 

3750 IF Ihn-0 THEM OUTPUT OF 1 1 a 1 ; 9 op # To l dS , En r* < ♦ ) , Ur , I r , A vr , A i r 
9790 IF IMn -1 THEN OUTPUT SFilal ;9op .ToldS # Srsf( ♦ > ,Fns 
9900 EMO IF 

9910 IF Iuf-I THEN OUTPUT OUfileiUu.Uo 
9920 IF Ira-1 THEN OUTPUT SRa f i ie ; Fna .Raw 

8920 IF ( In* I OR Ok- 1 > ANO Iplat-I THEN OUTPUT OPiot ;Qdp .Thetab 
9940 IF In-0 THEN 
3950 BEEP 

9950 INPUT "UTLL THERE BE ANOTHER RUN ( I -Y ,0»N )? " , Go_on 
9970 Nrun-J 

9990 IF Go_sn-0 THEN 9930 
9990 IF Go_on<>0 THEN Rapaat 
9900 ELSE 

9910 IF J <Nrun+ 1 THEN Rapaat 
3920 ENQ IF 
9930 IF In-0 THEN 
9940 BEEP 

99S0 PRINT USING “ 1 0X , " "NOTE: "“,ZZ/" data runs u era stcrad in fils ‘*“,10A"jJ 
1 ,Q 2 _?ilei 

9950 ASSIGN OFilal TO ♦ 

9970 OUTPUT 3Fila2;Nrun-l 
9980 ASSIGN OFilel TO Ol_fila* 

9990 ENTER 3Fi lal ;Qata3 .Ldtcl ,Ldtc2 . It t 

9991 IF Hunntu-I THEN GOTO 9000 

9993 ENTER 0Fi ial ; Aidtc I ,Aldfcc2 

9000 OUTPUT ®F i la2 ; OataS ,Ldtc l ,Ldtc2,Itt 

9001 IF Hunntu-1 THEN GOTO 9010 

9003 OUTPUT OF t Ia2 ; Aldtc 1 ,Aldtc2 
9010 FOR 1-1 TO Nrun -1 

3020 IF Ihn-0 THEN 

9030 ENTER 9Fi lal ;9cp .ToldS .Enf ( ♦> ,Ur ,Ir 

9040 IF Hunntu -1 THEN GOTO 9050 

9050 ENTER OFi la l ; Avr .Air 

9050 OUTPUT 0Fiia2i8op .Told*. Enf (*>, Ur ,Ir 

9070 IF Hunntu-l THEN SOTO 9090 

9090 OUTPUT 0Fi le2 ; Avr ,Air 

9090 ELSE 

9100 ENTER 3Fi le l ;9cp .To ld$ ,Enr ( * ) ,Fns 
9110 OUTPUT OFi Ia2;Bop , TcldS ,Enf ( ♦ ) .Fns 
9120 ENO IF 

9130 . NEXT I 

91 40 ASSIGN GFilel TO • 

31E0 PURSE "DUMMY" 

9160 ENQ IF 

9170 SEEP 

9130 PRINT 

9190 IF Iplot-I THEN PRINT USING " 1 OX , " "NO TE : " " , Z I X-Y pairs uera stared i 



plot 


data f i 


la " " 


10A" ; J-! 


,P_f l iaS 


9230 


ASSIGN 


SFilaZ 


ro * 




92 ! 0 


ASSIGN 


2P let 


TO • 




9220 


IF lur- 


l THEM 


ASSIGN 


}tjrile TO ♦ 


9230 


{ C T r a « 


1 THEM 


ASSIGN 


?Ra f : 1 a TO * 
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32 10 CALL Stats 
9C50 SEEP 

SZ50 INPUT “LIKE TO PLOT DATA ' ! ,0«N " ,Ok 

9ZTG IF Ok-*! THEN CALL PI 2 * 

9250 SUBEND 
3Z30 1 

3300 1 CURUE FITS OF F°OPERTY FUNCTIONS 
9310 OEF FNKcutT) 

33 Z 9 1 QFHC COPPER 350 TO 300 K 

9330 Tk-T+273.15 ' C TO K 

9343 K-434-. ! ! 2 • Tk 

9350 RETURN K 

9350 FNENO 

9370 OEF FNMul T ) 

93901 170 TO 350 K CURUE FIT OF VISCOUS ITY 

9390 Tk -T+273 . ! E » C TO K 

9400 Mu-EXPI -4*. 4635 + < 101 ! .47/Tk > )*1 .0E-3 

94 !0 RETURN Mu 

94Z0 FNENO 

9430 OEF FNCp( T > 

94401 ! 90 TO 400 K CURUE FIT OF Cp 
9450 Tk-T^273. 15 >C TO K 

9460 Cp-. 40 198+1 . 650075-3* Tk + 1 . 5 US4E-S ♦Tk * 2-5 . S7953E- 1 0* Tk - 3 

9470 Cp*Cp*!000 

9490 RETURN Cp 

9490 FNENO 

9500 OEF FMRho(T? 

9510 Tk-T+273.IS >C TO K 

9520 X-1 -< 1 .9*Tk/7S3.9S ) »K TO R 

9530 Ro-3S. 32+51 . 1464 I 4»X*( ! / 3 )+lB .4 1 901 S*X + ! 7. 475938 •X*'* .5+1 . ! 19929*X A 

9540 fio-Ro/. 052429 

9550 RETURN Ro 

9550 FNENO 

9570 OEF FNPr(T> 

9590 Pr-FNCp< T >*FNMu( T J/FNK ( T ) 

9590 RETURN Pr 
9500 FNENO 
9510 OEF FNK ( T ) 

95201 T;3S0 K WITH T IN C 
9530 K- .071 -.00025 I *T 
9540 RETURN K 
9550 FNENO 
9650 OEF FNT anh( X ) 

9570 P-EXPfX) 

9590 O-l/P 

9590 Tanh-< P-0 )/(P ± Q I 
9700 RETURN Tanh 
9710 FNENO 
9720 OEF FNTw/'U) 

9730 COM .'Cc/ C( 7 > . Isa l 
9740 T-C<0> 

9750 FOR I— ! TO 7 
9750 T-T^CI I ? *U * I 
9770 NEXT I 
9790 IF leal- I THEN 

9790 T=*f-S.74Z2934E-C a -T*' 9. 0Z77043E-3-T • < -9 . 32599 1 7E-5 ) ) 

9900 ELSE 

9910 T- Tj *- 9. 32999 7E-2 fc T Z . "5 1 99E-3-T*5.3S.35253E-5 ) 

9920 END IF 

9930 oc tjjrnj r 
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9S4C FNEND 

99£C DEF FN5eta( T > 

9960 P.op-FNRhcl T+ . 1 ) 

9370 Rpn-FNPho'7-. ! ? 

9990 Beta— 2/ ( Rop+Rom >• < Rcp-Rom > / .2 

9990 RETURN Bet- 

9900 FNEND 

99 1 0 DEF FNHfg'T ) 

9920 Hfg-1 . 374 l344E**-5-T*( 3. 309436 !E+2+T*1 .2165143 ) 

9930 RETURN Hfg 

9940 FNEND 

9950 DEF FNPsat ( Tc > 

9950! 0 TO 90 deg F CURUE FIT OF Psat 
9970 Tf-1.e*T=+32 

9990 Pa-5 . 945525+T f * ( . 1 5352092 + Tf •< 1 . 49409G3E-3*T f *9 . 6 1 5067 l E-6 ) ) 

9990 Pg-Pa-14.7 

10000 IF Pg; 0 THEN ■ + -PSI5, — in Hg 

10010 Psat-Pg 
10020 ELSE 

10030 Psat-Pg*29 . 92/14.7 
10040 END IF 
10050 RETURN Psat 
10050 FNEND 

10070 DEF FNHsmootht X ,9op ,1 sat ) 

10090 DIM A ( 5 , B ( 5 > , C ( 5 } , D ( 5 > 

10090 DATA .20525,. 25322,. 3 19048,. 55222 ..79909,1 .00258 
10100 OATA .74515, .72992, .73189, .71225, .68472, .64197 
10110 DATA .41092, .17726, .25142, .54906, .81916,1 .0945 
10120 DATA .7-1403, .72913, .72565, .696691 ,.665967, .51999 
10130 REAO A< O ,B( * > ,C( • ) ,D( • > 

10140 IF Bop<6 THEN I-9op " 

10)50 IF Bop-6 THEN 1-4 
10160 IF 9op-!0 THEN 1-5 
10170 IF Isat-1 THEN 
10190 Hs-EXP' A { I >* t 9( I >-LQG< X ) > 

10190 ELSE 

1 0200 Hs-EXP(C( ! >+QU >*L0G< X > > 

10210 END IF 
10220 RETURN Hs 
10230 FNEND 
10240 DEF FNPo ly( X ) 

10250 COM /Cply/ At 10 , 10 ) ,C( 10 ) ,B(S > ,Nop ,Iprnt ,0po ,Ilog ,Ifn,I join.Njoin 
10260 XI-X 
10270 PoIy-B( 0 ) 

10280 FOR 1-1 TO Nop 

10290 IF Ilog-t THEN Xl-LOG(X) 

10300 Pely-PeIy*B( I >-XI" I 
10310 NEXT I 

10320 IF I log-!. THEN Pc 1 y=EXP( Pc 1 v ? 

'0330 RETURN Pcly 
10340 FNEND 
10350 SU9 Pcly 

10360 DIM R( 10 > ,SM0 ) ,Sy( 1 2 • ,S. • 1 2 ' ,X .-. ( 1 00 J , V y ( 1 00 > 

10370 COM /Cply/ A( 10 , 10 > ,C( 10 ! ,9C 5 J ,N ,Iprnt ,0pc ,1 log ,! fn ,1 joir. .Njcm 
10390 COM / X >: y y / Xp( 25 > p( 25 > 

10390 FOP 1-0 TO 4 
10400 S( I .'-0 
10410 NEXT I 
10120 BEEF 

10430 INPUT “SELECT • 0- c I LE , 1 = *■!£ Y90APD ,2® PROGRAM ' “ , I rr 
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10440 I.— Irs*’. 

| 04 CA DCCC 

10460 INPUT "ENTER NUM8EP OP >.-> PAIRS", Np • 

10470 IP Ir- ! THEN 
1 0490 BEEP 

10490 INPUT "ENTER DATA c ILE NAME " t O_f i le* 

! 0500 BEEF 

10510 INPUT "LIKE TO EXCLUDE OATA PAIRS ( 1 -Y ,0-N P " , I ed 
10520 IP led- 1 THEN 
10530 BEEP 

10S40 INPUT "ENTER NUMBER OP PAIRS TO BE EXCLUDED" , Ipe* 
10550 END IP 

10560 ASSIGN ©File TO 0_?ileS 
10570 ELSE 
10590 BEEP 

10E9C INPUT "WANT TO CREATE A DATA FILE ( 1 -Y ,0-N ) 7 - t Yes 
10500 IP Yes- 1 THEN 
10610 BEEP 

10620 INPUT "SIUE A NAME FOR DATA FILE" ,0_f i leS 

10530 CREATE BOAT 0_fileS,5 

10640 ASSIGN ©File TO Defiles 

10550 END IP 

10560 END IF 

10570 BEEP 

10590 INPUT "ENTER THE ORDER OF POLYNOMIAL" ,N 

10690 FOR 1-0 TO N*2 

10700 Sy < I )-0 

10710 Sx< I )-0 

10720 NEXT I 

10730 IF led- 1 AND Irt-l THEN 

10740 FOR 1-1 TO Ipex 

10750 ENTER ©F:le;X,Y 

10760 NEXT I 

10770 END IP 

10790 FOR 1-1 TO Np 

10790 IF In-1 THEN 

10S00 IF Opo-2 THEN ENTER ©Fiie;Y,Y 

10510 IF Opc<2 THEN ENTER ©F:le; v ,X 

10520 IF Opo-1 THEN Y-Y/X 

10530 IF Ilog-1 THEN 

10940 IF Opo-2 THEN Xt-X/Y 

10950 X-LOG( X ) 

109S0 IF Opo-2 THEN Y-LOG(Xt) 

10970 IF Opc<2 THEN Y-L05CY? 

I 0S9O END IF 
10590 END IF 
10900 IF In*2 THEN 
10910 BEEP 

10920 INPUT "ENTER NEXT X-Y PAIR" ,X ,Y 
10930 IF y s «si y u c [• j OUTPUT © F 1 1 e * X Y 
10940 END IF 
10950 IF In .2 THEN 
1 0950 X ( I >-X 
1 09^0 Yy ( I )» Y 
10990 ELSE 
1 0990 X-X" ( I - i ) 

1 1 000 Y» Yp( I - J ; 

11010 END IF 
1 ! 020 0/ 0 >-Y 

I 1 0 J 0 c y ( 0 H s C y / 0)4/ 
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• ! 042 5< t >-:< 

1 J 0C0 Z. f 1 \aC,{ J \ + X 

11093 c 0 R J-1 ro N 
I I 3 72 P ( J > »R ( J - ! > * :< 

I ! 093 Sy ' J ) ■ S y < J J > 

11093 NEXT J 
! ! ! 02 FOR J-2 TO N • 2 
t I M3 3( J )mc i J - ! >*X 
11120 3a( J >-S.<< J >>S< J ) 

1 1 1 20 NEXT J 
1! 140 NEXT I 

MIS® IF Y-s-1 ANC In-2 THEN 
11160 9EEP 

11170 PRINT USING ” t OX ,00 , " " X-Y pat r 3 uars stored in file , 1 0A* ;Np ,0_f i leS 

1 ! 1 93 END IF 

11190 3*<0>-Np 

11200 FOR 1-0 TO N 

11210 C( I )-Sy( I > 

1 1220 FOR J-0 TO N 
t 1220 ACI ,J )-Sx< I + J ) 

1 1240 NEXT J 
I 1 2S0 NEXT I 
! 1250 FOR 1-0 TQ N- 1 
1 1270 CALL Oiv tde( I ) 

! 1290 CALL Subtract ( I* I ) 

11293 NEXT I 
1 1200 9(M)-C<N)/A<N,N J 
l 1210 FOR 1-0 TO N- l 
11320 8' N-l-I l-CIN-t-I ) 

11330 FOR J-0 TO I 

1 1240_ 8 ( M- 1 - 1 )-8( N-1 - I >-A< N-l-I ,N-J ) *9( N- J ) 

I 1 350 NEXT J 

1 13S0 8 ( N- 1 -I >-S< N- 1 -I ) / A ( N - 1 - 1 ,N - 1 - ! > 

1 1370 NEXT I 
1 1 390 1 PR INTER IS 701 
l 1 390 i PRINT S' ♦ ) 

1 1 400 i PRINTER IS 705 
11410 IF Iprnt-3 THEN 

11420 PRINT USING ’ 1 2X , “ "EXPONENT COEFFICIENT 

1 1430 FOR 1-0 TO N 

11440 PRINT USING " 1 5X ,00 ,5X ,MQ . TOE” ? I , 9 < I > 

1 1450 NEXT I 
1 1460 PRINT ” 

11470 PRINT USING " t 2X , " "DATA POINT X Y Y( CALCULATED ) OISCR 

EPAMCY” " " 

11490 FOR 1-1 TO Np 
1 1490 Yc-9(0) 

1 1500 FOR J-1 TO N 
11510 Yc-Yc*9( J ?*Xx( I r J 
1 1520 NEXT J 
1 1530 O-Yy < I )-Yc 

11540 PRINT USING " 1 EX , 20 , IX , U NO . 50E , 1 X > " ; I ,:</ ( I ) , Yy ( I ) , Yc ,0 
11550 NEXT ! 

11553 ENQ IF 
1 1570 ASSIGN ?F i i- TO • 

11590 SU9EM0 
11590 SU9 Oiv tde<MJ 

1 1 503 CON Cpl/.' A( 1 0 , 1 0 .* ,C C ! 3 > .9^ 5 } ,N , lorn t ,0po , I log , I fn , I j o in ,N jo in 
! ’ S 1 0 C 0R I =N T 0 N 
’ ! 520 Aq=aj l ,!* ' 
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1530 FOR J«M TC N 

!S40 A< I J >«A( I J >/Ac 

1650 NEXT J 

1ES0 C( I )*C< ! > / Ac 

1570 NEXT I 

1590 SUBEND 

1690 SUB Subtract ( K ) 

1700 COM /Cply/ A( 1 0 , 1 0 ) ,C( 1 0 ) ,B( 5 ) t N , Iprnt ,Opo , I log , I f n , I join ,Njoin 

1710 FOR I *K TO N 

1720 FOR J-K-1 TC N 

1730 A( I ,J >*A( K- 1 ,J >-A( I , J ) 

1740 NEXT J 

1750 C( I )*C ( K- 1 >-C( I > 

1750 NEXT I 
1770 SUBEND 
1790 SUB PI in 

1790 COM /Cply/ At 1 0 , 1 0 > ,C' 1 0 > ,B( 5 ) ,N , Iprnt ,Opo , I log , I f n , I jo in # N join 
1900 COM /Xxyy/ X* ' 25 > , Yy ( 25 ? 

1910 PRINTER IS 705 
1920 BEEP 

1930 INPUT "WANT TO PLOT Uo vs Uw? ( 1 «Y ,0-N ) " , Iuo 
1840 IF Iuo-0 THEN 
1950 BEEP 

I960 INPUT "SELECT (0*h/h0X sane tube , 1 -ht HF )/h( sm ) • , Irt 
1970 BEEP 

1980 INPUT "SELECT h/h RATIO ( 1 -FILE ,0-COMPUTED )" , Ihrat 
ie90 IF lhrat*0 THEN 
1900 BSE P 

1910 INPUT "WHICH Tsat ( 1 -6 .7 ,0—2 . 2 )" , I sat 

1920 END IF 

1 930 Xnin-0 

1940 Xr.ax-10 

1950 Xstep-2 

I960 IF Irt *0 THEN 

1 970 Ymin*0 

J 990 Ynax-1 . 4 

1990 Yst-p-.Z 

2000 ELSE 

2010 Ymm*0 

2020 Yna/.»!5 

2030 Ystep*5 

2040 END IF 

2050 ELSE 

2050 Opo*2 

2070 Yrnin^O 

2090 Yrax-12 . 

2090 Ystep-3 
2 ! 00 Xrin**0 
2110 Xrar =4 
2120 X s t ep» 1 
2130 END IP 
2140 IF I hr at * 1 THEN 
2150 Ym : n=0 
2 1 50 Y«*.\ * ! c 
2*70 Ystsp-3 
2 1 90 X r ; r = 0 
2190 Xm3.-=9 
! 2200 X s t - 2 
122 1C END IF 
12220 BEEP 
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PRINT " I N ; 3P ! ; I P ZZ02 3302 3322 5329;” 

133-10 PRINT ‘SC 2 !92 2 !G2:TL 3 2:” 

13350 3f.<-10O/C<naA-Xpun > 

13350 naA-Vrsin) 

13379 PRINT "PU 0,2 PQ" 

13390 PGR X a - X n i n TO Xrra,-; STEP '/ = t e 2 
13392 :<-< Xa-Xnm >*Sf* 

13320 PRINT •pA ,, .v • 0. v T; •• 

!33!0 NEXT Xa 

13230 PRINT "PA 100 ,0 ; P’J ; ” 

13330 PRINT "P'J PA 0,0 P0“ 

13340 FOR Ya-Y*in TO Yma* STEP Yatep 
13350 Y-< Ya-Ymn)*Sfy 
13350 PRINT "PA 0 , - ; Y , " YT ” 

13370 NEXT Ya 

13290 PRINT -PA 0,100 TL 0 3 " 

13390 FOR Xa-Xcun TO Xrtax STEP Xstep 

13409 X-< Xa-Xnin >*S f x 

13419 PRINT "PA" ; X , " ,190; XT" 

13430 NEXT Xa 

13430 PRINT "PA 100,100 P'J PA 100,9 PQ" 

12440 FOR Ya-Yntn TO Ynax STEP Ystsp 

134S0 Y-< Ya-Ynin>*Sfy 

134S0 PRINT "PQ PA 100 " Y "YT" 

13470 NEXT Ya 

13490 PRINT "PA 100,100 PU" 

13490 PRINT "PA 0,-2 SR 1.5,2“ 

12590 FOR Xa-Xnln TO Xnax STEP Xstap 
12519 X»( Xa-Xnin MS Px 
12530 PRINT "PA " ; X , " ,0? " 

13530 IF Iuo-9 THEN PRINT "CP -2 ,-! ?L9" ;Xa; " " 

13540 IF Iuo-I THEN PRINT "CP -1 .5 ,-1 ;LS" ;Xa; * " 

13550 NEXT Xa 

135S0 PRINT "PU PA 0,0" 

13570 FOP Ya-Yfiin TO Ynax STEP Ystep 
1 3590 IF A9S<Ya>< 1.5-5 THEN Ya-0 
13590 Y-( Ya-Ymn )*Sfy 
13500 PRINT "PA 0 , " ; Y , " " 

13510 IF Iao-0 THEN PRINT "CP -4 . 25 ; L3“ ; Ya; " " 

12539 IF Euo-1 THEN PRINT "CP -3 ,- . 25 ;L9" ; Ya ; " " 

13530 NEXT Ya 

1 3940 X labs lS-"Qi i P-rcsnt " 

13S50 IF Iao-0 THEN 

13550 IF Irt-0 THEN 

13S70 YLabel$-"h/H0X" 

13590 ELSE 

13S90 YLabelS-" h/hsmooth" 

13700 EMO IF 

13*10 PRINT "SR t . 5 ,3 ; P'J PA 50,-10 CP " ; -LEN ( X I ab s L 3 } / 3 ; " 0 ; L9 " ; X I at a L S ; " " 

13730 PRINT "PA -11,50 CP 0 , " ; -LEN( Y L ate 13 > /2 *5/5 ; "0 1 9 , 1 ; L9 " ; Y Labe i$ ; " " 

13720 PRINT "CP 0,0" 

13740 ELSE 

13750 PRINT " 3P0; SP3" 

13*50 PRINT "SR 1 . 3 ,3 . 4 ; PU FA -9,35;0I 9,1;LSU;PP l,0.5;L9o;PR -1,0.5;L9 (kW/m 

PRINT "PP -I , 0 . 5 ; SP l,t.S;L33;SR t.S,3;PR .5,.5;L9.:PP .S,0;L9K>" 

PRINT “PA 43,-10:01 ! 0;L2Y sPR .4 -!:L9u;PF 1 . 5:L2I i.' 3 J " 

PRINT " SP9 : SP 1 M 
5 NO IF 
Ipn-0 



1 37*0 
! 7^30 
13732 
’3900 
’3312 
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IS2G SEEP 

2830 INPUT "WANT TG PLOT DATA FROM ft FILE M «\ ,0-N > 7 " 

2940 Icn-0 

2950 IF 01. p* 1 THEN 

22EG SEEP 



2970 

2990! 

2990 

2900 

2910* 

2920 

2930 

2940 

2950 

2960 

Z970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 



INPUT "ENTER THE NAME OF THE DATA FILE " ,C_? i leS 
IF Iuc-0 THEN 
BEEP 

INPUT "SELECT ( 0-LINEAR , ! -LOG( X ,Y >" , 1 log 
END IF 

ASSIGN ®Fi I e TO D.fiief 
BEEP 

INPUT “ENTER THE BEGINNING RUN NUMBER" ,Md 
BEEP 

INPUT “ENTER THE NUMBER OF X-Y PAIRS STORED" ,Npa 

IF Iuo-0 AND Ihrat-0 THEN 

BEEP 

INPUT "ENTER DESIRED HEAT FLUX 



■4X L 
" 4X , “ " 3 
" 4X p * * 5 



END IF 
BEEP 

PRINTER IS 1 
PRINT USING “ 4X 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING "4X/-6 
PRINT USING "4X ,“? 
INPUT Sym 
PRINTER IS 705 
PRINT “PU DI" 

IF Sym-1 THEN PRINT 
IF Sym-2 THEN PRINT 
IF Syr.- 3 THEN PRINT 
Nn-4 

IF I log- 1 THEN Nn- 1 
IF MdM THEN 



.0 



Select a symbol:""" 

1 Star 2 Plus sign" " " 

3 Circle 4 Square""" 
Rombus " " " 

Right-side-up triangle"' 
Up-side-down triangle""* 



“SM*" 
" SM+ " 
" S Me “ 



13190 FOR I-! TO (Md-1 ) 

13190 ENTER eFil-;Xa,Ya 
13200 NEXT I 
13210 END IF 
13220 IF Ihrat-0 THEN 
13230 Ql-Q 

13240 IF I log* 1 THEN Q-LOG(Q) 

13250 END IF 

13260 FOR I-! TO Npairs 

13270 IF Iuc-0 AND Ihrat-0 THEN 

13290 ENTER 0File:Xa ,9' • ) 



3290 Ya*B( 0 ) 



13300 FOR K-! TO Nn 
13310 Ya-Ya A 9(K > - 0 “ K 
13320 NEXT K 
13330 END IF 

13340 IF Iuo-1 OR Ihrat-1 THEN 
13350 ENTER £File;X 3 ,Ya 
13350 IF Iuo=l THEN Ya-Ya.1000 
13370 END IF 

13390 IF Iuc-0 AND I Hr a * -0 THEN 
13390 IF I log-1 THEN Ya-EY Df > a ? 
13400 IF llcg*0 THEN Ya-Ol.'Ya 
13410 IF Irt-0 THEN 
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31 29 if Xa-0 THEM 

3130 /c=y 2 



13110 Ya-I 



i 3450 c L cc 
13450 Y a * Y a / Y o 
13170 END IF 
13490 ELSE 

13450 Ha«-FMHanooth(Q ,Xa ,Isat > 

13500 Ya-Ya/Han 

13510 ENO IF 

13530 ENO IF 

13530 X;<( I-l >-Xa 



13510 Yy< I-l >-Ya 
13550 X*< Xa-Xn 1 n > *5 f .< 
13550 Y-( Ya-Ynin ) *S f y 
13570 IF Sy«;3 THEN PRINT 
13530 IF Sy«<4 THEN PRINT 
13550 PRINT "PA" ,X ,V 
I3S00 IF Sy«>3 THEN PRINT 

I3S10 IF Syn-1 THEN PRINT 

13S20 IF Syn-5 THEN PRINT 

13530 IF Sym-S THEN PRINT 

13540 IF Sym-7 THEN PRINT 

! 3550 NEXT I 

13650 SEEP 



"3M” 

"SR 1 .4 ,2.4" 

"SR l .2,1 .S" 

"UC2 ,4 ,99 ,0 ( -9 ( -4 ,0 ,0 ,9 ,4 ,0, ; " 
"UC3.0.99 .-3, -5 ,-3 ,5 ,3,5 ,3, -5; " 
"UC0 ,5. 3,93, 3, -3, -5 ,0 ,3,3; " 

"'JC0 ,-5.3 ,99 , -3 ,9 ,6 ,0 ,-3,-3;" 



3570 ASSIGN SFila TO • 



3530 5N0 IF 
3599 PRINT "P*J SM" 



3700 9EEP 

3710 INPUT "WANT TO PLOT A POLYNOMIAL (1 -Y ,0-N )? " ,Okp 

3720 IF Okp-I THEM 

3730 9EEP 

3740 PRINTER IS 1 



3750 PRINT USING " 4X , " "Se lect Line type:""" 

3750 PRINT USING "SX,""0 Solid line""" 

3770 PRINT USING "SX,' , "I Qaahed""" 

3790 PRINT USING "SX,""2 ti ,5 Longer line - daah" "" 
3750 INPUT Ipn 
3900 PRINTER IS 705 



3910 9E5P 

3920 INPUT "SELECT ( 0-L INEAR , ! -LQG( X ,Y 5 ) " , I log 
3930 Iprnt"! 

3940 CALL °oly 
3350. IF luo-l THEN 
3950 9EEP 

3970 INPUT "DESIRE TO SET X Lower and Upper Llnit ( 1 -Y ,3-M 
3990 IF Ixlin-0 THEM 
3950 XU-0 
3900 Xu l =*7 



3910 ENO IF 



3929 IF I a 1 1 r * ! THEM 

3930 9EEP 

3540 INPUT "ENTER i Lower L:nit",Xll 
3950 SEEP 

3550 INPUT “ENTER X Upper Linit’, Xul 
3979 END IF 
3999 END IF 

3990 C D° < a =» < 1 1 TO <ul BTEP :<ateo.'2S 

4022 ^ ^ a r ^ n 4 J 

1010 ' a=F K !P« I /' X 3 ‘ 



, lx l in 
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14C2C IF Iuc e ! THEN = 1000 

! 4030 y*' Ya-Ynn .'• = f> 

14040 X-< Xa-Xnm )*S f / 

14050 IF Y:!0 THEN Y»0 

140SO IF Y 100 THEN EG TO !41S0 

1 4070 Pu-0 

14090 IF Ipn-1 THEN Idf-Icn MOD 2 
14090 IF Ipn-2 THEN Idf-Icn MOD 4 
14100 IF Ipn-3 THEN Idf-Icn MOO 9 
14110 IF Ipn-4 THEN Idf-Icn MOD IS 
14120 IF Ipn-5 THEN Idf-Icn MOD 32 
14130 IF Idf-t THEN Fu-1 
14140 IF Pu-0 THEN PRINT "PA* ,X , Y . *PD" 

14 150 IF Pu-1 THEN PRINT "PA" ,X , Y , "PU" 

14160 NEXT Xa 
14170 PRINT "PU" 

14190 GOTO 10920 
14190 END IF 
14200 BEEP 

14210 INPUT "UANT TO QUIT ( 1 -Y ,0-N >7 " , Iqui t 
14220 IF Iquit-1 THEN 14240 
14230 GOTO 12920 
14240 PRINT "PU SP0" 

1 4250 SUBEND 
1 4250 SU8 Stats 
14270 PRINTER IS 701 
14290 J-C 
14290 K-0 
14300 BEEP 

14310 IF Ip let — 1 THEN ASSIGN ©File TO P.fileS 
14320 BEEP 

14330 INPUT "LAST RUM No 7 ( 0-QUIT > - ,Nn 

14340 IF Nn-0 THEN 14700 

14350 Nn-Nn-J 

14350 S*-0 

14370 Sy-0 

14390 Sz-0 

U390 S>.s-0 

14400 Sy s-0 

14410 Scs-0 

14420 FOR 1-1 TO Nn 

14430 J-J+1 

14440 ENTER CF i 1 e ; 0 , T 
14450 H-Q/T 
144S0 S^-S^C 
14470 S/. 5-S>.s +C" 2 
14490 Sy*Sy*T 
14490 Sys-Sys+T^Z 
i 4 500 Sr-Sz-'-H 
U510 Sz2 = Scs j -H^ 2 
14520 NEXT I 
14530 Qave-Sh/Nn 
14540 Tave-Sy /Nn 
14550 Have-Sz/Nr. 

14550 5 de vp-SQR ( AE S ( ( Nr • S .■ z -5 ■ 2 > ( Nrt* t Mn- 1 ) ) ) ) 
1 4 570 Sd-vt-SOR' A95< ( Nr*S . * -S> '‘2 ) ( Nn* * Nr- 1 > ? ) ) 
14 59G Sdevh-SOP( ASS' (Nr*Szs-Sr '2 > "UnM Nr-! > ) ? ? 
14590 SH= 1 00*Sdevh/Have 
1 J 500 S q = 1 00 • S d e p / Q s v e 
I 45 10 S t - ! 00 • S de v t / Ta v - 
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! 15C0 IP K * ! THEM 1 15 59 
j 4930 po i mT ♦ 

1 4540 p p r m f USING "MX, " "0 A TA PILE:'*" !4A" ;File5 
I icc0 PRINT 

14550 POINT USING -fix.- -RUN Htube SdavH Qdp SdevO Thetab SdevT 

US70 K- I 

11590 PRINT USING ' 1 IX ,DQ ,2< 2X ,0. 3DE , ! X ,30 .20 > ,2X ,Q0.3D . 1 X ,30 .20" ? J ,Hava ,Sh ,Qa 

,Sq,Tave,Sl 

14599 GOTO 14323 

14700 ASSIGN ?F i L e ! TO ♦ 

14719 PRINTER IS 1 
1 4729 SURE NO 
14779 SUQ Coe f 

14740 COM /Cply/ A( 10,10) ,C( 10), 8(5) ,N, Iprnt ,Opc , I log , I fn , I join .Njoin 
! 4750 SEEP 

14750 INPUT -GIVE A NAME FOR CROSS-PLOT FILE" .CpfS 
14770 CREATE SCAT Cp f 3 ,2 
14790 ASSIGN GFiis TQ Cpf* 

14790 SEEP 

14900 INPUT "SELECT ( 0-LINEAR , 1 -LOG( X ,Y ) ) " , I tog 
14310 9EEP 

14320 INPUT "ENTER OIL PERCENT ' - 1 -STOP ) * ,9cp 

14830 IF 9op<0 THEN 14370 

14340 CALL Poly 

U9S0 OUTPUT 0Fiie:9op ,9( • > 

14960 GOTO 14810 
14970 ASSIGN 0File TO * 

14990 SU8EN0 

U990 SU8 UiisanICf ,C: > 

14900 COM /Uil/ 02 ,0i ,0o ,L ,Lu ,Kcu 
14910 DIM Enf ' 1 2 ) 

149201 UL ISON PLOT SUBROUTINE DETERMINE CF AND Cl 
14930 9EEP 

11940 INPUT "ENTER DATA FILE NAME", FileS 

14980 EEEP 

14950 PRINTER IS 1 

14970 PRINT USING "4X .""Select option:""" 

14590 PRINT USING "4X,"" 0 Vary Cf and Ct""" 

14990 PRINT USING "4X,"" I Fix Cf Vary Cl""" 

15000 PRINT USING "4X,"" 2 Vary Cf Fix Cl""" 

15010 INPUT "ENTER OPTION" , Ic f lx 

15020 PRINTER IS 701 

15030 IF Ic fix-0 THEN 15070 

15040 IF lcfix'0 THEM SEEP 

1E0E0 IF Icfix-1 THEN INPUT "ENTER Cf.Ca? 

1E050 IF Icfix-2 THEN INPUT "ENTER Cl" , C i 

1E070 PRINTER 15 1 

15090 INPUT "Uant To Vary Coe f f 1 -Y ,0-N ) " , Iscce f 
1E090 IF Icco-f-1 THEN INPUT "ENTER COEFF" ,R 
IE100 PRINTER IS 701 

15110 IF Icfix-0 OR Icfi/.-2 THEN Cfa-.304 
15120 IF I c f i ’ THEM Cfa-Csf 
i ci -0 C:a=»C: 

! E 1 40 <n=*.3 
1 E 1 EG c r-.3 
1 E ! 50 Jj-0 

I c | -rj o r ^2 . 

15190 IF I c c oe f - I THEN °r’° 

’ E 1 90 PR INTER I S ! 
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) 5200 PRINT Qc ,D i ,kcu 
15210 ASSIGN @Fi I * TO File* 

’ 5220 ENTER GFileiNrun .Data* .Ldtcl .LdteZ ,1 tt 
’5230 Rw-Do*LCS( Dc/O: >/( 2*Keu ) 

15240 Sx-0 

15250 Sy-0 

15250 Sxy-0 

15270 S>.2-0 

15280 Sy2-0 

15290 FOR 1-1 TO Nrun 

15300 ENTER 0Fi le ; Bop ,To 1 d£ ,Emf ' • ) ,Fms 

15310! CONVERT EMF’S TO TEMPERATURE 

15320 FOR J-l TO 5 

15330 T( J >-FNTvsv(Enf( J) ) 

15340 NEXT J 

15350 Tsat-( T( 1 ) + T(2> >*.5 
15350 Tavg-T(S) 

15370 Grad-37 . 9S53 4 - . 1 04399* Tavg 
1 5390 T drop -Em f ( 7 )• 1 . E^S / ( 1 0 . # Gr ad > 

15390 Tavgc*T( 5 )-Tdrop • . 5 

15400 IF ABSCTavg-Tavgc >>.01 THEN 

15410 Tavg-( Tavg+Tavgc )+ .5 

15420 GOTO 15370 

15430 END IF 

15440! 

15450! Compute properties of water 
15460 Kw-FNKw(Tavg) 

15470 Muwa*FNMuw( Tavg ) 

15490 Cpw-FNCpu(Tavg) 

15490 Prw-FNPrw( Tavg ) 

15500 Rhow-FNRhow( Tavg ) 

155101 

155201 Compute properties of Freon-114 

15530 Lr t d-T drop / LOG ( ( T ( 5 >-Tsat )/ ( T( 5 )-Tdr cp-Tsat ) ) 

15540 IF Jj-0 THEN 
15550 Tw-Tsat +Fr *Lmt d 
15550 Thetab-Tu-Tsat 
1 5570 J * - 1 
15590 END IF 
15590 TF-( T w+T sat )• .5 
15500 Rho-FNRho( Tf > 

15510 Mu-FNMu(Tf ) 

15620 K-FNK ( T f ) 

15530 Cp-FNCp ( T ? ) 

15540 Be t a-FNBe t a ( T f ) 

15550 Hf g-FNHf g( Tsat ) 

15550 Ni-Mu/Rho 
15570 Alpha-K/( Rho*Cp > 

15590 Pr-Ni /Alpha 
15590’ 

1 5700 1 Ana lysis 
15710’ COMPUTE MOOT 
15720 A-PI*<Qe"2-0i A 2 >/4 
t 5^30 P«P I • Oc 

!5“40 Mdot -3 . 9557E-3+Fms • '3.5l955E-3-Fms*( 9.2200SE-5-F «-•' 1 . 23S9SE-7-Fms *4 .2'99“ 
5-10))) 

15750 Q-Mdo t *Cpu* Tdrop 

I 5“G0 1 COMPUTE NATURAL-CONVECTIVE HEAT - T D ANSFE? COEFFICIENT 
1 E - "*?© ' FOP UNENHANCED END'S) 

15790 Hbar-190 
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15790 Fe-< Hbar* D /( Kcu*A ' ? .5*Lu 
15900 Tanh-FNTanh( Fe > 

1591C Theta»Thetab*Tanh. /C e 

'.9920 Xx*( S.9! •9eia*Theiab*De*3*Tanh/< Fe*Ni * Alpha > >*’ . !££££7 
15930 Yy-( 1 +( . 559/Pr >" < 9/ 1 6 ) r ( 9/2? ? 

15940 Hbarc-K/Oo*'. .£* .3S7*X*/Yy > A 2 

15950 IF A9S( (Hbar-Hbarc )/Hbar>; .001 THEN 

15950 Hbar-(Hbar+Hbarc >• .5 

15970 SOTO 15790 

15990 END IF 

159901 

159001 COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENOS 

15910 Q1 -<Hhar-P*Kcu*A r . 5*Thetab *Tanh 

15920 0c»0“2*0l 

15930 As-P ! *D2*L 

159401 COMPUTE ACTUAL HEAT FLUX 

15950 Qdp-Qc/As 

1 5960 IF Icfix-0 OR IcfixM THEM Cs f- 1 /Cf <1 . /Rr ) 

15970 Thet ab-Cs f /Cp*Hf g* ( Qdp/( Mu*Hf g )♦ ( .01 4/« 9.01 •Rho ? >* .5 >“< 1/Rr >*Pr-l 

15990 Ho-Qdp/Thetab 

15990 Omaga-Ho/Cf 

1G000 Ue-0/<PI*Do*L*Lwtd ) 

16010 Vw“Mdot/< Rhow*PI*0i *2/4 ) 

16020 Rsu-Rhou*Vu*Oi /Muua 
16030 Tui-Tw+Q*Ry/( PI *0o«L ) 

1 6040 6ama-Ku/Di # Reu A . 8*Pru*< 1 /3 . )•( Muua/FNMuu( Twl ) >* . 1 4 
160501 PRINTER IS 1 
16050 Yw-( 1 . /Uo-Rw )*0nega 
16070 Xu-0nega*0o/< 5ana*Di ) 

16090 Sx-S>:+Xu 
16090 Sy-Sy+Yu 
16100 S»:y“Sxy+ Yu # Xu 
16110 Sx2*Sx2+Xw*Xw 
1 £ 1 20 Sy2*Sy 2+Yu* Yu 
15130 NEXT I 
1 S 1 40 ASSIGN OFile TO * 

l E 1 50 Sx * Sy-Nrun* Sx y ) / ( Sx * Sx-Nrun* S *2 > 

16160 C-<Sy-Sx*M )/Nrun 

16170 IF Icfix-0 OR Icfxx-3 OP Icfix-4 THEN 

16190 Cic-1/M 

16190 Cfc-1/C 

16200 END IF 

16210 IF Icfix-1 THEN 

16220 Cic-1/M 

16230 Cfc-Cf 

1S240 END IF 

16250 IF Icfix-2 THEN 

16260 Cic-Ci 

16270 C f c » ! / C 

1 £290 END IF 

1S290 IF ABS( (Ci-Cic )/Cic ): .001 OR A8S( CCf-Cfc >/Cfc K .001 THEN 
1 £300 Ci-< Ci+C: c )• . 5 
16310 Cf »< Cf+Cfc ) * . 5 
15320 PRINTER IS l 

16330 PRINT USING M 2y."" Cs f - H M ,MZ . 3DE ,2X , " " Ci * " " ,MZ . 3DE " ; Cs f ,C i 
i £ 34 0 PRINTER IS 701 
16350 SOTO 15210 

1 C 3C0 CMH 7 C 

16370 PRINT 

15390 ° Pi N T E P IS 70! 
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1E390 PRINT USING "ZZV. “ " CP Ci"** 

15400 PRINT USING “3X ""ASSUMED " " MZ.30E 3X MZ.3CE";Cfa Cia 

15410 PRINT USING " 9X /’“CALCULATED " " ,M2 . 3DE , 3X ,MZ . 3DE " ; Cs f ,C t 

1S430 PRINT 

15 430 SunZ -Sy2-Z *1 • Sxy -3 *C‘Sy 2 •S.<:*Z*M*C*S^ i NrurrC - : 

15440 PRINT USING " 1 0X , “ " Sun of Squares » " “ , Z . 3QE " : Sur 3 
15450 PRINT USING " ! 0X . " "Coef f ic Lent * " ” ,0 . ODD” : Rr 
IS 450 SUSENO 
15470 DEF FMMuu< T ) 

15493 A«247.9/< T+133. IS ) 

15490 Mu-2.4E-5*I0*A 
ISE00 RETURN Mu 
15513 FNEMQ 
15533 OFF FNCcw(T) 

15533 Cqu-4.21 1209S9-T*' 2 . 2S32SE-3-T ♦ < 4 . 4Z3S 1 E-5+2 . 7 1 4Z9E-7*T ) > 

15540 RETURN Cpu*!303 

15550 FNEMO 

16550 0EF FNRhou(T) 

16570 Ro-999. 52946^1 ♦< .3126 9- T*<5. 43251 3E-3-T ♦ I .234147E-5) > 

1 5590 RETURN Ro 
16590 FNENO 
16600 DEF FNFrui< T ) 

15610 Prw-FNCpu< T >*FNMuu( T )/FNKw< T ) 

15620 RETURN Pry 
16630 FNENO 
I 6640 QEF FNKu< T ) 

16650 X-< T*273.1S>/273. 15 

16660 Ky— .92247+X*< 2.939S-X*< 1 .9007-XM .52577-. 07344 *x ) ) > 

16670 RETURN Ku 
16690 FNENO 
15690 SUB Plot 

16700 COM /Cply/ A< 1 0 , 1 0 > ,CM0 > ,9<S ) ,Noq , Iprnt ,Qqo , I log , I fn ,1 join ,M jo in 
1 ST! 0 DIM 9s<3> 

15723 INTEGER Ii 
16730 PRINTER IS 1 
16740 ldv»0 
15750 SEEP 

16760 INPUT "LIKE DEFAULT UALUES FOR PLOT ( 1 -Y ,0-N >7 ‘ , I dv 
16779 Opc*0 
15793 SEEP 

16790 PRINT USING "4X ," "Select Option:""" 

16900 PRINT USING "SX,""0 q versus delta-T""" 

16910 PRINT USING "6X,""1 h versus delta-T""" 

16923 PRItJT USING "6X,""2 h versus q""" 

16930 INPUT Opo 
15940 9EEP 

16953 INPUT "SELECT UNITS < 9*S I , 1 -ENGLISH > " , Iun 
I6SS0 PRINTER IS 705 
15973 IF I dv \ > 1 THEN 
15990 9 EEP 

19393 INPUT "ENTER NUMBER OF CYCLES FOR X-AXIS", C,; 

! 5903 SEEP 

15310 INPUT "ENTER NUMBER QF CYCLES FOR Y-AXIS" ,Cy 
15923 BEEP 

16930 INPUT "ENTER MIN X-UALUE (MULTIPLE OF !0>",.<n:n 
16940 SEEP 

IS9E0 INPUT "ENTER MIN Y-UALUE (MULTIPLE OF !0>",Ymm 
15950 ELSE 
1 5979 IF Opc-0 THEN 
15390 Cy-3 
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THEN 



IS 390 C ’ 3 
I 70Q0 X.n : n* . ! 

I 70 in 7 n * n «| 0 Q 
j 7070 cmq r c 
j 7.370 rc n p c ^ i 
! ^040 Cy-2 
! 70S0 C*»2 
! 7QS0 Xn:n- . ! 

17070 Ynin-100 
17030 END IF 
17030 IF Opo-2 THEM 
17100 IF Iun-0 THEM 

171 19 Cy-3 

17120 04-2 
17130 Xn in- 1000 
! 7 1 40 Yn i n- l 00 
17150 ELSE 
171S0 Cy-3 
17170 04-3 
17190 Xnin-100 
17190 Ymn* I 0 
17200 END IF 
17210 END IF 
17220 END IF 
17230 9EEP 

17240 PRINT " IN ?SP 1 ; IP 2300 ,2200 , 9300 ,5800 ? " 

17250 PRINT "SC 0,100,0 J00;TL 2,0:* 

172S0 S fx» 1 00/C* 

17270 Sfy-100/Cy 
17290 BEEP 

17290 INPUT -WANT TO BY-PASS CAGE * 5 ( I -Y , 0 -N ) " , Ibyp 
17300 IF Ibyp- l THEN 19540 
17310 PRINT "PU 0 .0 PO" 

17320 Nn-9 

17330 FOR 1-1 TO C4+1 

17340 Xat-Xfnin*l0' k ( I-t 1 

17350 IF I-Cx+l THEM Nn-l 

17350 FOR i-l TO Mn 

17370 IF J-l THEN PRINT "TL 2 0" 

17390 IF J-2 THEN PRINT "TL I 0" 

17330 Xa-Xat * J 

17400 X-LSTI Xa/Xnin )*S fx 

17419 PRINT "PA" ; X , " ,0; XT;" 

17420 NEXT J 
17430 NEXT I 

17440 PRINT "PA 100,0; PL 1 ;" 

1 7450 PRINT "PU PA 0.0 PO" 

! 7450 Nr. -9 

17470 FOR I-l TO CyH 
17490 Yat-Ynin* 1 0' ( I-! > 

17490 IF I -Cy > 1 THEM Nn-l 

17500 FOR J-l TO N n 

17510 IF J-l THEM PRINT " TL 2 0" 

17520 IF J-2 THEN PRINT "TL 1 0" 

17530 Ya-Ya t • J 

17540 Y -LG T < Y 3 / Yn 1 n ) • S f y 

I 7550 PRINT "PA 0 , * ; Y * Y T “ 

«-cc 0 NEXT J 
1 ~'C”n ricy r r 

1 “*590 DD INT “°A 0 ! 00 r L 0 2 " 
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i 7Egr> MnsQ 

17600 FC C ! * ! TO Cx A l 

1 7E 1 0 Xat-Xoun* 10" I- 1 1 

17520 IF I »C^ A ! THEN Nn-1 

1 7630 FOP J-! TO Nn 

17540 IF J-l THEN PRINT M TL 0 2 " 

1 7650 IF S'. 1 THEN PRINT "TL 0 I " 

1 755G Xe-XaW 
1 7S7G X-LGT ( Xa/XniR )*5 f.*. 

17530 PRINT "PA" sX ,100? XT’* 

176SG NEXT J 

17700 NEXT I 

17710 PRINT "PA 100,100 PU PA 100,0 PD" 

17720 Nn-9 

17730 FOR 1-1 TO Cy+1 

17740 Yat-YnirrlG"' I-! J 

1 7750 IF I-Cy+1 THEN Nn-1 

17750 FOR J-1 TO Nn 

17770 IF J-1 THEN PRINT "TL 0 2" 

17790 IF JM THEN PRINT "TL 0 1" 

17790 Ya-YaW 

17900 Y-LGT ( Ya/Ymi n > ♦ S f y 

17810 PRINT "PD PA 100,",Y,"YT" 

17820 NEXT J 

17930 NEXT I 

17940 PRINT "PA 100,100 PU" 

17950 PRINT "PA 0.-2 SR 1.5,2“ 

17950 Ii-LGT (Xnin) 

17970 FOR 1-1 TO Cx+1 
17990 Xa-Xni n* 1 0* ( 1-1 ) 

17990 X-L6T (Xa/Xnin )*Sfx 
1 7900 PRINT "PA" ;X ," .0; " 

17910 IF 1 1 •— 0 THEN PRINT "CP -2 ,-2 : L0 1 0 ; PP -2 ,2 ; L8 " ; I i ; " " 
17920 IF !i>0 THEN PRINT "CP -2 ,-2 ; L9 1 0 ; D P 0 ,2 :LB" ; I : : " " 

17930 Ii-Ii+1 

17940 NEXT I 

17950 PRINT "PU PA 0 ,0" 

17950 I l-LGT (Ynm) 

17970 Y10-10 

17990 FOR 1-1 TO Cy+1 

17990 Ya-Ymin* 10" < 1-1 ) 

18000 Y-LGT < Ya/Ynin >*S f y 
19010 PRINT "PA 0 , " ; Y , * " 

19020 PRINT "CP -4 ,-.2S:L9!0;PR -2 .2 ; L8" ; 1 1 ; " " 

19030 Ii-Ii+I 
1 9040 NEX T I 
19050 SEEP 

19050 INPUT “WANT USE DEFAULT LABELS ( 1 -Y ,0*N J? " , I dl 
19070 IF Idi'. >1 THEN 
19090 BEEP 

19090 INPUT "ENTER X-LA9EL" ,X labels 
19100 9EEP 

19110 INPUT "ENTER V-LABEL” YlabelS 

19120 END IF 

12130 IF Ocs\2 THEN 

19140 D PINT "SP 1 2;PU C A 4 0 - 1 4 • ’* 

19150 PRINT "L9' T :PP -1.5,3 C D CQ 1.2,0 PUsPP . 5 , -4 ; LSuic ; PP .5 
’ E 1 EC C PINT ”LB-T;°P .5 -! : oc .5 1 : ” 

121^0 IF lan=0 THEN 
’9190 P P I N T " L 9 1 
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i c i go ELSE 

JS20C PRINT "LS) IP)” 



9210 END IF 
£220 END IF 



19230 IF 0po-2 THEN 
15240 IF Iun-0 THEN 

15250 PRINT “SR 1.5,2;FU PA 40,-U;LBq (U/m;SR 1,1.5;PR 0.5,1;LB2;5R 1.5,2:PF 
0, E 1 ;L9 >" 

19250 ELSE 

19270 PRINT “SR t.5,2:PU PA 34,-14;LBq (Btu/hr;PR .5,.5;LB.:PR .5, -.5;' 

19290 PRINT "LBf t ; PR .5,1;SR 1,1.5?LB2;SR 1.5,2;PR .5,-1 ;L9 ) : " 

19290 END IF 
19300 END IF 
15310 IF Qpo-0 THEN 
19320 IF Iun-0 THEN 



5330 PRINT "SR 1 . 5 ,2 ;PU PA - 1 Z , 40:01 0,1;LBq (U/n;PR -1,0.S;SR 1,1.5; LE2 : SR l 
5 2 ; PR 1 . 5 ; LS ) " 

5340 ELSE 

0350 PRINT "SR 1 .5 .2 ;PU PA -12,32501 0,l;LBq (Btu/hr;PR -.5 , .5;LB. ;PR .5,.5; M 
9350 PRINT "LBf t « SR 1,I.5;PR -t,.5;LB2;PR 1,.5;SR 1 .5 ,2 sLB >“ 

9370 ENO IF 
8380 END IF 
9390 IF Opc>0 THEN 
9400 IF Iun-0 THEN 

94101 PRINT “SR 1 .5 ,2 ;PU PA -12,39:01 0,!;LBh (U/m<PR -1,.5;SR l,1.5;L92;SR 1 
5 ,2 * PR .5 , .5: * 

0420 PRINT "SR 1 .2 ,2.4; PU PA -12 ,37*01 0,1<LBh:PR l,0.5;L8o;PR -1,0.S;LB (U/m 



9430 PRINT "PR -1,.S;SR 1,t.5:LB2;SR !.5,2:PR .5,.5;LB,:PR .5,0;LBK>“ 

9440 ELSE 

9450 PRINT "SR l.5,2;PU PA -12,29:0! 0,1;Leh (9tu/hr;PR -. 5 , . 5 ? LB . ; PR .5, .5;" 

9450 PRINT "LB ft ;PR -1,.5;SP 1,?.5;LB2;SR 1.E.2;PR .5,.5;L5.;PR .5,.5;LBF>" 

9470 END IF 

9490 END IF 

e490 IF IdI-0 THEN 

9S00 PRINT “SR 1 .5 , 2 ; Pl l PA 50.-15 CP -LENC X 1 abe 1$ >/2 : “0 ; L9" : X label S “ 

9510 PRINT "PA -14.50 CP 0 . “ : -LEN( VlabeiS J/2-5/S : "01 0 . 1 ;LB“ ; Ylabe! S ; " ’* 

9520 PRINT "CP 0,0 01" 

9530 END IF 
9540 Ipn-0 
5550 XII-1 .5+5 
8550 Xui — 1 .E+6 
9570 Icn-0 
9590 Ifn-0 
5590 Ijs:n-I 
9500 BEEP 

9G10 INPUT “UANT TO PLOT OATA FROM A FILE ( 1 -V ,0*N J7 " ,Q1 
9520 IF 01-1 THEM 
9530 BEEP 

9540 INPUT “ENTER THE NAME 0 C THE OATA FILE” ,D_f 1 1 e$ 

9550 ASSIGN £Fii- TO D f i!eS 
9550 BEEP 
95“0 BEEF 

9590 INPUT " ENTER THE' BEGINNING RUN NUMBER" ,Md 
9E90 BEEP 

5700 INPUT " ENTER THE NUMBER O c >'-V PAIRS STORED" .Npairs 

OT10I DCCD 

2720 1 INPUT "CONNECT DATA L'I T H LINE < ! »v C*N1 7 " Ici 

B~30 9EE r 

9740 PP INTER IS 1 



1 14 



13750 


PRINT 


USING 


MX 


"Ss 


lact a symbol: 


i Q“C0 


OO TMT 


USING 


" cv 


’* 1 


Star 2 Pius sign 


19770 


CO TMT 


IJC T * J 0 


• C V •• 


" 3 


Circla 1 Square 


i arar* 


CO TMT 


US I NG 


" C V •• 


•• c 


Rcnbus 


19790 


PRINT 


USING 


•*sx / 


" 5 


Right -a ide-up triangle" 


i 9909 


P P. I N T 


USING 


•• c y " 


" 7 


Up - s t da -down triangle 


19910 


INPUT 


Syr 








i 5920 


PRINTER 15 7 


05 






1 9930 


PRINT 


"PU 01 


" 







13940 IF 9yn*1 THEN PRINT " EM*" 
! 9930 IF Syn-2 THEN PRINT 
1 99S0 IF Syr-3 THEN PRINT •‘SNo" 
1 9970 IF Md, ! THEN 
1 9990 FOR I ■ 1 TO f ttd- ! > 

1 9990 ENTER 0File; /a ,Xa 

19900 NEXT I 

19910 END IF 

19920 FOR I-J TO Npazrs 

19930 ENTER 0Fii-;Ya ,Xa 

19940 IF I - ! THEN Ql-Ya 

189S0 IF I-Npatra THEN Q2-Ya 

1 9960 IF Opo-1 THEN Ya-Ya/Xa 

19970 IF Opo-2 THEN 

19990 Q-Ya 

19990 Ya-Ya/Xa 

19000 Xa-Q 

19010 END IF 

19020 IF Xa<X ! 1 THEN Xii-Xa 
19030 IF Xa>Xul THEN Xul-Xa 
19040 IF Iun-1 THEN 

190E0 IF Cpo<2 THEN Xa-Xa*1.8 

19050 IF Opo:0 THEN Ya*Ya*.!75l 

19070 IF Opo-0 THEN Ya-Ya*.3J7 

19090 IF Opc-2 THEN Xa-Xa*.3l7 

19090 END IF 



9100 X-LGT ( Xa/Xntn >*S ?.< 



19110 
19120 
13130 
19140 
19150 
19150 
19170 
19190 
19190 
! g^ 0 Q 
19210 
1 9220 
1 9230 
! 9240 
1 9250 
19250 
! 92"0 
I 9230 
1 9290 



Y-LGT< Ya/Yrun>*Sfy 
K j-0 

CALL Synb(X ,Y .Syn.Iel ,KJ ) 
GOTO 19270 



IF Syn>3 THEN PRINT 
IF Syn<4 THEN PRINT 
IF IcL-0 THEN 
PRINT "PA * ,X , Y ” 
ELSE 

po r i\| T * pA •’ x y •• on ■ 
END IF 

IF Syn 3 THEN PRINT 
IF 3yn-4 THEN PRINT 
IF Syr-5 THEM PRINT 
IF Syn-5 THEN PRINT 
IF 3 y n - 7 THEN PRINT 
NEXT I 
PRINT "PLT 



-sir 

"SR 1.4 ,2.4- 



'* SR I .2 ,1 .5" 

"UC2 ,4 , 9 3 , 0 , - 9 , - 4 ,0,0,9 ,4 , 0 : " 
"UC3 ,0,99 ,-3.-5 .-3 ,5 .3 ,S ,3.-5: - 
"UC0 ,5. 3, 99, 3,-9, -5, 0,3,9:" 

"UC3 ,-5.3 .99 ,-3 ,9 ,5 .0 ,-3 ,-9: * 



9EEP 



19300 INPUT "WANT TO LABEL 7 ' I -Y ,0-N ) - , I lab 
19310 IF Iiab-1 THEN 
! 9320 F°INT "SP0;SR2" 

'9330 9EEP 

19340 I c Klab-O r HEN 
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! 3350 < L 3 fc - S 
» 93SQ V lab-95 

13370 INPUT "ENTER INITIAL X ,Y LOCATIONS" , XLab ,Ylab 
19390 Xit-Xiab-5 
19590 ^tt-ZLab^S 
! 3 -ICO PRINT " 9R 1,1.5" 

19413 PRINT "SMsPA" .Xtt ,Ytt ,“L9 - Heat Fils" 

19430 Yt t-Y 1 1-3 

19433 PRINT "PA" ,Xt t ,Yt t , "L9 Oil Flu* Mars" 

19440 IF Syr?- 1 THEM PRINT "EM*" 

19450 IF Sym-2 THEM PRINT "SM*" 

19453 IF Syn-3 THEN PRINT "SMo" 

19470 K lab* I 
19490 ENO IF 
1 9493 K j — 1 

19500 CALL Symb(Xlab .Ylab .Syn.IcL ,Kj > 

19513 PRINT "SR 1,1 .5:SM" 

19523 IF Syn'4 THEM PRINT "PR 2,3" 

19533 SEEP 

19540 INPUT "ENTER BOP" .Sop 

19553 IF Sop< 1 3 THEM PRINT "PR 2 ,0 ; L9 " ? Bop ? " " 

19553 IF Bop >9 THEN PRINT "PR .5 ,3;L9" ;8op f " " 

19573 Ihf-3 

19533 IF 01-02 THEN Ihf-I 

19593 IF Ihf »0 THEN PRINT "PR 4,3;L9Inc" 

19533 IF Ih f- 1 THEN PRINT "PR 4,3;L90ec" 

19513 PRINT "PR 2 .SjLS" ; 0_?ils$; ** 

19523 PRINT "SP3;SP!;SR 1.5,2" 

19533 Y lab-Ylab-5 
1 9543 ENO IF 
19G50 BEEP 

1 9563 ASSIGN OFile TO * 

19570 Xll-Xil/1 .2 
19593 XuL-Xul ♦ 1 .2 
1 9G93 1 SOTO 9040 
19733 ENO IF 
19710 PRINT "PU SM" 

I 9723 9EEP 

19730 INPUT "WANT TO PLOT A POLYNOMIAL (l -Y ,3-N >? " ,So_on 
19740 IF So_on«1 THEN 
j g7cg occp 

197S0 PRINTER IS 1 

19773 PRINT USING "4X, ""Select Line type: 

19793 PRINT USING "SX.""3 Solid lire""" 

19793 PRINT USING "5X,""1 Oaah-d""" 

19933 °RINT USING "5X , " "2 . . .5 Longer line - daah 

’9913 INPUT Ipn 
19920 PRINTER IS 705 
1 9923 SEEP 

19940 INPUT "SELECT < 0-L IN , I -LOS« X ,Y > >" , I log 

1 9953 Iprnt-1 

199S3 CALL Poly 

19970 IF Ifn-0 THEM 

19990 9EEP 

19990 INPUT " ENTER NUMBER OF C !LES TO JOIN? " ,Njc in 

I OQQ0 c^n TC 

19910 I • c i 0 

19920 IF I F r Njo m ’HEN Ij-in-l 
1 9933 IF IFr-3 OP tjcm-i r HEM 
! 99 10 FOP ! ,a 0 : 0 I 
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i aac/* a al r , -.a-, r . r . ^ 

•-*-'***■ w - ' ‘ J — - = * J w.Lj. 

! gacQ mcvt r. 

I ag7Q t p«» r f «i ( 

t gga0 c\in rc 



; 9990 rc Nj^in-Pn THEN 
O90CG FOR !j*G TO 3 



00010 Q( t j ) s q 5 ( i j ?/Mjo tn 

00000 9s( : j >-o 

00030 NEXT Ij 
00040 I fn— 0 
0GGS0 ELSE 
200E0 GOTO ! 9SG0 
00070 END IF 
00090 9EE? 

00090 INPUT "ENTER V LOWER ANO UPPER LIMITS" , / L 1 ,/ul 

00100 FOR Xx-0 TO Cx STEP C/-./O00 

00! ! 0 Xa-Xn in* 1 0"Xx 

20100 IF Xa< XI 1 OR X a > X u l THEN 00390 

20130 Icn-Icn+! 



20140 

20150 

00150 

20170 

20190 

20190 

00200 

29210 

20220 

20230 

00240 

20250 

00050 

20270 

00290 

00290 

20300 

20310 

20320 

20330 

20340 



Pu-0 

IF Ipn-I THEN Idf-Icn MOO 2 
IF Ipn-2 THEN Idf-Icn MOO 4 
IF Ipn-3 THEN Idf-Icn MOO 9 

IF Ipn-4 THEN Idf-Icn MOO IS 

IF Ipn-5 THEN Idf-Icn MOO 29 

IF Id?-! THEN Pu-1 

IF Op o-0 THEN Y 3 -FNP 0 1 y ( X 3 ) 

IF Opc-2 ANO I log-0 THEN Y 3 -X 3 /FNP 0 ly ( Xa ? 
IF Opc-2 ANO I log- 1 THEN Ya-FNPo ly < X 3 ) 

IF Opo-1 THEN Ya-FNPcly(Xa> 

IF Ya< Yn tn THEM 20390 

IF Y 3 \ Y ! 1 OR Y 3 > Y u i THEN 20390 

IF Iun- 1 THEM 

IF OpoO THEN X a » X a ♦ ! .9 

IF Opo. 0 THEN /a-Ya*. 175! 

IF Gpo-0 THEN Ya-Ya*.317 
IF Opo-2 THEN Xa-Xa* .317 
ENO IF 

Y-LST ( Ya/Ymm >*S f y 
X-LGT ( Xa/Xmi n > *S f x 



20350 IF Y<0 THEN Y-0 

20350 IF YM 00 THEN GOTO 20390 

20370 IF Pu-0 THEN PRINT "PA" ,X , Y , "PO " 

20390 IF Pu-I THEN PRINT "PA" ,X , Y , "P»J* 

20390 NEXT Xx 
20400 PRINT ,,D U“ 

20413 GOTO 19600 
20420 ENO IF 
20430 9EER 

09440 INPUT "WANT TO PLOT REILL 1 ' S DATA? ( I -Y ,0-N ?" , Ir ly 

09459 IF Qpc-0 OR 0pc-l THEN 

09460 XI 1-3 



09470 Xu 1-00 
00430 END IF 



00490 IF Opc = C TU EU 



09500 XL 1-19000 
09519 Xu 1 - 1 G00G0 
09509 ENO IF 



09509 IF Irlv=! TU EU 
09540 MI-09 
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2OE50 Yu l -70 
2GE50 SEEP 

20E70 INPUT "ENTER LQUER AMO UPPER Y-LIMITS c OP PLOTT IMG " . Y 1 1 , Yu I 
23E9G FOR x*»0 TO C.*. STEP C.-. •000 
2 05 SO Xa-Yn i n * ! 0 ' X • 

20500 IF v a .XI I OR Xa. Xul THEN 20730 
20510 Xl-LOG(Xa) 

20520 IF Cpo-0 THEM V l — 2 . S402937E- 1 + < I *U . 3720 t E I -X l *2 . S 1 347S7E- I > 
20520 IF Qpc-1 THEN Y l --2 . 5402937E- 1 +X i *( 3 . 9720 1 5 ! -X I *2 .5 ! 24797E-I > 
206*10 IF Opc-2 THEN Yl — 3 . 707390 1 E-1 +X l *' 9.72591 90E-1 1*6. 9925942E-3 ) 
20650 Ya-EXPI Y1 ) 

20550 Y-LGTC Ya/Ymm )*SFy 
20570 X-LST' Xa/Xrun >*Sf>, 

20590 Ipu-0 

20690 IF Y-;VI1 THEN Ipu»t 

20730 IF Y> Yu I THEN GOTO 23730 

20710 IF Ipu-0 THEN PRINT "PA" ,X , Y , "PO " 

20720 IF Ipu- 1 THEN PRINT "PA" t X , Y , "PU" 

20730 NEXT Xx 
20740 PRINT "PU* 

20750 ENO IF 
20760 SEEP 

20770 INPUT "WANT TO PLOT ROHSENOU CORRELATION? < I -Y ,0-N ) " , Irohs 

20790 IF Irohs- I THEN 

20790 YU-15 

20900 Yul -93 

20910 9EEP 

23920 INPUT "ENTER Taat ( Dag Cr.Taat 
20920 Ca f- . 0040 
23940 SEEP 

20950 INPUT "ENTER Car < OE c -0 . 004 > " ,Ca f 
20950 T f-Taat ■*•2 

20970 FOR Xx-0 TO Cx STEP C*/200 
20990 Xa-X«in* 1 0" Xx 

20999 IF Xa< X i l OR Xa,V<ul THEM 21179 
20900 Xl-LOG(Xa) 

20910 IF Opo< 2 THEN Tf-TaaUXa/2 
20920 Rho-FNRho( T f ) 

20920 K-FNK ( T ? ) 

20940 Mu-FNNu(Tf) 

20950 Cp-FNCp ( T f ) 

20950 Hf g-FNHf g( Taat > 

20973 N l -Mu/Rho 
20990 Pr-Cp ♦Mu/K 

29993 Onaga-Caf *Hfg/Cs*< < .014/(9.91 *Rhc > > * . 5/ < Mu *H fg > r ( 1 . /3 J *Pr " t . 7 

21000 IF Opc-0 THEN Ya-< X a /Onega ) *2 
21010 IF Opo-l THEN Ya-' Xa/Gr?eg3 > '3/Xa 
21023 IF Qpo-2 THEN Ya-Xa" < 2 . /3 >/Orega 
21323 IF Opc-2 THEN 

2 1 040 Tfc-Tsat *Xa/Ya* . E 
21050 IF ASS ( T f-T f c ). .01 THEN 
2 1 350 T*-< T f X T Fc >♦ .E 
21070 GOTO 23920 
21093 END IF * 

21093 ENO IF 

21100 Y-LST(Ya/Ynin?*Sfy 
2 1 1 13 /-»LGT f /a /'run 1 * S ^ 

21 ’39 I p u = 0 

211-0 IF f • f \i THEM I n u -= j 

2 ! > 10 [c / v„i rucrvj I ! ’ ~0 
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ten IF 1^*0 THEM F°INT ” p*" y / "PQ” 

2\ I C 0 rc U,J 3 I JUCM DOTMT "OA" ;< ✓ "OIJ' 1 

2!! 73 ME XT X.< 

2 ! ! 93 PRINT " PU “ 

2 ! ! 90 END IF 
2 ! 200 9EE° 

2?2?3 INPUT "UANT TO OUIT 4 . 1=V 3=N ?’,[qt 

2)220 IF !qt a ! THEN 212-10 

21223 GOTO 135=3 

21240 PRINT "PU PR 3 3 3PG” 

21250 5U9ENO 

21253 SU9 Svwbl X ,Y f 5yji f Ic! # K j ) 

21270 IF Syn. 3 THEN PRINT "SM" 

21290 IF Syri<4 THEN PRINT "SR 1.4 ,2. 4” 

21290 Yad-0 

21230 IF K j ■ 1 THEN Yad*..9 
21310 IF Ic I *0 THEM 
21323 PRINT "PR" ,X ,Y+Yad , " M 
2 1 330 ELSE 

21340 PRINT “PR" ,X , Y + Yad , "PO* 

21350 ENO IF 

21360 IF Sym;3 THEM PRINT "SR 1.2, 1.5" 

21373 IF Sytt-4 THEN PRINT "*JC2 ,4 .39 ,0 ,-9 ,-4 ,0 ,0 ,9 ,4 ,0; • 

21390 IF Sym-5 THEN PRINT "UC3 ,0 .99 .-3 ,-S .-3 .5 ,3 t S ,3 ,-S; " 

21390 IF Sym-S THEM PRINT ’*UC0 .5.3 ,99 ,3 , -9 ,-S ,0 ,3 ,9 1 " 

21403 IF Syn-7 THEN PRINT "UC0 ,-S .3 ,99 ,-3 ,9 ,S ,3 , -3 ,-9 ; " 

21410 IF Kj-1 THEN PRINT "3M:PR 3, -.9" 

21420 SU8EN0 
21430 SUS Purg 
21440 9EEP 

21453 INPUT "ENTER FILE NRME TO 9E DELETED" ,F U sS 
21460 PURGE File* 

21 470 GOTO 21440 
21490 SUB ENG 
2U90 SUB Tdcn 
21530 COM /Cc/ C( ? ) , leal 
21513 DIM Enf( 1 ) 

21523 ORTR 3.10096351,25727.94369,-767345.9295,79025595.91 
21530 ORTR -9247495599. 5. 97E99E+1 1 ,-2.951925^13,3.940795+14 
21540 -HERO Cl ♦ ) 

21550 BEEP 

21560 INPUT "GIUE R NRME FOR FILE TO BE CRERTED" ,F i le* 

21570 9 EEP 

21580 INPUT "SELECT TUBE ( 0-UH . 1 -HF ,2-UT ) " . 1 1 t 
21590 9EEP 

21600 INPUT " SELECT THERMOCOUPLE TYPE ( 0-NEU , 1 -OLD > " , leal 
21E10 IF 1 1 1 * 2 THEN 0i*.3127 
21620 CRERTE BOAT File*,! 

21530 ASSIGN 0Fiie TO = :ieS 
21=40 OUTPUT 0F 1 1 e ; 1 1 1 
21553 J-0 
21560 9EEP 

2 1 973 1 INPUT "ENTER MONTH , GATE AMO TIME f MM r OC : HH : MM : SS > " .Da t 
2 1 590 ’ OUTPUT T G9 ; “ TO” ; Da * • $ 

2 1 593 1 OUTPUT 709; "TO" 



2 1 "’00 1 ENTER *09: DateS 



21*13 


PRINTER IS 


21*20 


0 O T M T 


21*30 


= ° I n r 


2 1 7 4 0 


°° I NT 



; n t h 



da *2 and tine: 
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2 1 "SO 


°9IMT USING " 


i 0.x , " " 


Fn 3 


Tin 




r«, Mm 2 


r ^ ^ 


2 l ""SO 


re u*sifj THEN 














2!? -, 3 


PRINTER IS 73 ! 












2 1 793 


C O I KJ J 














21793 


PRINT " 


Month, da t 


a and 


t in-: 


: " : OateS 




2 1 9GG 


t p 1 1 1 —3 THEN 


DO [NT 


USING " 


! OX , " " 


T:iba 


T' - - : 


U • a land Snoot h 


21313 


r c ru-,] THEM 


PRINT 


USING " 


l OX " 


Tube 


Type: 


High c Lu.*. 


2 1 923 


IF Itt*2 THEN 


PRINT 


USING " 


t jjy •• " 


Tube 


Tw^. • 
■ / ^ • 


Turbo-9 


21933 


PRINT 














2 1343 


PRINT USING " 


13X " 


Fna 


Tin 




Tev Uu " 2 


Tdrop 



C ! 953 PR INTER I S ! 

2 ! 959 K « » 

2 1 970 ENO IF 
21993 BEEP 

2 1 993 INPUT "ENTER FLOUMETER REDOING" .Fna 
21903 OUTPUT 705; "AR AF3 A L4 UR!" 

2 1 9 ! 3 FOR L-0 10 4 
21923 OUTPUT 739; "AS SA" 

21933 IF L>3 AND L<4 THEN 22310 
21943 S-0 

21953 FOR 1-3 TO 9 
21950 ENTER 709;E 
21973 S-S+E 
21993 NEXT I 

21993 IF 1-3 THEN Enf ( 0 >-ASS( S/ 13 > 

22303 IF L-4 THEN Enf < 1 >-A2S( 5/ 1 0 > 

22313 NEXT L 

22323 OUTPUT 739; "AR AF33 AL33 UP I " 

22330 OUTPUT 739; "AS SA" 

22340 Etp-3 

22353 FOR 1-3 TO 9 

22350 ENTER ?39;Et 

22373 Etp-Eip+Ei 

22390 NEXT I 

22390 Etp-Etp/10 

22133 T in-FMTvav< En f ( 1 > > 

22113 Tav-FNTvsvIEnf < 0 > > 

22123 Grad-37 . 9953+ . 104399*rin 

22133 Mdot-3 . 9S57E-3+Fns ♦(3.51 9S5E-3-Fma •< 9 . 9230SE-S-Fna ♦( 1 . 23S99E-7-Fna *4 , 3 1 997 
c_ » 3 ) ) ) 

22M0 Uw-Mdot/< l003*PI*Ot *2 )M 
22153 Tdrop-Etp*! .E+S/( 10-Grad ) 

221 S3 PRINT USING " 1 0X , 3' 00 .00 ,4X ) , 1 X , 2 .00 , 4X ,MZ . 40 " ; Fna , Tin , Tev ,Uw "2 , Tdrop 
22170 PEEP 

22190 INPUT "WANT TO ACCEPT THIS OATA SET? ( l -Y ,3-N ) " ,Qk 

22193 J-JM 

22233 IF Qk -3 THEN 

22210 J-J-l 

22223 GOTO 21993 . * 

22230 ELBE 

22240 OUTPUT 0FtIe;Fna ,Ersf( •) ,Etp 
22250 PRINTER IS 731 

22250 PRINT USING ’ 1 3X ,3( 00.00 ,4X > , 1 X ,2.00 , IX ,MZ . 10" :Fna .Tm.Tev ,Uw '2 , Tdrcp 
222' 7 3 PRINTER IS ! 

22293 9E”P 

22233 INPUT ''JILL THERE 3E ANOTHER OATA SET? ( 1 = X ,0-M ? '* ,Go_on 

22303 IF Gc THEN 2 1990 

22313 ENO IF 

22323 ASSIGN ?F : 1 - r 0 * 

22 333 -o INTER IB ""01 
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::mo print 

22350 PRINT USING “ » Qy '• '* M Q T E : ” ” II "" data ar- stored m file "" , ! SrV' jJ F 



IIZG2 PRINTER IS ! 



33373 SU2END 
33333 3US Ucprt 
33303 PRINTER IS ! 

22400 SEEP 

22410 INPUT "Enter Uo Fila Mane'.FiIaS 
33433 SEER 

22439 INPUT "Number* of Data Puna” ,Nrun 

22440 INPUT "Co You Want a Plot Fi la?< ! *Y ,0«N }” , Iplci 
33433 SEER 



ZZ4E0 IF I-'lot 5 *! then 

ZZ 470 INPUT "Give Plot Pile Nana" ,P_f • lei 
33433 CREATE SCAT P files 4 
22490 ASSIGN ©Plat TO P.flls* 

22S00 END IF 
22EI9 PRINTER IS 701 
22520 PRINT 
22530 PRINT 

22S40 PRINT USING - 1 0X , " - Uai&r Ual Uo 

225E0 ASSIGN ©Fils TO File© 

22550 IF Iplet-I THEN ASSIGN ©Fllal TO P_?llaS 
22570 FOR I-! TO Nrun 
22590 ENTER ©Filst'Ju ,’Jc 

225S0 IF Iplot-1 THEN OUTPUT ©Fiia!:Uu,Uo 
22S00 PRINT USING “ 1SX ,0.00 t 5X .MZ.30E* :Uw ,Uo 
225 !0 NEXT I 



22520 ASSIGN ©File TO • 

22520 ASSIGN ?Fila» TO * 

22540 PRINT USING " 1 0X , “ "NOTE : " " ,ZZ data sets ara stored in f i la " 
n ,F: 1 •« 

22SE0 IF Ip lot- 1 THEN 

22S60 PRINT USING " ! 0X , " “MOTS : " " , Z Z , * X-/ Paira ara atorad in fils * 

n,R_f Ua3 



22570 END IF 
225S0 PRINTER IS I 
22SS0 SUSENO 
22700 SUB Saiect. 
22710 COM / 1 dp / Ido 
22720 SEEP 
22730 PRINTER IS l 



22740 


PRINT 


USING "4X," 


•• 


Sal 


i act option:""" 


227S0 


PRINT 


USING "SX," 


- 


0 


Taking data or ra-procea a mg pravioua 


227S0 


PRINT 


USING "SX," 


" 


i 


D Icttmg data on Leg-Log """ 


22770 


PRINT 


USING " 5 X , 




-* 


Plotting data on Linear 


22790 


PRINT 


USING "SX." 


•• 


3 


Make oroaa-piot ccef p t file 


22^90 


PRINT 


USING "SX." 


•• 


4 




22900 


OR TMJ 


USING "SX." 


•• 


5 


°urga 


22910 


PRINT 


USING "SX." 


- 


S 


T-Orcp oorract ion" " " 


22920 


PRINT 


USING "SX." 


“ 


“7 


Print Uo c 1 1 a 


229Z0 


oo r 


USING "SX " 


” 


Q 


Modify X-Y ft la""" 


22240 


po i m r 


USING "EX ’’ 


■* 


Q 


Neva 


22250 


PRINT 


USING "SX " 




19 


Conb/F i/uc 


22250 


IMP U T 


r 









data 



22370 

2^320 

22920 

22900 



rc t j*. -«3 rucf| CALI 
IF I do-' THEN CALL 
rc I-jo-Z r HEN CALL 
IF Id" ”3 r HEN CALL 



Ma i r. 
° lo t 
P L i n 



, I S A ” ; Nru 
, l 5A” ; Nru 
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H Z *3 ! *2 IF Id-^i THEM CALL M^ir* 

22922 IF I dp -5 f u c n rai_i_ p-jrg 
32330 IF I dp*E THEN CALL Tdcn 
ZC940 IF Idc a T THEN CALL 'Jeer* 

229E3 IF I dp =9 THEM CALL X,wod 
22950 IF I dp-9 THEM CALL Ms vs 
22973 IF I dp -13 THEM CALL Conb 
22390 SUE END 
22 990 3»jB X y r c d 
23030 PRINTER IS 1 
23310 SEEP 

23023 INPUT “ENTER FILE NAME" .FiieS 
22020 ASSIGN ®Fil-l TO F;i-S 
23340 9EEP 

23353 INPUT “ENTER NUMBER OF X-Y PAIRS " ,Np 
23053 9EEP 

22373 INPUT “ENTER NEU FILE NAME " , F i 1 523 
22383 CREATE BOAT Fi 1-23,5 
23090 ASSIGN 0File2 TO F;le2S 
23100 BEEP 

23110 INPUT “ENTER NUMBER OF X-Y PAIRS TO BE DELETEO* ,Ndel 
23120 IF Ndei-0 THEM 23153 
23130 FOR 1-1 TO Ndel 
23143 BEEP 

23153 INPUT “ENTER QATA SET NUMSER TO BE OSLETEO” ,Nd( I ) 

231 50 NEXT I 

23173 FOR J-t TO Np 

23190 ENTER ®Fil-1 :X # Y 

23193 FOR 1-1 TO Ndel 

23233 IF Md< I )-i THEN 23240 

23213 NEXT I 

23223 OUTPUT 3File2;X.Y 

23230 PRINT J ,X ,Y 

23243 NEXT J 

23253 PRINTER IS 701 

23250 ASSIGN SFilel TO ♦ 

23270 ASSIGN ®Fiie2 TO ♦ 

23290 SUBENO 
23293 SUB Move 
233001 FILE NAME: MOVE 
23310' 

23323 DIM 9op(56 > ,A(SS ) ,9(SS ) ,C(SS > .0(65 > ,E(5S ) ,F(SS ) ,G(SS > ,H< 56 > , J( 66 > ,K(SS ) ,L( 
56 ? ,M( 56 ? 

23330 DIM Toid3<56HI41 ,N( 66 > ,Ur< 65 ) , Ir( 55 ) 

23340 BEEP 

233S3 INPUT “OLQ FILE TO MOLE" ,02_f i I -S 

23350 ASSIGN ®F;ie2 TO D2_ f ; 1-S 

23370 ENTER OF; i-2 ;Nnjn ,Oate$ ,Ldtc 1 ,Ldtc2,Itt 

23390 FOR I- 1 TO Nr*jn 

23393 ENTER OF i le2 ; Boot I ) . To IdS' p 

23400 ENTER ®F;1-2;A< I J ,8< I) ,C( I ) ,0' I > , E < I > ,F' I ) ,G( I > ,H( I > , J ( t > ,K( I ) ,LC I > ,M( I > ,N 
C I ) 

22413 ENTER OF • 1 -2 ; Ur ' I > , lr* I ) 

22420 NEXT I 

22433 ASSIGN ®F;l-2 T 0 * 

22443 SEEP 

234S3 INPUT “SHIFT DISK AND HIT CONTINUE" ,01: 

93450 SEEP 

23 4~2 INPUT “INPUT BOAT 3125“, 3;e- 
23493 CREATE SCAT 02 * : L-S Sir- 
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23430 assign sf i i * i ro Q 2 _fii-$ 

235G0 OUTPUT 5F i I e ! : Nrun ,0a t *5 ,Ldt c 1 ,Ldts2 , 1 1 1 

23510 FOR r-l TO Nrur 

23520 OUTPUT OF 1 ! * ! ; 9cp< I ; TsidSI I > 

22530 OUTPUT ?Fi l- 1 : A< I J ,3' I > ,C( I ) ,0< I ) ,E< I ) t F' I ) ,S< I > .H< I ) , J' I ) ,K< I ) ,L< I ) ,M< I > 
N< ! ) 

23510 OUTPUT OF i Is I :UrU > , Ip' I ) 

22550 NEXT I 

23550 ASSIGN OFli-l TO ♦ 

225"*0 f RENAME "TEST" TO 02 ft!*S 

23590 SEEP 2000 ,.2 

27cag SEEP 4000 .2 

23503 SEEP 4030 [.2 

22510 PRINT -DATA FILE MOVED* 

23523 SUSENO 
23530 SUB Comb 
23540* FILE NAME: C0M9 
22550* 

23550 DIM EmM 12) 

23673 .SEEP 

23680 INPUT *0L0 FILE TO FIXUP" ,02_f l Is* 

23693 ASSIGN *Flls2 TO 02_fil*S 
23730 01_filsS-*TEST* 

23710 CREATE SOAT OJ.flieS.20 

23720 ASSIGN OFHel TO Ol.flleS 

23720 ENTER SFi le2 :Nrun .OateS .Ldtc 1 ,Ldtc2 . It t 

23740 Mpunm-20 

23750 IF K-0 THEN OUTPUT OFils! jNrunm ,0at-3 .Ldtcl ,Ldts2 .1 tt 
23753 FOR 1-1 TO Nrun 

23773 ENTER OF 1 1-2 ; Sop , To I dS ,Enf ' ♦ > ,Vr . Ir 
23790 OUTPUT SFl 1*1 ?8op , ToldS .Emf ( • > ,Vp , Ip 
22790 NEXT I 

23930 ASSIGN *Fli*2 TO • 

22913! RENAME "TEST* TO D2_flle« 

23823 BEEP 2300.. 2 
23930 SEEP 4000 '..2 
23940 SEEP 4300,. 2 
23850 8EEP 

23953 INPUT "UANT TO AOO ANOTHER FILE ( 1 -Y ,3-N )?* ,Oka 
23873 IF Oka-1 THEN 
22993 K-1 
22990 SEEP 

23900 INPUT "GIVE NEU FILE NAME*. Nf lie* 

23910 ASSIGN SFiIe2 TO Nftl*9 
23923 GOTO 23730 
23920 END IF 

23940 ASSIGN OFll-2 TO • 

22950 SU8EN0 
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